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High-Temperature Hypersonic Flow’ 
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Summary 


Two important features of hypersonic flow over slender or thin 
bodies are the displacement effect of the boundary layer and the 
large downstream influence of leading-edge bluntness. The 
present paper contributes new theoretical and experimental re- 
sults on this problem. 

The interaction of the two effects is treated theoretically by ex- 
tending the basic shock-layer concept. In the outer inviscid 
flow, a model consisting of a detached shock layer and an entropy 
layer is introduced to account for bluntness. In the boundary 
layer, the approximate solution is found to be governed by a local 
flat-plate similarity. Under the assumption of a strong bow 
shock and a specific heat ratio close to unity, a theory is developed 
For flat-plate afterbodies, the theory 
yields a solution agreeing with blast-wave theory at one limit and 
strong-ipteraction theory at the other. Within the framework 
of the present theory, the problems involving angle of attack are 
also analyzed. Complementary to the above study, a hypersonic 
similitude involving strong shocks, but not requiring y close to 
The solutions obtained in the leading approxima- 


for an arbitrary thin body 


one, is given 
tion provide a natural comparison with experimental data cor- 
related on the basis of this similitude. 

Flat-plate experiments in air, conducted in the C.A.L. 11 X 
15-in. hypersonic shock tunnel under cold-wall conditions, in- 
cluded measurement of surface heat-transfer distributions and 
schlieren studies for zero and nonzero angle of attack. Steady 
laminar heat-transfer rates were measured by means of thin-film 
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resistance thermometers at air test-flow Mach numbers around 
12, free-stream Reynolds numbers from 1.4 XK 10% to 1.8 & 10° 
per in., and leading-edge Reynolds numbers from 3 to 1.5 & 104 
For most of the experiments, airflow stagnation temperatures 
ranged from 1,900°K to 2,000°K, giving wall-to-stagnation 
temperature ratios of about 0.15. The range of test conditions 
at this stagnation temperature encompassed the limiting cases of 
dominant bluntness and dominant viscous-interaction effects 
Heat-transfer distributions were also measured on a sharp plate 
for air stagnation temperatures ranging from 2,000°K up to 
$000°K 

The experimental data are quite well correlated in terms of the 
foregoing theoretical similitude variables characterizing com 
bined effects of boundary-layer displacement and bluntness 
The correlations obtained suggest that for the present experi 
mental conditions, at least, the hypersonic viscous similitude is 
valid even with leading-edge bluntness. Quantitative agreement 
with the zero-order theory, discussed in detail in the paper, 1s 


generally fair. 


Symbols 


speed of sound in the free stream 

constants 

(uw Ou /OyY)w/(1/2)p.. L*, the skin-friction coefficient 

(x OT/Oy),,/p.U(H, — H q/p.U(H, — H,), 
the surface heat-transfer coefficient 

specific heats of the gas at constant pressure and at 
constant volume, respectively 

the constant of proportionality in the linear vis- 
cosity-temperature relation u/p Cisse. 

w( Tx )/p(T.. )\(7 Tx ) 

the drag of the blunt leading-edge or nose 

the specific internal energy of the gas 

symbols representing certain functions of the vari- 
ables involved 

the specific enthalpy 

h + (1/2)(u? + v?), the total specific enthalpy 

vth-order Bessel’s functions of the first and second 
kind, respectively 

Dy/(1/2)p.. Ut, the leading-edge (nose) drag coef 
ficient 








L = 


m saad 


Subscripts 


b 
B 
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M*ke(t/x), a parameter controlling inviscid tip- 
bluntness effect 

a reference length 

a constant associated with the streamline in the en- 
tropy layer 

U/a,,, the free-stream Mach number 

the pressure 

Prandtl number pC,,/x 

local surface heat-transfer rate 

the gas constant in the equation of state p 

the Reynolds number p,, Ux/u, 

the thickness of the blunt leading-edge or nose 

the temperature 

To(1 + 3(T,,/7T ) |/6, a reference temperature for de- 
fining C 

the velocity components parallel to the x and y axes, 
respectively 

the value of u in the free stream 

u — U, the streamwise perturbation velocity 

rectangular coordinates in directions parallel and 
normal, respectively, to the free-stream velocity, 


with the origin at the nose 


the y ordinates of the inner edge of the entropy 
layer, the outer edge of the entropy layer, the 
shock, and the body surface, respectively 

the variables related to Y, and x, respectively; the 
relations between z and Y,, and between ¢ and x, 
are different for different cases 

the angle of attack of the flat plate; positive a refers 
to compression or wedge flow over the surface in 
question 

the variable x.K, 

C,/C,, the specific-heat ratio 

K.Ma/2x2, a parameter governing the bluntness- 
displacement-incidence interactive effect 

and displacement 


the boundary-layer thickness 
thickness, respectively 
(y¥ — 1)/(v +1), (vo — 1)/(v¥a + 1), respectively 
a typical flow-deflection angle which may be taken 
as dY,/dx 
(H — H,)/(H. — H,,) 
a function whose derivative 0®/0n is u/U 
thermal conductivity of the gas 
the viscosity of the gas 
certain variables 
two spatial variables used in the boundary-layer 
analysis, see Eq. (2.1) 
the gas density 
a body-thickness parameter used in the similitude 
consideration 
M3\/C/Re,,M8V/C/Rez, 
(0.664 + 1.73(H,,/H,,)|x, a parameter governing 
the boundary-layer displacement effect 
the constant exponent in the relation ~ « 7@ 


pertaining to the inner edge of the entropy layer 

pertaining to the Blasius solution 

at the outer edge of the entropy layer 

based on the reference length (x = L) 

pertaining to the condition behind the strongest 
part of the shock at the nose 

behind the shock (at x) 

based on the nose thickness (x = /) 

at the body surface or wall 

at the outer edge of the boundary layer 

pertaining to the free-stream stagnation condition 
(for the temperature) 

pertaining to the free-stream condition 

based on the reference temperature 7 
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(1) INTRODUCTION 


i bw EFFECTS OF leading-edge or nose bluntness and 
boundary-layer displacement are important in the 
hypersonic flight of thin bodies at high altitude. In 
recent years, these phenomena have received consid- 
erable attention in both theoretical and experimental 
studies. Theoretical studies to date have considered 
only limiting cases where either one or the other effect 
Thus, on one hand we have the so-called 
leading-edge 


is dominant. 
blast-wave theory,'~* for dominant 
bluntness eTect, which has been criticized on various 
grounds.*' In many cases the blast-wave phenomena 
appear to be subject to the influence of viscous effects 
in an obscure way, tending to discredit the theory. On 
the other hand, for dominant boundary-layer displace- 
ment effect we have the rather well-established theories 
of strong and weak interactions for a flat plate, as well 
as treatments of related problems. '!~°6 

Experimentally, studies of hypersonic flow over flat 
plates with bluntness'®*—* have been carried out 
mainly for adiabatic-wall conditions, in which case 
boundary-layer displacement effects are considerably 
larger than expected for the cold-wall conditions of 
hypersonic flight. Most of the experimental investi- 
gations have been concerned with measurement of 
shock-wave shapes and or surface-pressure distribution 
(the studies of references 51, 53, 55, and 57 are excep- 
tions). Test Mach numbers in helium have ranged 
from about 10 to 24, and in air from about 5 to 9.6. 
Generally speaking, the experiments have provided use- 
ful insight into leading-edge bluntness and boundary- 
layer displacement phenomena and have served to 
assess the limitations of theory. 

The present paper gives the results of a theoretical 
and experimental study of leading-edge bluntness and 
boundary-layer displacement eifects for the general 
situation in hypersonic thin-body flow where both 
phenomena may be equally important. In Seetion 2, 
a theoretical approach to the general problem is de- 
veloped using the concept of a detached shock layer for 
the inviscid flow and local similarity for the boundary 
layer. The leading approximation for the flow model 
adopted is developed in detail; its simplicity permits 
solutions to numerous cases of interest. The solutions 
obtained reveal the combined effects of bluntness, vis- 
cous interaction, and body angle in hypersonic flows over 
thin bodies. In addition, the classical hypersonic 
similitude of Tsien and Hayes is extended to include 
bluntness and displacement effects. 

On the experimental side, Section 3 describes studies 
in the C.A.L. 11 X 15-in. shock tunnel of hypersonic 
airflows at Mach numbers around 12 over sharp and 


blunt leading-edge flat plates. For most of these 
studies, the stagnation temperature was close to 


2,000°K. The corresponding ratio of wall-to-stagna- 
tion temperature was about 0.15, giving highly cooled, 
laminar boundary layers. Shock-wave shapes were 
determined by schlieren photography, and detailed, 
surface-heat-transfer distributions were measured by 
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BOUNDARY-LAYER 


thin-film thermometry. Leading-edge Reynolds num- 
bers were varied from 3 to 1.5 X 104 for stream Rey- 
nolds numbers from 1.4 X 104 to 1.8 X 10° per inch. 
The results at zero angle of attack encompass the 
limiting cases of dominant bluntness and dominant 
viscous-interaction eifects. Angle-of-attack effects 
were studied at low and high leading-edge Reynolds 
numbers. Heat transfer was also measured on a sharp 
plate for stagnation temperatures up to 4,000°K. 

In Section 4, correlations of the experimental data 
are made in terms of the extended viscous similitude 
for bluntness, and the degree of correlation obtained is 
discussed. In addition, the experimental data are 
compared with predictions from the theoretical solu- 


tions devel»ped in Section 2. 


(2) THEORY* 


(2.1) Scope and Assumptions 


In this study, a continuum flow model consisting of 
a laminar boundary layer adjacent to the body and an 
inviscid outer flow field will be assumed. Since the 
Mach numbers considered are high, and the bodies 
thin, the lateral extent of the flow field bounded by the 
bow shock is small. Our analysis may therefore be 
carried out within the scope of the boundary-layer 
theory'*° and the (inviscid) hypersonic small-disturb- 
ance theory.**—** Within this framework, two comple- 
mentary studies are made. First, and the major part of 
our analysis, is the development of an approach which in- 
cludes the interaction of boundary-layer displacement, 
tip bluntness, as well as body thickness effects. To 
achieve this end, certain approximations and simplifica- 
tions will be introduced. The second study is con- 
cerned with hypersonic similitude when both boundary- 
layer and tip-bluntness eifects are present, and with 
simplification of the similarity law for flows with strong 
shock waves. For this similitude consideration, a nuin- 
ber of the assumptions introduced in the first study are 
not required. 

Underlying the present approach are the two basic 
concepts of local boundary-layer similarity and a de- 
tached shock layer. Both of these concepts are de- 
sirable for a development consistent with the so-called 
Newtonian shock-layer approximation.**—** The con- 
sistency between local similarity and the shock-layer 
approximation has, in fact, been discussed earlier." 
The concept of a detached shock layer, which includes 
the so-called ‘‘free layers’’*:*? as special cases, is a neces- 
sary modification of the shock-layer theory to allow 
treatment of tip bluntness. 

A similar treatment of the tip-bluntness effect has been 
given previously by Chernyi.*> The essential idea of 
a detached shock layer is implicit in his integral method 
for hypersonic inviscid flows. While Chernyi’s method 
and the present analysis are equivalent in the idealized 


* Most of the detailed analysis in this section was done by Mrs. 
Angela Chang of the Aerodynamic Research Dept. 
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Fic. 2.1.(a) Top: Illustrations of the stream-line pattern of the 
flow field around a slender body with small blunt nose. (b) 
Bottom: Division of the inviscid-flow region about a slender after 
body into a shock layer and an entropy layer. 
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Fic. 2.2. Determination of the maximum perturbation vel- 

ocity | Au Us — u over a flat-plate afterbody in an inviscid 


hypersonic flow 


limit y — 1, they differ, however, in the higher-order 
approximation for y not close to unity. It may also 
be pointed out that the inviscid tip-bluntness effect has 
been considered recently by Guiraud.’* The idea of a 
detached shock layer is employed, and certain equations 
obtained are very similar to ours.* Guiraud’s study does 
not lead to explicit solution of the problem, however. 
Throughout the present analysis an ideal gas with 
constant specific heats will be assumed ahead of and 
behind the shock, although y may change across the 


* We thank Mr. Guiraud for calling our attention to reference 


66 in a recent private communication. 
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shock. In order to arrive at a simple form of the local 
similarity for the boundary layer, the Prandtl number 
is taken to be unity, the body temperature is assumed 
to be uniform, and a linear viscosity-temperature rela- 
tion is adopted. The analysis presented in this paper 
will deal mainly with steady, plane flows. 

The three requirements which control the accuracy 
of the analysis may be expressed as follows: 

(a) The disturbance velocity Au/U = [(u/U) — 1], 
as well as the square of the flow angle, in the inviscid 
flow field must be negligibly small in comparison to 
unity. 

(b) The bow shock wave must be sufficiently strong. 

(c) The specific heat ratio y must be sufficiently close 
to unity. 

The first requirement implies exclusion of the blunt- 
nose region itself from the present enquiry, and we will 
deal only with its downstream effect. The consistency 
of these requirements with the degree of approxima- 
tions involved will be examined later. The zero-order 
theory, which is to be studied in detail, pertains to the 
leading approximation under the last two requirements 
cited, and may therefore be compared with the classical 
(Newtonian-plus-centrifugal-force) result of Busemann 


for slender bodies.”?>~** 


(2.2) An Approach to the Problem of Tip Blunt- 
ness and Boundary-Layer Displacement 


In the following, the local similarity and its consist- 
ency with the shock-layer approximation will be dis- 
Then the theory of a detached shock layer 
Finally, the zero-order theory for 


cussed first. 
will be considered. 
the complete problem wil) be developed and applied 
to particular cases of interest. 


The Local Similarity and the Shock-Layer Theory 


The idea of assuming a local similarity as a means of 
estimating heat-transfer rates on blunt bodies at high 
speed was first proposed by Lees” and later utilized by 
a number of workers.*** Lees observed that the 
thermal boundary layer is rather insensitive to pressure 
gradient, particularly if the boundry layer is cooled. 
For most engineering applications it suffices to ignore 
the pressure-gradient term in the equation of motion, 
in which case the problem is reducible to that of a flat 
plate.* For the latter problem, a self-similar solution 
exists if the wall-to-stagnation temperature ratio is 
constant along the plate. 

It is essential to recognize that in the present study, 
where only thin bodies are considered, the direct 
pressure-gradient effect is unimportant not primarily 
because of a cold wall, but because the direct effect of 
pressure gradient on the hypersonic boundary layer 
under consideration is at most of the order 


e= (7 — 1)/(¥ +1) 


* The similarity under consideration may be more appropri- 
ately termed the “local flat-plate similarity.” 
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The condition of small ¢€ is a requirement in existing 
shock-layer theories for inviscid hypersonic flows.*~% 
It is natural to treat the present boundary-layer prob- 
lem by an analogous perturbation procedure which uses 


€ as a controlling parameter. The argument used for 


the local-similarity method in blunt-body applications 
is avoided in the present study. 

The foregoing considerations become evident if the 
variables x,y in the boundary layer are transformed as 


follows: 


sag f C p dx 
0 bd. L 
VR (f° p \( C p -) wits 
= teL 
” : 0 Pp L 0 P« io 


where Re, = 9.UL/ ta 
Assuming a linear viscosity-temperature relation 


ge, = C.1/7. 


and a unit Prandtl number, the system of differential 


equations is reduced to 


9 in = = 
29 inn + $¢,,, 2&|®,P:, #,f,,| a 


dp (* ge 
2 dx 2 
[( hh pax)/o* | 7 + 
ae 
1—- —je-@,?| | 
T y 


20,, + f0, — 24 ». — 0, + 


d 
’,0 in( 1 -7)| = (0 
dé To 


and the boundary conditions to 


6’ = ¢%,=0=0 at n=0 | (2.1: 

6 =0=1 at n = wf hte 
where 

’, = u/U 6 = (H — H,)/(H. — Ay) 


In arriving at these equations, quantities of the order 
of 1/eM/? and Au;/U higher than the remaining terms 
have been neglected. From the above formulation, it 
is quite clear that the leading approximation solution for 
small ¢ will be independent of &, and is therefore gov- 
erned by the Blasius equation 


2b 5!” + Dyby"” = 0 Loa) 
Pp = ,/(0) = Og’(~) —-1 = 0 fF ie 
with Og = ,’. The problem of determining the 


higher-order approximation is equivalent to finding a 
correction to the classical flat-plate solution for a slight 
nonuniformity of the outer inviscid flow. The non- 
homogeneous term appearing on the right in the first 
of Eqs. (2.1) reveals clearly the role of the wall-to- 
stagnation temperature ratio in controlling the degree 
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BOUNDARY-LAYER 


of approximation of local boundary-layer similarity. 
Changing 7,,/7,) from unity to zero may typically re- 
duce the correction term by 75 percent. 

Neglect of the pressure-gradient term in the bound- 
ary-layer momentum equation has also been employed 
by Lees'® to simplify the problem of strong-shock 
boundary-layer interaction. For this problem, how- 
ever, a self-similar solution exists even with the pres- 
sure gradient retained. The theoretical justification 
for neglecting the pressure gradient in cases in- 
volving arbitrary pressure distribution on the basis of 
the shock-layer theories has been pointed out only 
quite recently." 

In the leading approximation for small ¢, the surface 
heat-transfer rate, the skin friction, and the displace- 
ment thickness of the boundary layer are determined 
explicitly from the pressure distribution. For plane 


flows, 


) * » dx\—" 
MCy ~ 0.332%, (! Y( ‘ =) 
Pu J0 Pa 4 


Cy ~ 2Cun 


Us (0.064 + 1.73 = > xX (2.3) 
d ] ~ € 00 40 T x 


( x pb a) ( an 
J, po L p 


where x, = M*VC/Rez 


These results, in view of the approximations involved, 
will be subject to errors of order « and Au;/U. The 
usual errors inherent in the boundary-layer equations 
as well as errors from possible slip, temperature-jump, 
and vorticity-interaction effects, will presumably be 
no larger than those incurred by the approximation 
us ~ U9. 1.2326 The magnitude of Aus;/U will be ex- 
amined later. 

Before proceeding to the analysis of the inviscid flow 
region, the choice of a reference temperature for the 
determination of the constant in the linear viscosity- 
temperature relation will be discussed. The assump- 
tion of such a viscosity law may seem to be unrealistic 
since the temperature variation within the hypersonic 
boundary layer is quite large. The approximation may 
be improved considerably by choosing the constant C to 
represent correctly the value of viscosity at an appro- 
priate reference temperature. Instead of choosing the 
wall temperature as a reference, as was proposed earlier 
by Chapman and Rubesin“” and adopted later by others 
for hypersonic applications, we recommend, for the 
present study involving boundary layer over slender 
bodies, a reference temperature 7, defined by 

Ts = T.[1 + 3(Tw/To)]/6 (2.4) 
C= Ce = [u(Tx) u(T..) (T/T) 
This choice of 7% is particularly meaningful when local 
flat-plate similarity holds, as in the present problem 
For then we may apply Crocco’s integral relation, and 
the reference temperature 7} given above is seen to be 


whence 
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simply an average temperature across the hypersonic 
boundary layer over a thin body, i.e., 


"1 du "1 =6du 
Tx = 1 = 1 

JV u JV Us 
The skin friction and the heat-transfer rate on a flat 
plate under uniform pressure based on 7% of Eq. (2.4) 
agree exceedingly well with results based on the Young 
formula® at the higher Mach number range, provided 
the Prandtl number is not far from unity 


The Bluntness Effect and Detached Shock Layer 


The necessity for allowing the shock layer to detach 
from the afterbody of a blunt leading edge is evidenced 
by the breakdown of the standard shock-layer theory 
around the shoulder region where the pressure would 
accordingly fall below zero and the shock-layer thick- 
ness would become infinite. The essence of our analysis 
is that the pressure at the base of the detached layer 
will not be put equal to zero, a priori, as was done for a 
“free layer,’’*? but is to be determined as part of the 
solution to the problem. 

It is essential to realize that, within the inviscid, 
hypersonic-flow region bounded by the shock and the 
body surface or the outer edge of the boundary layer 
downstream of a small blunt nose, there is an inner core 
of considerable extent containing gas particles which 
have come through the comparatively forward, and 
hence stronger, portion of the shock (see Fig. 2.1). It 
is appropriate to term this region the “entropy layer,” 
since its specific entropy is much higher than that pre- 
vailing near the shock above. The temperature within 
this layer is accordingly much higher, and the density 
much lower, than the corresponding values near the 
shock; the gas within the entropy layer has been heated, 
so to speak, by the strong forward portion of the curved 
shock. It is natural to consider the remaining, and 
relatively thin, high-density part of the flow field behind 
the shock as a shock layer in the usual sense, except for 
the fact that it is no longer attached to a rigid body. 
The equivalent body surface for this shock layer can be 
taken, over the relatively slender, downstream portion 
of the flow field, to coincide with the outer edge of the 
entropy layer, y = Y,(x), which, in view of the particle- 
isentropic condition, must also be a streamline [see 
Fig. 2.1(b)] Across this streamline boundary, the pres- 
sure must be continuous 

The flow structure within a hypersonic shock layer 
has been treated quite extensively by Freeman” and 
others.*!—** Our present inquiry will deal with the 
entropy layer, the outer edge of which is the base of the 
shock layer. To investigate this, we need a knowledge 
of the pressure at the base of the shock layer in terms of 
Y,. The leading approximation for this relation, ac- 
cording to shock-layer analysis for a plane flow, is 


simply 


Pe = p.~U(Y,'2 + Y.Y,"’) (2.5) 
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which represents the Newtonian-plus-centrifugal (Buse- 
mann) pressure formula simplified for thin bodies. 

In view of the high temperatures attained within the 
entropy layer, the local velocity u could be reduced con- 
siderably from its free-stream value LU’. Since the as- 
sumption of «4; ~ U has been made in the boundary- 
layer analysis, and will be used later in the treatment of 
the entropy layer, an assessment of the magnitude of 
Au/U is important and will be made prior to further de- 
velopment. 

For the purpose of examining the magnitude of Au/U 
in the entropy layer, we consider first the particular case 
of a flat plate under the influence of the strong blast- 
wave effect produced by its blunt tip. For this case, 
the surface pressure distribution is known’ and we 
can determine the maximum value of the perturbation 
velocity (for arbitrary x) from the particle-isentropic 
condition and the Bernoulli relation. Maximum 
values of |Au/U| for y = 1.20, 1.40, and 1.667 are pre- 
sented in Fig. 2.2 as a function of x/kt. The results in- 
dicate that sufficiently far downstream of a blunt nose 
(x/kt 2 10), Au/U remains sufficiently small to permit 
the approximation u ~ U. 

Since in most cases 


Po < ps = [2¥/(y + 1)]M?e" 


where @ is the local shock angle (which gives also the 
typical local flow angle), from the entropy behind the 
normal shock and by the Bernoulli equation, a con- 
servative estimate for the perturbation velocity may in 
the general case be taken as 


JAu/U| < (1/2)@2-2/ 


where the exponent 2(y — 1)/y has the values 1/3, 
4/7, and 4/5 for y = 1.20, 1.40, and 1.667, respec- 
tively. 


We may now state in more explicit terms the re- 
quirements which will control the degree of accuracy 
of our analysis: 

(a) (1/2)(0)?"~/7 must be negligible in comparison 
with unity. 

(b) 1/6? must be sufficiently small in comparison 
with unity. 

(c) ¢ must be sufficiently small in comparison with 
unity. 

Obviously, condition (a) is essential for the assump- 
tion wu = U at the outer edge of the boundary layer and 
within the entropy layer as well. Condition (5) calls 
for a sufficiently strong shock, and (c) and (b) together 
will permit treatment by the shock-layer concept. By 
these requirements, the admissible ranges of the shock 
(and flow) angle 0, and «¢, are necessarily restricted: To 
satisfy the small-perturbation assumption, @ must be 
small, but not so small that the shock cannot remain 
strong. In addition, the shock-layer approximation 
requires ¢ to be close to unity, but not so close that the 
small-perturbation assumption will be violated. The 
ranges for 6 and e can, of course, be enlarged by increas- 
ing M and by including (consistently) terms of higher 
orders in e€ and 1///°6?. 


Our analysis will be based on the fact that the density 
within the entropy layer is low compared to that in the 
free stream, 1.e., p < p It is essential to examine 
the compatibility of this fact with the requirements set 
forth above. Since the density depends not only on 
the pressure but also on the entropy, we should distin- 
guish two situations. In one, the streamlines over most 
of, or at least a considerable part of, the entropy layer 
come from the immediate vicinity of the blunt nose. 
In other situations, the streamlines from the blunt- 
nose region reach only a small fraction of the entropy 
layer. The rest of the streamlines come from the 
slender but still the more forward and stronger portion 
of the curved shock. In the first situation, the specific 
entropy is at the same level as that prevailing at the 
nose, and the order of magnitude of the density may 
therefore be represented for particle-isentropic flow by 


y+1/p\'% = 
Este = (*) = of (6)? | 
p y¥— 1 \dn € 


Comparing this with requirement (a@) which implies 
that 6°’” < ¢ for e < 1, we see that the local density is 
indeed small in comparison with the free-stream value 
when (a) is satisfied. In the second situation, the 
density distribution given by the strong blast-wave 
theory'~*:*!:*? should be valid, since streamlines over 
most parts of the entropy layer now come from the 
slender portion of the shock. The fact that the blast- 
wave solutions reveal the structure of a detached shock 
layer and a low-density entropy layer is well known, and 
the ratio p/p. in this situation may, in fact, be repre- 
sented (assuming a flat-plate afterbody) by 


p/p» = O[(1/ee*] 


which is small for all practical purposes. Presumably, 
the second situation may be realized only for an ex- 
tremely small blunt nose (hence large «/t) together with 
an extremely high Mach number. In either situation, 


the supposition 
p/p. <1 (2.6) 


within the entropy layer can be well justified. Its 
error is presumably no larger than that incurred by the 
assumption u ~ U, that is, requirement (a). 

An immediate consequence of the low density in the 
entropy layer is that the variation of pressure across it 
is small and negligible. We can examine this from the 
equation of motion governing the momentum compo- 
nent in the direction transverse to the layer. The frac- 
tional variation in pressure across the entropy layer is 
readily seen to be 


—- = 0[ 2 (45 *)|-0(4) «1 (2.7) 
Pe p Y, p 


and is therefore negligible according to our require- 
ments. 
In view of the above conclusion, we can now employ 


the simplification 


p ~ p(x) (2.8) 


uc U 
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for study of the flow pattern within the entropy layer. 
The equations stating continuity and particle-isentropy 


( a) Pm 2) ) v 
u v == 
Ox Oy , , Oy 


( re) 4 > ) | l ( fe) 4 > ) ; 
1 v np= .— v 1p 
Ox Oy , Y Ox Oy ae 


Eliminating the density p and making use of the sim- 
plification on uw and p, these two equations give an 
equation for v, orv/u. This equation can be integrated 


are 


to yield 


a V,"(x) + (Y, — y) d In p. (2.9) 
u ¥ dx 
for which the boundary condition at the inner edge of 
the entropy layer (also the outer edge of the boundary 
layer), ie.,v/u = Y,’(x) at y = Y,(x), has been satis- 
fied. Putting dy/dx for v/u, we can integrate Eq. (2.9) 
once again and obtain an equation for the streamline 


[lv — Vi(x) | p.(x) = m (2.10) 


where m is a constant differing for each streamline. 
This equation should hold throughout the entropy 
layer. Applying it at the streamline boundary y = 


Y,(x), we arrive at 
[Y.(x) — Y,(x)]" p.(x) = constant (2.11) 


This simple integral relation between the pressure and 
the cross-sectional area of the entropy layer is basic 
to our subsequent development. 

The pressure /, is that at the base of the shock layer, 
Y = Y,, and is given by shock-layer theory once 
Y,(x) is specified. This relation for p, from shock- 
layer theory can be used to eliminate p, in Eq. (2.11), 
and we can write down, at least in principle, an ordi- 
nary differential equation involving only Y,. The 
shock-layer theory gives the Busemann formula, Eq. 
(2.5), as a leading approximation. The pressure-area 
relation consistent with this approximation is therefore 
[from Eq. (2.11)] reducible to 


(Y, — Y,)(Y.Y.’)’ = constant 


In carrying out a higher order analysis (in ¢ and 
1/6"), we should not only employ a pressure formula 
correct to the order desired, but also work consistently 
with a more accurate pressure-area relation based on 
Eq. (2.11). The first-order ¢ correction to the Buse- 
mann formula has been given explicitly by Freeman in 
reference 30. 

The constant on the right-hand side of Eq. (2.11) 
is evidently related to the drag of the small blunt 
nose which produces the downstream blast-wave ef- 
fect. Consider the balance of momentum change in 
the streamwise direction, assuming a supersonic flow. 
With the help of the equation of continuity and making 
use of the requirement (a), i.e., « ~ U, we can write 
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Dy x sh ay p 
“ + { PoY,'dx = i] dy + 
= +0 x ey 


Y rik >) 
f - dy — f a 'e 
yn 20 y-1 


where Dy is the drag contributed by the pressure dis- 
tribution around the small blunt nose excluded from 
the range of integration for the integral on the left.* 
The second terms on both sides of Eq. (2.12), as well 
as the third term on the right, belong to the orders « 
and 1/M/°@* higher than the others, and can be dropped 
in the leading approximation (a more detailed discus- 


to 
to 


sion is given in reference 65). To the same degree of 
approximation, we can replace the location of the shock 
front by the base of the shock layer. We can now put 
p = p(x) under the integral sign of the first integral on 
the right, and reduce Eq. (2.12) to 


Dy/2 = (p./(y — DY. — Y%) 


that is, 


eD y/p.U? (2.13) 


This is the pressure-area relation pertaining to the 
zero-order theory previously given, but with the con- 
stant now determined.t| To this approximation it 
seems that consideration of the momentum-integral 
relation not only determines the constant of Eq. (2.11), 
but also furnishes an alternative, and apparently 
simpler, derivation of the pressure-area relation. For 
the higher approximation, however, we may find it 
convenient to use Eq. (2.12) strictly as a subsidiary 
condition. The constant in the pressure-area_ rela- 
tion, Eq. (2.11), can be determined by passing to the 
limit x —~ 0 in Eq. (2.12). 

We may note that, on the basis of Eq. (2.13), the 
pressure p,, i.e., (Y,Y,")’, can never become negative. 
This is to be expected since the shock layer in our treat- 
ment is no longer required to follow the body surface. 
However, Eq. (2.13), as well as Eq. (2.11), have not 
excluded the possibility of an attached shock layer, 
and also contain the free layer as a limiting solution. 
We can see this readily by examining Eq. (2.13) for 
the nearly limiting situation of an extremely small 
nose. The constant on the right of Eq. (2.13) in this 
case then reduces to a very small but positive number. 
Now, for p, to be finite, it is clear that the shock layer 
must remain essentially attached to the body, i.e., 
Y, ~ Y,. As soon as p, approaches zero, the entropy 
layer will, accordingly, begin to thicken and the shock 
layer will detach to form a free layer. 

* The nose drag Dy has been defined as the total drag of a 
symmetrical nose, namely, y = +Y¥,(x). In the case without 
top-bottom symmetry, D,y/2 should be taken as the contribution 
of one part of the surface which begins at the apex, or the stagna- 


tion point, of the nose. 

+ For the axisymmetric case, the relations corresponding to Eq. 
(2.11) and Eq. (2.13) are, respectively, [Y,2 — Y,?]"p, = con- 
stant, and [Y,2 — Y,?][Y.’2 + (Y.Y."/2)] = 2eDy/mp.. U’. 
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Our general approach to the inviscid leading-edge 
problem is also similar to a treatment given previously 
by Chernyi,‘® as pointed out earlier. In Chernyi’s 
treatment, which is essentially a momentum-integral 
method, use is made of u ~ U and p = p,(x). The 
basic equations employed, written here for steady plane 


flows, are 


AEROS 


PA 
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to the constant on the right of our Eq. (2.11), certain 
knowledge of the entropy distribution near the nose 
seems necessary, according to Guiraud’s analysis. 
Hence, this term cannot be determined until the com- 
plete problem is solved. In our analysis, however, the 
constant is found with the momentum-integral relation, 
Eq. (2.12), as an auxiliary condition, and is thus de- 


termined by the nose drag (further discussions are 


Dn + [. prY,'dx = Po 1. m- 4 given in reference 65). 
2 J +0 y= 4 In what follows, we shall apply the results of the 
‘ (2.14) above analysis to the problem involving both tip 
U?Y,(Y,')? bluntness and boundary-layer displacement effects. 


Pr. = poU*(Y,’* + YY," 


It is evident that, as far as the leading approximation 
for small ¢ and 1/.1/°6" is concerned, the above equation 
and Eq. (2.13) obtained by us are essentially equiva- 


The zero-order theory for small ¢ and large .1/°6° will 
be developed in detail. 
The Zero-Order Theory 


For small ¢, the boundary layer has been shown to 
be governed by a local similarity; assuming, in addi- 






























































lent. When the higher order approximations are con- tion, that the shock wave involved is strong, the simple Frc 
sidered, we must compare Chernyi’s equation, 1.e., pressure-area relation Eq. (2.13) for the entropy layer 
Eq. (2.14), with our Eq. (2.11). We may note, in should also hold. Combining these results, the leading- 
Eq. (2.14), that implicitly the outer edge of the entropy edge problem of hypersonic flows subject to the tip- i 
layer has been taken to be the shock itself, and the bluntness and boundary-layer displacement effects can J 
= : 3 - - . u 
Busemann pressure formula has been used. Hence, it be treated in a rather straightforward fashion. Eq 
is doubtful that Eq. (2.14) would yield a valid correc- To apply Eq. (2.13) to this problem, the inner edge ot 
tion to the leading approximation. At any rate, the of the entropy layer, Y,(x), will be taken as an equiva- , 
two treatments will be seen to be quite different when lent body surface which should account for the dis- the 
higher order approximations are considered. placement effect of the hypersonic boundary layer. 
The entropy layer has also been studied recently by Inasmuch as the assumption u; ~ UL’ is valid,'*'® the [J 
Guiraud in references 7 and 66 in which an equation outer edge of the boundary layer may be represented 
quite similar to our pressure-area relation, Eq. (2.11), approximately by the displacement thickness 6x be- 
has been given. To determine the term corresponding cause of the extreme low density resulting from viscous 
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angle of attack. 


heating. In this instance, we can replace Y, by 
Y,.(x) + 6%(x). By local similarity, 6, is given by 
Eq. (2.3) explicitly in terms of the pressure distribution, 
which is now provided by the Busemann formula. 

The equation governing our zero-order theory is 
therefore 


[Y, — VY, — be |1(Y.Y.’)’ = ekt/2 
TAs ry? (2.15) 
€ (0.664 + 1.73 pes VY, 
To Mm (9777 


jy = 


After solving Eq. (2.15) for Y,, we can determine the 
pressure from the Busemann formula and the heat- 
transfer rate as well as the skin-friction from Eq. (2.3). 
The rudimentary manner in which the three flow re- 
gions—boundary layer, entropy layer, and shock 
layer—interact is contained in the ordinary differential 
equation given above. In the absence of tip bluntness, 
we have only the shock layer and the boundary layer. 
The first of the above equations then becomes 


oe _ if + bx 
and Eqs. (2.15) reduce to 


2(Y.— Ye)[V Y.¥Y.']' = 


€ (0.664 ‘225 7) Xi VL 


: (2.15a) 
To Mw 


The first of Eqs. (2.15) is singular at the origin and ad- 
mits the solution Y, ~~ x?/*. Thisis the correct behavior 
when the blast-wave effect of the leading edge is strong,”:* 
and confirms that near the leading edge the displace- 
ment effect of the boundary layer gives way to the 
more powerful effect of leading-edge bluntness. For 
cases which do not involve bluntness, Eq. (2.15a) 
reveals singular behavior of another type: as x —> 0, 
Y, ~ x*/4, This is the behavior to be expected when 
the displacement effect predominates.’ '®'* With 
these singular solutions as leading terms, solutions to 


Eq. (2.15) can be obtained by a forward integration 
from the origin. These singular leading terms in the 
asymptotic solutions for Y,, as well as the pressure and 
the heat transfer, will be examined subsequently for 
flat-plate afterbodies. 

Implicit in the above discussion is the assumption 
that the afterbody shape is sufficiently regular at the 
leading edge except for its abrupt beginning. It may 
also be pointed out that the blast-wave solution VY, ~ 
x*/8 is not the only mathematical solution to Eq. 
(2.15) near the origin for «ek # 0. For example, in the 
special case of Y,, = 6« = 0, a complete solution to Eq. 
(2.15) can readily be obtained. However, the only 
particular solution compatible with the small-distance 
requirement (a) turns out to be Y, « x*/* (this point 
is discussed further in Appendix I of reference 65).* 

As is obvious from the requirements (a)—(c), the 
error in the zero-order theory is 


O [e + 1/M%? + (1/2)(0)?7- 97) 


If terms of the next orders are consistently included in 
the analyses of the entropy layer and of the boundary 
layer, in the manner previously described, the first 
two error terms may presumably reduce to e* and 
1/M ‘64. 

In view of their frequent occurrence in the later de- 
velopment, it is convenient to introduce two new vari- 
ables 

K, = M*ek(t/x) ’ 


; / » (2.16) 

x. = €[0.664 + 1.73(7,/T») |M*V C/Re,S 
where Re, is a local Reynolds number based on the free- 
stream properties, namely, p..Ux/u Except for the 
factors « and e [0.664 + 1.73(7,,/T,)], these two vari- 


* The existence of a complete solution to Eq. (2.15) in the case 
of Y,, = 5« = O was called to the writer’s attention by Mr. J. P. 
Guiraud of ONERA, France, in a private communication 
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ables are those employed previously for the bluntness 
and the displacement problems.'~* 

In what follows, the zero-order theory will be applied 
to the problem of a flat plate in hypersonic flow subject 
to the three effects of bluntness, displacement, and 
angle of attack. In order to see clearly the manner in 
which these effects will act in combination, the simpler 
problems involving only two simultaneous effects at a 
time will first be singled out. 


The Flat Plate With Bluntness and Displacement 
For this problem, Eq. (2.15) can be reduced to 

s(ez’)’ — Ves’ = 1 (2.17) 

where the new variables are 

z= 8M(x2/K.)(Y./Ke), {= 16(y/K,*")° 

Note that both A,.v(= M’%*ekt) and x.2/K,( = e[0.664 + 

1.73(1,/T) |?M°C/kRe,) are independent of x, hence 

both z and ¢ are proportional, respectively, to Y, and 

x. Let @ = 2°/2 and \X = dd/dg. Eq. (2.17) may 

then be written 


V 26 (dd/dd) — VX = 1 


which can be integrated by a separation of variables to 
give 


$ = 4[Vd — (0/2) + (92/3) — In (1 + Vd) P 


This result provides the correct behavior at the origin, 
namely, Y, ~ x*/' asx — 0. To obtain a relation be- 
tween ¢ and ¢, we note that 


d¢é = (d¢/do)do = (1/d)(dd/dd)dr 


and the solution of Eq. (2.17) can finally be brought to 
the parametric form 


s = 2V2[Vd— (A/2) + AVD/3) — | 
In (1 + Vd)] | 
c= (1 + Vay — (22/901 + V8 + | 
911 + Vr)? — (46/3) + + (2.18) 
Vy) + (10/3) In (1 + Wd) = | 
4V In (1 + Vd) + | 
2[In (1 + Wa)]2 + (76/9) | 


This result for Y,, which may also be taken as the 
shock location Y,, is presented in Fig. 2.3(a) as 
M(x2/K)*(Y./K.) vs. (x./K,?/*). The corresponding 
surface pressure and the surface heat-transfer rate are 
given in Figs. 2.3(b) and 2.3(c) as (K,/x.?)?(Pu/YoPo) 
and ¢[0.664 + 1.73(T/To) (K./x2)°(M°C,,/V ya). 
We can see that Y, varies as x*/4as x > ©, and as x?/8 
asx—(. These two limits correspond to situations in 
which displacement and bluntness dominate, respec- 
tively, agreeing with strong-interaction theory in the 
first case and strong blast-wave theory in the second. 
One important feature here is that all parameters and 
spatial variables (including wall-to-stagnation temper- 
ature ratio and the specific-heat ratio) are grouped to- 
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gether to form one single variable x,/K,?/*.. An analo- 


gous variable, which controls the combined effect of 


bluntness and displacement, has been previously 
pointed out by Lees and Kubota in a discussion assum- 
ing a weak nose effect.”, From Fig. 2.3 we see that 
transition from the strong-blast-wave to the strong- 
interaction regimes occurs in the range 10~! < y,/K2 
< 2. 

The appearance of + and V in the nondimensional 
forms of pressure and heat-transfer rate introduced 
above result from the use of the dynamic pressure 
p.U? = y.p.M? in connection with the Busemann 
formula Eq. (2.5). Retaining the factor y 
replacing it by unity, is in keeping with the convention 
Furthermore, it 


, instead of 


of existing shock-layer theories.” 
permits e to represent conditions behind the shock 
even when e., is not small. 

It is of interest to examine more closely the asymp- 
totic solutions. They can be written, for large and small 
[x./K,?/*] respectively, as 
MY ,/x = 1.21X'/?,  1.31K,'/8 
Pu/ Yap» = (W3/2)x.,  (1/18)"8K 2/3 
«(0.664 + 1.73(T~/T») ](M°C,,/V yx.) = 

0.219X/2, 0.119K 2° 


(2.19) 


The appearance of x, in the above strong-interaction 
solutions indicates clearly the roles of « and the wall- 
temperature level. In view of the factor [0.664 + 
1.73(7,,/T,)| in x,, the boundary-layer displacement 
effect will depend strongly on the wall temperature. 
An interesting observation is that if we simply add the 
two asymptotic solutions for p,,/p.., the sum differs from 
the exact solution in the transition region by less than 
12 percent. This seems to lend support to the empirical 
linear-combination laws proposed previously by certain 
investigators, “”:>* at least when « is small. 

We may now compare the zero-order solution in the 
two limits, i.e., Eq. (2.19), with the corresponding self- 
similar solutions which do not assume y close to 
unity.*: 119.27 Some of these results are shown in Fig. 
2.3. We may conclude from these comparisons that as 
an approximation for the surface pressure and skin 
friction when y = 1.4, the zero-order solution may not 
be adequate, and a higher order solution (in ¢) appears 
desirable. For the surface heat-transfer rate as well 
as the shock shape, however, the zero-order theory 
works reasonably well, in spite of the fact that y in this 
case is actually not too close to unity. 


The Inclined Flat Plate with Displacement 


For this problem without bluntness, Eq. (2.15a) can 
be transformed to 
(¢ = ¢)(W22')’ = 1 (2.20) 
where the variables z and ¢ now stand for 
z= 4M'la 3V, ly 2x, ¢ = 4[M°a?/x,]? 


The angle of attack of the plate, a, is positive for com- 
pressive flow over the plate. The negative sign in 
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front of ¢ in Eq. (2.20) applies toa >0. In view of the 
strong-shock requirement, our solution for a < () per- 
tains only to the situation when the boundary layer is 
thick enough to induce a strong bow shock. This 
implies that the region of validity is confined to certain 
upstream portions of the flow field, and to a rather small 
range of negative a. Again, there is only one variable 
governing the problem. In this case, it is MJ*a*/x,. 
Eq. (2.20) appears to be quite different from that under- 
lying the analysis of Paiand Shen.'!".” The difference 
arises perhaps mainly from the different flow models 
assumed for the outer flow. 

For the case of a > (0 the solution of Eq. (2.20) for 
large M/*a?/x,, i.e., large x, will approach that for a 
wedge without displacement effect, i.e., 

Y,= art... (2.21) 

For small argument of the variable, the strong-interac- 

tion solution results, as given previously in Eq. (2.19), 
and 

2 x 3 


Y. = ys. 
: Xe wi 


) ax + o4* (2.22) 


This development near the sharp leading edge is valid 
for both a > 0 and a < 0, and reveals a pattern of 
asymmetry. The curves for a > (0 presented in Fig. 
2.4* give a smooth transition from the regime of strong 
interaction to that of a usual hypersonic wedge flow. 
Curves for a < () in the same figure show a gradual 
deviation from the Stewartson-Lees solution as we go 
downstream. For the shock shape and the pressure, 
the range of transition for a > () appears to occur in 
10-2 < M?a?/x, < 10 whereas for the heat-transfer rate, 
the value corresponding to the pure wedge-flow solution 
is closely approached as early as \/*a®/x, = 2. From 
the agreement of our zero-order solutions with the self- 


* The numerical results of Fig. 2.4 were obtained by digital- 
machine computation (IBM 704), using an interval of 100 points 


per unit of (¢)"*. 


tack. 


Results based on the zero-order theory of a flat plate under the combined effect of leading-edge bluntness and angle of at- 


similar solution for y = 1.40 in the two limits (see Fig. 
2.4), we may anticipate the zero-order results to be 
reasonable for heat-transfer, pressure, and shock shape, 
provided the surface temperature is substantially lower 
than the stagnation temperature. 


The Inclined Flat Plate With Tip Bluntness 


For the inviscid problem of the blunted flat plate at 
incidence, we have from Eq. (2.15) the simple differen- 


tial equation 


bo 


(zg ¥ ¢)(z2’)’ = 1 (2.23) 
where the variables z and ¢ now stand for 


z= 2a*!./ekt, ¢ = 2!a/*x/ekt 


In this case, we again use the same convention for a. 


Near the leading edge, we find from Eq. (2.23) the 
result 
Y,/t = (3/2)?/*(ek)/9(x/t)?/® + (1/6) @ (v/t) +... 
(2.24) 


which holds for both positive and negative a and reveals 
an asymmetrical perturbation of the blast-wave solu- 
tion. At the far-downstream end, where the effect of 
the surface incidence dominates, we will develop the 
asymptotic solution only for a > 0. Eq. (2.23) then 
yields for large ¢ 


iI 


c+14+ (1/Vo[As2V 9) + BY,.2V8)] + 


a] 


(2.25) 


¢+1+ (A/¢/*) cos (2V¢ + B) 


In view of the first and second terms we see that the 
shock will tend to the wedge surface, but never quite 
reach it, as we go downstream. More interesting, 
however, is the oscillatory behavior of the solution. The 
corresponding behavior of the pressure is 

VP, 2 —_ oo’)! — 
Pw/(Y¥oPoM*a*) = (22) = 
(2.26) 


1 — (A/¢#/4) cos (2V¢ + B) + 
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Comparing Eqs. (2.25) and (2.26), the oscillatory decay 
of the tip-bluntness effect will be more readily revealed 
in the pressure distribution than in the shock shape. 
Some results obtained from numerical integration of 
Eq. (2.23) by digital machine computation are pre- 
sented in Fig. 2.5.* The oscillatory decay of the tip- 
bluntness effect is indeed exhibited by the pressure dis- 
tribution for a > 0. 

The physical reality of this decay is doubtful. The 
result of the zero-order solution in this respect is evi- 
dently connected with the passage to the mathematical 
limit, y > 1. Apparently, the oscillation represents a 
repeated reflection or skipping of the re-attaching shock 
layer as it reaches the vicinity of the surface of a very 
slender afterbody at infinite Mach number. The 
phenomena may hence be interpreted as a form of 
glancing reimpingement of the shock layer. From the 
fact that this oscillation is a decay phenomenon, its de- 
tail will probably be substantially altered by higher- 
order € corrections. For y = 1.40, we may not antici- 
pate any oscillation in the pressure decay over a blunted 
wedge, since the amplitude of the fluctuation shown is 


never large. 


The General Case of the Flat-Plate Problem and Similitude 
Based on Zero-Order Theory 


For the flat-plate problem involving all three effects 
(x., K., and Ma) we may choose to use, for example, 
the transformation 


z= 8M(x2/K)*(Y./(Kx), £ = 16(x./K,?/*)$ 
and Eq. (2.15) becomes 
(g — Te)(z2")! — V2e’ = 1 (2.27) 
where IT is the parameter (independent of x) 
T = K.Ma/2x2 (2.28) 


Fig. 2.6 gives the results of integration for the surface 
heat-transfer rate. Shown are curves corresponding to 
five positive values of T, namely, Tl = 0, 10~*, 10—}, 1, 
and 10. The curve for Tl = 0 has been previously 
given in Fig. 2.3(c). From these results we can see the 
role of the wedge angle in controlling the combined 
effect of bluntness and boundary-layer displacement. 
The slight oscillatory mode revealed in the curve at the 
upper left corner of the figure for ! = 10 corresponds 
to the oscillatory pressure decay discussed previously 
for the inviscid wedge problem. 

It may be pointed out that the zero-order theory has 
been applied thus far to problems involving the rela- 
tively simple case of the flat-plate afterbody. Eq. 
(2.15), meanwhile, is applicable to an arbitrary slender 
afterbody shape. In any case, a solution can be ob- 
tained by a “forward integration.’”’ With Eq. (2.15), 
an interesting aspect of the problems of optimization 
from the viewpoint of lift, drag, or heat transfer, which 

* An interval of 100 points per unit of ¢1/3 was used in the com- 


putation. 
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take into account both the bluntness and the displace- 
ment effects, may also be explored. 

From the flat-plate problem, we see that the zero- 
order solution will be governed by only one variable 
and one parameter I’. Viscous hypersonic flows over a 
blunted inclined plate will accordingly be similar for the 
same I’. For thin bodies other than flat plates, inspec- 
tion of Eq. (2.15) shows readily that the parameters 
governing the similitude are T and (y,/A,?/*);. The 
correlations based on this similitude will be subject, of 
course, to the same errors as Eq. (2.15). Implicit in 
the above results are also the assumptions of the linear 
viscosity-temperature relation, unit Prandtl number, 
as well as uniform 77/7). Thus, for a more accurate 
correlation under more general conditions, the simili- 
tude based on our zero-order theory may not suffice. 
As a part of the study complementing the foregoing 
analysis, an examination of the hypersonic similitude 
under less restrictive conditions will be given in Section 


92 


a=... 


(2.3) Similarity of Hypersonic Flows Involving 
Leading-Edge Bluntness and Boundary Layers 


For both experimental investigations and practical 
applications, laws of similitude for correlating flow- 
field properties under different free-stream conditions 
and body dimensions are desirable. The similitude 
arrived at previously in Section 2.2 is simple but restric- 
tive. In the following, the more general viscosity- 
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temperature relation 4«7™ is assumed, the Prandtl 
number is taken constant but otherwise arbitrary, and 
the body-surface temperature is not necessarily uni- 
form. A perfect gas with constant specific heats is 
still employed but the specific-heat ratio is no longer 
required to be close to unity. 


Basic Requirements for the Validity of Hypersonic Bound- 
ary-Layer Similitude Involving Leading Edge Bluntness 

The (inviscid) hypersonic similitude which takes 
into account the effect of slight blunting has been given 
by one of us.** In arecent paper, Hayes and Probstein** 
considered the similitude involving boundary layers. 
From these two studies a viscous hypersonic similitude 
involving bluntness is not difficult to infer. Hayes 
and Probstein caution, however, that the similitude in- 
volving bluntness may not be extended to the viscous 
case, except under a condition so restrictive as to make 
the result uninteresting. In view of this, an examina- 
tion of the basic requirements which underlie hyper- 
sonic similitude seems proper at this point. 

Hayes and Probstein observ that the “‘equivalence 
principle,’ that is, the hypersonic small-disturbance 
theory, necessarily fails in the vicinity of the blunt 
nose itself, and it would appear that the layer of fluid 
of increased entropy nearest the body may not obey 
the law of similitude. It would also follow that the 
hypersonic boundary-layer similitude cannot exist. 
The major conclusion of our following discussion is that 
a full similitude capable of correlating the density and 
the temperature fields near the outer edge of the bound- 
ary layer is not possible, but for pressure and velocity 
fields, as well as for the density and the temperature 
within the boundary layer, the similarity laws do apply. 

Let us first examine whether similitude can exist in 
the inviscid flow region. On the basis of our discussion 
in the earlier part of Section 2.2, we see that under re- 
quirement (a), the approximation « ~ U holds through- 
out the inviscid flow region including the entropy 
layer. The applicability of the requirement (a), and 
the degree of approximation involved, have been ex- 
amined. From the fact that u ~ U we can conclude 


M (x.?/K.)?(Y./K.), 


Xe™ 


T/T) = f (x/L) 


with Yo Tis, 


This form may also be compared with the similitude 
reducible from Eq. (2.15) discussed previously. We 
see that the simplicity gained by the zero-order theory, 
which assumes a very strong shock and a very small e, 
consists essentially in the elimination of the parameters 
vy, Tw/To, and x,7/K,. 

There are two important limiting cases of the simili- 
tude given above. One is the case of negligible bound- 
ary-layer displacement effect, the other that of a very 
strong shock wave. For the first, the interaction param- 


(K./x.7)?(Pu/YVoPo) 


i K, 3 . ; 
P (0.661 + 1.73 7 ) x || (M*Cy/V 7 ) 
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that the pressure p and the velocities v and u do obey 
the “equivalence principle’’—hence, the hypersonic 
similitude.* For the density and temperature in the 
entropy layer, however, the ‘‘equivalence principle’”’ 
may not be strictly applicable since these quantities 
may depend on the entropy determined upstream at the 
nose shock. 

For the boundary layer, we observe that with the wall 
conditions specified the solution will be completely de- 
termined, in general, by the three quantities at its 
outer edge: the pressure p;, the total enthalpy //; 
(=H/..), and the velocity u;. Now, all these quantities 
obey the law of similitude of the inviscid region; in 
particular, we have u ~ U, and these are exactly the 
conditions required for the similitude of the hypersonic 
boundary layers in the sense of reference 44. We may 
note that the actual magnitude of the error in the heat- 
transfer rate, the skin friction, and the displacement 
thickness incurred by the approximation u; ~ U will 
not differ a great deal from |Au;/U'|. In fact, if we 
assume local similarity (which is appropriate for the 
present application) the corrections to these quantities 
will be, respectively, —1/2, —*/2, and +!/2 of | Au;/U 

Therefore we conclude that with the exception of the 
density and temperature fields in the entropy layer, the 
classical hypersonic similitude of Tsien and Hayes can 
be extended to include both bluntness and boundary 
layer, so long as the flow regions immediately down- 
stream of the leading edge are excluded and the flow 
Mach number is high, i.e., 


M> 1 (2.29) 


In the absence of the tip bluntness, the first of these re- 
quirements may, of course, be replaced by a less restric- 
tive condition related to the error due to slip, tempera- 
ture-jump, or vorticity effect. 

The inclusion of the tip-bluntness effect introduces 
an additional parameter \/*kt/L* into the boundary- 
layer similitude discussed in reference 44. The simi- 
larity law correlating the flow quantities of interest in 
the experiments to be described may be written as 


=f F : ( * ) Pe ike a , YY, @, Pr | (2.30) 
L’ \K2*), K. | 


eter M?+*/+/Re, will become merely a scaling fac- 


tor, and we have 
Mébx/L, M*®Cy, M°C, 
(Met+e/V Rez) f(x/L; Mr, M*kt/L, y, «, Pr) 
A discussion of the second case follows. 
Hypersonic Similitude Involving Strong Bow Shocks 
For this case, we require the shocks to be very strong 


* For the axisymmetric case, it will be meV k(t L) 
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so that terms of order 1//°6? are negligible in com- 
parison with unity. From the viewpoint of practical 
application, consideration of this special form is not 
trivial. Since the essential difference in basic require- 
ments between this similitude and the zero-order theory 


M(x.*/K,)*(Y;/K.x), 


9 


is the condition on y, the zero-order theory furnishes a 
convenient basis for comparison with data correlated 
by this similitude. 

Under the assumption of a strong shock the Rankine- 
Hugoniot shock conditions may be simplified. Using 
a suitable affine transformation, we arrive at a form of 
similitude comparable with Eqs. (2.29) and (2.30): 


(K./x2)°(Pu/ Yabo) | 


| x Xe : 
Tw\ K.} ao | :( = ) iT. 7,0, Pr] (2.31) 
€ | 0.008 -f 1.43 é ) | (M°C,y/V y ) ; L Ke . L ie 
( Xe 


with Y,,/7L, To/T = f(x/L) 


Comparing this with the more general similitude of 
Eq. (2.30), we see that in the above form of similitude 
the dependence on the parameter (x2/K,) or, what 
amounts to the same thing, kRe,/M'***, becomes very 
weak as a result of the shock being very strong. On 
the other hand, comparison with the similitude based 
on the zero-order theory shows that as a consequence 
of relaxing conditions on Pr, w, and e, we have now to 
include as invariants the parameters Pr, w, y, as well 
as 7,,/7). The dependence on y and on 77/7) in the 
strong-shock similitude cannot be too strong, however, 
since both are connected with the first-order e€ effect. 


(3) EXPERIMENTAL STUDIES OF 
FLATPLATE FLOWS 


The foregoing theoretical study of Section 2 has 
treated the general problem of combined tip-bluntness, 
viscous-interaction, and body-angle effects in hypersonic 
flow about slender bodies. In particular, specific 
solutions have been given in Section 2.2 for the general 
flat-plate problem. In addition, hypersonic viscous 
similitude has been extended in Section 2.5 to include 
leading-edge bluntness effects. The present section 
will describe the accompanying experimental program 
which was designed to provide experimental data for as- 
sessment of the extended similitude of Section 2.3, and 
for comparison with predictions of the theoretical solu- 
tions of Section 2.2. Detailed results of the experi- 
ments are discussed in Section 4. 


(3.1) Outline of Experiments 


The experimental program involved the measure- 
ment of surface heat-transfer distributions and shock- 
wave shapes in hypersonic airflows about sharp and 
blunted flat plates under cold-wall conditions such that 
the (laminar) boundary layers were highly cooled. 
These studies were made in the C.A.L. 11 X 15-in. 
hypersonic shock tunnel at free-stream Mach numbers 
in the range 10.5 to 12.5. The range of leading-edge 
and unit free-stream Reynolds numbers employed was 
large enough to encompass the extremes of dominant 
bluntness and dominant viscous-interaction effects. 
Instrumentation has also been developed for measure- 


ment of surface-pressure distributions, and surface- 
pressure studies are currently underway. 

Most of the experiments were conducted with air- 
stagnation temperatures in the range 2,000 + 100°K, 
at Mach numbers in the range 11.3 to 12.5. With a 
plate surface temperature of close to 300°K, the cor- 
responding wall-to-stagnation temperature ratios were 
within the range 0.15 + 0.01. For these conditions, 
tests were made at zero angle of attack for leading-edge 
Reynolds numbers from about 3 to 15,000. Stream 
Reynolds numbers per inch here ranged from about 
14,000 to 180,000, and leading-edge thicknesses from 
0.0002 to 0.203 in. 

Angle-of-attack effects were studied at low and high 
leading-edge Reynolds numbers over the incidence 
range —10 to +15°. 

Heat-transfer measurements were also made over a 
range of air-stagnation temperatures from 2,000°K 
to 4,000°K for a sharp plate at zero angle of attack. 
The ratio of wall-to-stagnation temperatures in these 
tests ranged from 0.15 down to 0.075. 


(3.2) Shock-Tunnel and Test Conditions 


All experiments were conducted in the C.A.L. 11 X 
15-in. hypersonic shock tunnel. A detailed description 
of this tunnel and the measurement techniques em- 
ployed in its operation is given in reference 59. In the 
initial experiments of the present study, conducted 
with helium as the driver gas, 14 ft. of 3 in.-i.d., low- 
pressure-driven tube was used with a 14 ft., 3.5 in.-1.d. 
high-pressure driver. In the later stages of the pro- 
gram, an additional 14 ft. of driven tube was added, 
thereby doubling the total driven-tube length. The 
additional tube length served to increase the useful test 
time obtained with hydrogen as the driver gas. It may 
be noted here that doubling the driven-tube length re- 
sulted in no detectable change in measured heat-transfer 
rates or pressures in the tunnel test section with helium 
driver-gas operation. This fact indicated negligible ef- 
fects on the test fiow from possible mixing at the 
tailored interface. 

The model tests were conducted in a two-stage nozzle 
described in reference 59, which was coupled to the end 
of the driven tube. The first stage of this nozzle is 
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Fic. 3.1. Second-stage 15°-wedge nozzle used for tests 


convergent-divergent and expands the flow in a horizon- 
tal plane to a 1.5 X 12.25-in. cross section. Prior to 
entering the second stage, the flow is turned 10” up- 
ward by means of a deflection wedge in order to centri- 
fuge out diaphragm particles. Final expansion in the 
vertical plane then takes place in a sharp-entrance, 
15°-wedge nozzle which terminates in an 11 X 15-in. 
test section. Fig. 3.1 shows a schematic drawing of 
the 15°-wedge nozzle and illustrates the general test 
arrangement and model locations. Further details of 
model location are discussed in section 3.3. 

The tunnel was operated under tailored-interface 
conditions with helium or hydrogen at room tempera- 
ture used as the driver gas. With tailored-interface 
operation, the driver-gas initial conditions are adjusted 
to avoid shock reflection at the driver-driven gas inter- 
face. Theoretical tailoring with room-temperature 
helium and hydrogen results for shock Mach numbers 
in air of about 4+ and 6.5, respectively. The corre- 
sponding (stagnation) temperatures behind the re- 
flected shocks are close to 2,000 and 4,000°K. The 
range of stagnation pressure, fy, behind the refected 
shock employed in the experiments was from 1,000 to 
10,000 psi. Experience has shown the quality of noz- 
zle reservoir conditions attained behind the reflected 
shock to be negligibly affected by small departures from 
theoretical tailoring. With the 2S-ft. driven tube, the 
duration of useful test flow was about 3 and 6 msec. 
for hydrogen and helium driver gases, respectively, 
following the initial period of 0.5 to 1 msec. required 
to attain steady, test-flow conditions. 

Shock-tube conditions measured in each experiment 
included initial air pressure in the driven tube, initial 
shock-wave speed, and the nozzle reservoir pressure py» 
attained behind the reflected shock wave. The nozzle 
reservoir conditions, including stagnation enthalpy, 
pressure, and density, are calculated from measured 
shock speeds with the use of real-gas, equilibrium-flow, 
shock relations. The measured pressure fp, obtained 
as a function of time by means of SLM piezoelectric 
and Norwood strain-gauge transducers, provides a 
positive and direct check on the calculated reservoir 
conditions. Over a range of values of py from 1,000 to 
10,000 psi, the average difference between measured 
and calculated values of pp is about 5 percent. 

Calibration studies of the test flow in the 15°-wedge 
nozzle have included simultaneous measurement of 
static pressures at the tunnel sidewalls and flow pitot 
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pressures across the tunnel by means of a multiple- 
probe rake, flow angularity studies with a multiple- 
wedge rake instrumented for heat transfer,*! and meas- 
urement of surface pressures on a sharp-edged plate at 
large compressive incidence such that interaction ef- 
fects are negligible. Static pressures at the sidewalls 
have been measured with Altec-Lansing microphones 
and also with barium-titanate piezoelectric transducers. 
The latter have also been employed in pitot-pressure 
rakes.* As regards test Mach number, the calibration 
studies for helium/air operation, 7) ~ 2,000°K, indi- 
cate an average value on the tunnel centerline at station 
A (Fig. 3.1) of 12.6 + 0.3 at a stagnation pressure po 
of 10,000 psi, 12.3 + 0.3 at 5,000 psi, and 11.3 + 0.2 
at 1,000 psi. Typical calibration results for Mach 
number variations with time at several positions along 
a diagonal in the plane of station A (Fig. 3.1) are shown 
in Fig. 3.2 for reservoir pressures of 1,030 and 4,350 psi. 
The open symbols indicate Mach numbers at the diag- 
onal positions shown as calculated from the ratios of 
the measured pitot (p)’) to reservoir (fo) pressures 
assuming isentropic, equilibrium flow. The local pitot 
pressures in this case were measured with a pitot rake 


* Most of the pressure instrumentation employed for calibra- 
tion studies was developed by the Hypersonic Wind Tunnel De- 
partment of the C.A.L. Experimental Facilities Division. Grate- 
ful acknowledgement is made to K. D. Bird and J. F. Martin, 
Head and Assistant Head, respectively, of that Department, for 


their cooperation in making this instrumentation available 
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Fic. 3.2. Typical calibration results for Mach number (17) sur- 
vey at Station A. Helium/air operation, 7) ~ 2,000° K. 








368 JOURNAL OF THE AEROS 


employing barium-titanate transducers. Also shown 
in Fig. 3.2 is the flow Mach number as calculated from 
the ratio of measured wall static pressure to reservoir 
pressure, the wall static pressure being obtained with a 
modified Altec-Lansing microphone located on the 
tunnel centerline at station A. 

Measurements of local pitot pressure and flat-plate 
heat transfer at different longitudinal stations on the 
15°-nozzle centerline have indicated very small longi- 
tudinal gradients over a distance of three to four inches 
Differences in heat transfer 
The maximum 


at the test location used. 
were not measurable over this distance. 
Mach-number gradient at the test location is estimated 
to be less than 0.1 per inch. It will be noted in Fig. 
3.1 that the test location is included within the theoreti- 
cal, uniform-flow test rhombus as calculated from the 
method of characteristics with allowance for boundary- 
layer development. 

Information on the flow angularity with respect to 
the 15°-nozzle axis has been obtained from surveys with 
a multiple-wedge rake instrumented for heat-transfer 
measurement.®' More recently, further information 
on the flow angularity has been obtained from the re- 
sults of extensive heat-transfer studies with a yawed 
slender cone.” The results of these tests indicate the 
centerline flow at station A to have an angularity of 
—().8 degrees in the pitch plane (downward in Fig. 3.1) 
and 0.6 degrees in the yaw plane. 

For individual experiments, the test-flow properties 
at station A were calculated from the test Mach number 
M, the reservoir pressure fp, and the stagnation en- 
thalpy //.,, assuming isentropic expansion with ther- 
modynamic equilibrium maintained. At the stagna- 
tion temperature of 2,000°K and the relatively high 
stagnation pressures employed (exceeding 1,000 psi), 


RUN 588 He /AIR 
M, = 3.9 
MODEL B: oc = 15° COMPRESSION y= 0° 


PACE SCIENCES 


MAY 1961 


oxygen dissociation is quite negligible and molecular 
vibration is the only inert degree of freedom excited. 
The assumption of a test flow in thermodynamic 
equilibrium for these conditions is considered valid on 
the basis of calculations of typical vibrational relaxa- 
tion effects to be expected in the first expansion stage of 
the 11 X 15-in. tunnel nozzle. These calculations 
were carried out for pure nitrogen at a stagnation tem- 
perature of 2,000°K over a range of reservoir pressures 
appropriate to the present experiments. The nitrogen 
molecules were assumed to be harmonic oscillators ini- 
tially in complete thermodynamic equilibrium at the 
reservoir conditions. In subsequent expansion through 
the nozzle, a Boltzmann population of vibrational 
states was assumed to persist, characterized by a lag- 
ging vibrational temperature. This model for relaxa- 
tion from an initial Boltzmann distribution is consistent 
with the findings of Montroll and Shuler®* for Landau- 
Teller collisional transition probabilities. For the par- 
ticular geometry of the first expansion stage, the results 
of the calculations show vibration-relaxation effects to 
be relatively small and probably unimportant at 
reservoir pressures of 1,000 psi, and quite negligible at 
higher pressures. 

At the stagnation temperature of 4,000°K and the 
lower stagnation pressures employed, an appreciable 
amount of oxygen is dissociated. Departure from 
equilibrium of oxygen dissociation in the nozzle expan- 
sion must be considered at these conditions, as the 
enthalpy involved is of the same order as that in nitro- 
gen vibration. As part of a separate program of non- 
equilibrium-flow research underway at C.A.L., calcula- 
tion studies of nonequilibrium nozzle airflows involving 
oxygen dissociation have been made on the basis of 
various simplified kinetic models of air.*-®* For a stag- 
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Typical drum-camera oscilloscope record. 
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Flat- plate model “A” 


nation temperature of 4,000°K and the present pres- 
sure levels, these studies indicate the effects of dissocia- 
tion nonequilibrium on test-section conditions to be 
relatively small. While some freezing of atomic oxygen 
may be expected to occur, the amount of enthalpy 
frozen in is not sufficiently large to substantially affect 
the thermodynamic and gasdynamic state of the flow. 
This is particularly true at the higher stagnation- 
pressure levels of the present tests, i.e., around 5,000 
psi. 

As regards the possibility of air liquefaction in the 
nozzle flow, the lowest static stream temperature at 
station A was around 70°K, attained with helium as 
the driver gas and a stagnation temperature of around 
2,000°K. This static temperature was well above 
the equilibrium liquefaction limit at the pressures in- 


volved. 
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A typical drum-camera record from the multichannel 
Hathaway oscilloscope employed for recording is shown 
in Fig. 3.3. The record shows the time histories of 
reservoir pressure py») and sidewall static pressure p 
at station A, in addition to surface temperature change 
with time at various locations on a flat-plate model at 


station A. 


(3.3) Heat-Transfer Models and 
Instrumentation 


Two flat-plate models were employed in the experi 
ments. Figs. 3.4 and 3.5 show the pertinent details 
of the two models which will henceforth be designated 
as A and B. Both models are made of The 
leading-edge wedge angles are 30° for A and 20° for B 
These relatively large angles permitted location of a 
heat-transfer gauge within */; in. of the leading edge 
(see below) and also facilitated honing of the leading 
for the sharp-edge 


steel. 


edges to a thickness of 0.0002 in, 
tests. 

For heat-transfer 
fitted with thin-film resistance thermometers of plat- 
inum deposited by air spray on a thin Pyrex insert, 
as indicated in Figs. 3.4 and 3.5. The insert of Model 
A contains a central row of 15 elements. That of Model 
B contains elements at several spanwise stations as well 
as the central row in order to provide a check on span- 
The Pyrex inserts and steel were mated 


measurements, each model was 


wise effects. 
very carefully, particularly at the leading edge, to 
minimize the step. 

The heat capacity of such thin-film thermometers is 
sufficiently small that the films record instantaneous 
surface temperature change with time within very 


small errors.*? Local heat-transfer rates were obtained 
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Fic. 3.5(b). 


by numerical (IBM 704) solution of the one-dimen- 
sional heat-flow equation with the measured surface 
Further de- 
tailed information on the thin-film circuitry and record- 
ing methods employed with the 11 X 15-in. C.A.L. 
shock tunnel is given in reference 60. 


temperature-time histories as input data. 


Model calibration included determination of the 
product pC,« for the Pyrex plate inserts, as well as de- 
termination of the resistance-temperature curve for 
each film. The product pC,« is required in the calcula- 
tion of heat-transfer rate from 
temperature change. 


measured, surface- 
Calibration of pC,x with the 
electrical pulse method® gave values of V pCox in close 


agreement with the value 0.0745 Btu/ft?—°F pre- 





Flat-plate model ‘‘B.”’ 


viously determined for Pyrex by Vidal.*? This value 
was used in reduction of the data. 

The resistance-temperature curve for each film was 
determined by a film resistance measurement at two 
temperatures with the assumption that the slope was 
constant. The assumption of constant slope has been 
shown to be good over a very wide range of film temper- 
ature change. In the present tests, the surface 
temperature rise in the time period of interest did not 
usually exceed 10°F. Film resistance was checked at 
frequent intervals. Only the slope of the resistance 
temperature curve enters the determination of heat- 
transfer rate. 
were found to be exceedingly stable. 


Over a period of about a year the films 
With Model A, 
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for example, the maximum change during this time in 
the resistance-temperature curve slope of any one ele- 
ment was about 3 percent, and for most elements it 
was less than 1 percent. 

Regarding film durability, the major loss of films 
occurred during the first few runs with the newly in- 
strumented models. This experience is common with 
thin-film instrumentation. Model A _ proved very 
durable, with only one of 15 elements being lost. On 
Model B, however, 12 of the 20 central elements were 
eventually lost as well as several of the outboard ele- 
ments. 

Model A was employed for the initial, sharp leading- 
edge tests and for all bluntness tests. For the bluntness 
tests, the leading-edge material was progressively re- 
moved to a maximum leading-edge thickness of about 
0.010 in. Further blunting up to 0.203 in. was done 
by cementing on thin slabs with sharp 90° and 120° 
corners maintained by subsequent light grinding and 
hand finishing [see Fig. 4.2(d)]. The leading edge of 
Model B was not blunted. Leading-edge thicknesses 
f were measured on an optical comparator. For the 
smallest thicknesses, 0.0002 to 0.0096 in., comparator 
measurements were made of impressions of the leading 
edge taken with soft solder across the model span. The 
thickness of 0.0002 in. represents the best estimate ob- 
tainable for the sharpest leading-edge condition and is 
to be considered as an upper limit on ¢. 

Both models were supported by stings from a cir- 
cular-are sector as indicated in Figs. 3.1, 3.4, and 3.5. 
At zero angle of attack, the instrumented surface of 
both models lay in the tunnel centerplane. The leading 
edge of A was located at station A. The leading edge 
of. B was about one inch upstream of station A. At 
angle of attack, the leading edge of A remained fixed 
in the tunnel centerplane while that of B moved a 
maximum distance of ().5 in. off the centerplane at a = 
15°. 

Angles of attack (a), roll (zero), and yaw (zero), 
were mainta:ned to about +0.10 degrees to the geo- 
metric axes of the tunnel. 
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Fic. 4.1. Typical variation and reproducibility of measured heat 
transfer. 
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Schlieren Apparatus 

The schlieren system employed was the conventional, 
off-axis single-pass type with 12-ft focal length, 16-in. 
diameter, parabolic mirrors. A short-duration barium- 
titanate spark was the light source. A magnification 
somewhat greater than one was used to increase detail. 
The collimated schlieren beam was very carefully 
aligned with the model in roll and yaw with the aid of 
sighting wires located at extreme spanwise stations 
Calibration photos which included a length calibration 
grid were taken prior to each run. 

Shock-wave coordinates relative to the plate were 
measured from the photographic negatives on an optical 


comparator. 


(a) 
t=.0002" 
M=12.3 
Re, =18.8 

Re, =7.2x10"/1" 


(b) 
t». 0096" 
12.3 
Re, = 7h 


Re, #7.4x10"/ i" 


(c) 
t=. 049" 
Me12.3 
Re .=3770 
Re,-7.7x10"/ 1m 





i 
Ake 
4 i 


& 





(d) 
t».203" 
M=12.3 
Re.-15000 
Re, =7.4x108/18 





Fic. 4.2. Schlieren photographs of shock-wave shapes (a = 0°, 
Ty) ~ 2,000°K). 


(4) EXPERIMENTAL RESULTS 


(4.1) General 


Some general aspects of the experimental results will 
be discussed prior to considering the detailed findings. 

Fig. 3.3 shows typical oscilloscope recordings of 
plate-surface temperatures. With constant surface 
heat-transfer rates, the temperature at any point in- 
creases as the square root of time. With the small ratio 
of wall-to-stagnation temperatures of the present tests, 
and the small rise of surface temperatures in the short 
test-flow period, the changes in surface temperature 
have negligible effect on the local heat-transfer rates. 
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Also shown in Fig. 3.3 is the output from a simple 
analog computer network which solves the one-dimen- 
sional heat-flow equation.” The input signal for the 
analog network is the output voltage from a thin-film 
element. The output signal of the computer, shown 
in Fig. 3.3, is the local heat-transfer rate g. Such 
analog networks have recently been incorporated into 
all recording channels. 
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(4.2) Schlieren Results 


Schlieren studies were carried out at zero angle of at- 
tack for the nominal conditions described in Section 
3.1. At the test Mach number of around 12 and 7) ~ 
2,000°K, the range of Re, covered was such that 8 


-_"" 
interaction effects would dominate down to small values 


varied from large values for which viscous- 
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Fic. 4.3. Summary of shock-wave shapes (a = 0. 7, ~ 2,000°K) 


In general, at least two tunnel runs were made for 
each set of experimental conditions. Fig. 4.1 shows a 
plot of local heat-transfer rate g versus distance down- 
stream from the plate leading edge for two successive 
runs with closely equal stagnation temperatures and 
pressures. The reproducibility here is perhaps some- 
what better than that usually obtained. The average 
scatter obtained in g was about +5 percent. 

The distribution of g in Fig. 4.1 is typical for all con- 
ditions studied. The magnitude and distribution of ¢ 
indicate a laminar boundary layer. Changes in gq in- 
dicative of transition were not observed. Early experi- 
ments with Model A having a length of 4 in. showed 
somewhat high values of g at the most downstream ele- 
ments, followed by a very rapid drop in g. It was 
suspected that these peculiarities were caused by some 
type of trailing-edge influence, and they were substan- 
tially eliminated on lengthening the model to 7 in. 

Tests with Models A and B for the sharpest leading 
edges (0.0002 in.) under the same tunnel conditions 
gave the same results within the normal scatter. This 
indicated small effect due to the difference in the two 
wedge angles of 30 and 20 degrees, respectively. As re- 
gards spanwise variations of g, tests with Model B over 
a wide range of reservoir pressures showed no signifi- 
cant spanwise variations across the central 1.5 in. of 
the test surface. 


for which inviscid bluntness effects would dominate. A 
valuable aspect of the schlieren studies was the direct 
assessment they provided of the degree to which the 
experimental flows met requirements (a) and (d) of the 
theory (Section 2.2) as expressed in terms of the local 
shock-wave angle. 

Typical schlieren photographs for a range of Re, at 
nearly constant Ke, per inch are shown in Fig. 4.2. 
The outer region of the boundary layer is evident in (a) 
and (b). With increased blunting, the density de- 
creases adjacent to the surface. For sufficient blunt- 
ing, the boundary layer is not detectable at the low- 
density levels obtained. The background striations in 
the photographs of Fig. 4.2 are caused by the relatively 
poor schlieren quality of the test-section windows em- 
ployed. These windows produced negligible image dis- 
tortion. 

A comparative summary of typical shock-wave 
shapes deduced from photographic negatives is given 
in Fig. 4.3. Here the behavior of the shock shape is 
seen to be qualitatively similar to that observed by pre- 
vious investigators (e.g., reference 49) for hypersonic 
flow with adiabatic-wall conditions. For sufficiently 
large Re;, the flow is nearly inviscid and the shock shape 
becomes dependent only on K,. At sufficiently small 
values of Ke, the shape is determined essentially by 
mutual interaction of the shock wave and plate bound- 
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ary layer. The effect of the boundary layer for small 
Re, may be seen by comparing shapes H and I in Fig. 
4.3 which were obtained for the same leading-edge 
thickness ¢ but for different free-stream Reynolds num- 
bers per inch. With decreased Reynolds number the 
boundary-layer displacement thickness and hence, the 
effective body thickness is increased, which pushes out 
the shock wave. A significant difference between the 
present studies and previous adiabatic-wall investiga- 
tions is the fact that in the present experiments the 
boundary layer was highly cooled. The large cooling 
thins the boundary layer in comparison with the adia- 
batic-wall case, and thereby reduces the displacement 
effect. 

For sufficiently large values of Re, the hypersonic 
similitude for inviscid, perfect-gas flows with bluntness** 
should apply. The similitude indicates that the shock 
shape should then correlate as 


Y,/t = kM2f(x/M°kt, y) 


even for the weak-shock case. Correlations of this type 
have previously been observed for helium test gas re- 
Fig. +.4 shows correlation of the present air re- 
t versus x/t at the three largest 


Also shown for comparison 


sults.* 
sults in the form Y 
values of Re, employed. 
is the limiting theoretical prediction from the present 
zero-order theory of Section 2.2 for small 8, and a pre- 
viously obtained strong blast-wave result.* The cor- 
relation of the experimental shock shapes is good. The 
agreement with the zero-order prediction, which as- 
sumes a strong shock wave and y close to one, is within 
At small x/t, the large fiow de- 
At large 


about 10 percent at best. 
flections cause departure from the theory. 
v/t, departure results because of progressive deteriora- 
tion of the strong shock condition required by the 
theory, i.e., 1/?(d Y,/dx)? > 1. 

A general correlation of all experimental shock shapes 
in terms of the variables x{A,~*MY,/x and B = 
given by the strong shock similitude for a per- 
In reducing 


ia 
fect gas (Section 2.2) is shown in Fig. 4.5. 
the schlieren data to this form, stream conditions were 
obtained as discussed in Section 3.2. Values assumed 


4 . 0096 


GLAST wave 
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Correlation of shock-wave shapes for large Re: (a = 


0°, Tyo ~ 2,000° K). 
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Fic. 4.6. General correlation of shock-wave shapes for arbitrary 


B(a = 0°, JT ~ 2,000°K 


for y and the leading-edge drag coefficient k were 1.40 
and 2, respectively. The correlation is fairly insensi- 
tive to the value of k. The value of C used in evaluat- 
ing x, was based on the reference temperature 7% 
suggested in Section 2.2. Correlation of the data in 
terms of variables applicable for a perfect gas flow is 
reasonable in the present situation. With the large 
boundary-layer cooling present, the maximum tempera- 
ture in the boundary layer was only of the order of 500 
to 700°K (Jo ~ 2,000°K). The inviscid flow field 
external to the boundary layer was several hundred 
degrees cooler than this except in the immediate stagna- 
tion-point region with large bluntness Assuming 
equilibrium flow about the plate (no vibrational lag), 
real-gas effects would be quite small where flow de- 
flections are sufficiently small to satisfy the requirement 
of the theory—i.e., away from the stagnation zone. 

In addition to requiring the usual conditions for 
hypersonic small-perturbation theory, 1J > 1 and 
6< 1, a correlation in Fig. 4.5 also requires the condi- 
tion of a strong shock wave. The lack of correlation of 
certain data points can, in general, be attributed to 
failure of the experiments to meet these conditions. In 
order to emphasize this aspect of the experimental 
data, the data points have been divided into two groups 
on the basis of whether or not they meet the criterion 


2/M < dY,/dx < 0.30 
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Re: (a = 0). 


These arbitrary limits on d Y,/dx result from the condi- 


tions 
(1/2)(dY,/dx)?—/" < 1/4 (a) 
M?*(dY,/dx)? > 4 (b) 


considered to be reasonable minimum requirements cor- 
responding to (a) and (0) in Section 2.2. The filled 
symbols in Fig. 4.5 are those for which dY,/dx falls 
within the above bounds. It can be seen that these 
points correlate reasonably well. The remaining points 
violate condition (b) far downstream and condition (a) 
close to the leading edge. In particular, it will be 
noted that all the data points for the experiment with 
t = 0.0002 in. and Re, = 14.3 fail to meet condition 
(b) above. 

Also shown in Fig. 4.5 for comparison with the ex- 
perimental data is the zero-order result of the present 
theory. In addition to requiring the above conditions 
to be satisfied, the zero-order result requires small e, 
or y close to unity. The agreement with the experi- 
mental data is seen to be fairly good over a wide range 
of 8 where the assumptions of the theory are most valid. 
In particular, the largest discrepancy between experi- 
ment and theory occurs for ¢ = 0.0002 in. and Re, = 
14.3, where no data points satisfy the strong-shock limit 
(b) above. 


(4.3) Heat-Transfer Results 


For purposes of correlation and comparison with the 
present zero-order theory for a perfect gas, a local heat- 
transfer coefficient Cy = g/p.U(H., — H,) was used. 
Values of other parameters employed in reduction of 
the data were y = 1.4, k = 2,and C = C(Tx). 
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A general characteristic of the correlations to be dis- 
cussed (and also of Fig. 4.5) is the fact that they are 
not too sensitive to the drag coefficient k, leading-edge 
thickness /, or to the stream Mach number J/ and the 
constant C for fixed stagnation conditions py and JZ... 
Changes in J at fixed py and H.., of course, also produce 
corresponding changes in p..U and Re,. This relative 
insensitivity also applies to a considerable degree to 
the measure of agreement between theory and experi- 
ment. The insensitivity results because changes in the 
parameters concerned generally affect both ordinate 
and abscissa in such a way as to displace points in a 
direction more parallel to the correlation curve than 
normal to it. As a consequence of this behavior, pos- 
sible small errors in the test Mach number at given /) 
and /7., would have very little effect on the present 


correlations. 


Small Re;, a = 0 

The present similitude study (Section 2.3) indicates 
that for the strong shock-wave condition, \/*Cy for 
a = () should correlate in terms of the usual viscous- 
interaction variable x = MV C/V Re, at large values 
of 8B = x,K.~?/*, for given values of y, H,,/H., etc. 
Correlation for a weaker shock will show an additional 
dependence’on K,x.~*. If K.x,~? is sufficiently small, 
however, even the weaker-shock case should correlate 
in terms of x. In terms of the present experiments 
where the strong-shock condition is met only to a mar- 
ginal degree, this means that a correlation of J/*Cy in 
terms of x could be expected only for the lowest values of 
Re. 

Typical experimental results in the form J/*Cy 
versus x are plotted in Fig. 4.6 for a range of stagnation 
temperatures up to 4,000°K. The values of Re, range 
from 2 to 200, approximately. Also shown are theoreti- 
cal predictions for the strong-interaction region given 
by Lees,'® Whalen,” and the present zero-order theory, 
as discussed in Section 2.2. The zero-pressure-gradient 
prediction of A. D. Young (reference 63, page 422) has 
been evaluated for the typical experimental conditions 
M = 12.3, H,/H. = 0.14, Pr = 0.72, and w = 0.80. 
The latter prediction is quite insensitive to this particu- 
lar choice of conditions. 

In general, the data of Fig. 4.6 for 2,000°K stagna- 
tion temperature correlate reasonably well. With the 
data scatter obtained, it was difficult to detect a posi- 
tive influence of increasing Re, for values of Re, less than 
about 300. From the strong-shock similitude, the rela- 
tive importance of viscous-interaction and bluntness 
effects depends only on the variable 8, provided the 
shock is strong. For the highly cooled boundary layer, 
bluntness effects at a given station x are relatively larger 
than in the adiabatic-wall case, since large cooling thins 
the boundary layer and reduces the severity of viscous 
interaction. This was first pointed out by Lees and 
Kubota,’ and is evident from examination of the effect 
of H,,/H.. on 8 as discussed in Section 2.2. If a value 
for Re, of 300 is taken as that indicating a measurable 
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bluntness effect in the present experiments, the corre- 
sponding value of 8 at x = 1 in. downstream from the 
plate leading edge is about 0.4. For this value of 8, the 
present zero-order theory predicts an increase in W*Cy 
of about 10 percent due to leading-edge bluntness. 

The data of Fig. 4.6 for stagnation temperatures be- 
tween 2,000 and 4,000°K provide an extended check 
of the similitude as well as a check of possible real-gas 
effects. Equilibrium real-gas phenomena reduced the 
Mach number of these tests to the range 10.5 to 10.9. 
The high-temperature data correlate well with those 
for low temperature, indicating that the wall-to-stagna- 
tion temperature ratio and real-gas effects have small 
net influence on the heat transfer. This behavior is 
not unexpected, since in the similitude of the zero- 
order theory the dependence on //,///., disappears 
completely. Further, real-gas effects in the flow about 
the plate are minimized for the sharp leading-edge 
condition. Such real-gas effects as occur may be ex- 
pected to be correlated in terms of the heat-transfer 
coefficient used which is based on a total enthalpy 
rather than a total temperature difference. 

As regards the agreement of theory and experiment 
in Fig. 4.6 for small Re,;, the present zero-order theory 
lies below the experimental data. The difference is 
attributed principally to failure of the experiments to 
meet adequately the strong shock condition required 
by the theory, and also to the fact that the experimental 
values of y and Pr are not too close to 1. Effects due 
to velocity-slip, temperature-jump, and shock-wave 
heating, which are not accounted for in the theory, are 
considered to be small. The shock is strongest close 
to the plate leading edge where the slope is largest. In 
this region, the data points in the largest range of x 
(from 15 to 27) show the best agreement with theory, 


being some 15 to 20 percent high. This difference 
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becomes larger with increasing distance downstream as x 
and the shock slope decrease. The agreement with 
Whalen’s more exact calculation” based on the self- 
similar solution and y = 1.4, 7/7) = 0.2, and Pr l 
is excellent at the largest values of x 

Comparison of the data in Fig. 4.6 with the zero- 
pressure-gradient result of A. D. Young suggest that 
the data may lie in a transition region whereby J/*°Cy 
tends from the variation proportional to x at low x 
(Young) to the variation proportional to x** at high 
x (Lees, etc.). It was desired to obtain low Mach 
number (low x) data to determine whether such data 
would fair in smoothly to the zero-pressure-gradient 
theory at sufficiently low x. However, tests could not 
be carried out at sufficiently low Mach numbers for this 
purpose in the 11 X 15-in. tunnel. As an alternative, 
it was decided to reduce the data obtained for small ¢ 
at large positive angle of attack (compressive flow), 
and high density, in terms of local 1/*Cy and x values 
based on conditions behind the leading-edge shock wave 
as computed for inviscid flow. At sufficiently high 
density behind the shock, the viscous interaction effects 
are small and the zero-pressure-gradient result should 
apply. 

The filled symbols of Fig. 4.7 are the positive angle-of- 
attack data (7) ~ 2,000°K) plotted in this form, i.e., 
M,*Cy, versus x,, where state (1) is that behind the 
leading-edge shock for inviscid fiow. The data for 
a = Oand 7) ~ 2,000°K of Fig. 4.6 are also included 
as the open, unflagged symbols. The local values of 
X = x, go down to about 0.4 for a = +15 At this 
end, the experimental results do appear to fair in quite 
well with the zero-pressure-gradient theory. Likewise, 
the data for a > () blend quite smoothly with those for 
a = (), although at lower a, where interaction effects 
are present, a rigorous correlation on this basis is not 
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Fic. 4.8. Correlation of experimental heat transfer for large Re: 
(a = 0°, Ty ~ 2,000°K). 


to be expected. The good agreement of the experimen- 
tal data with the classical, zero-pressure-gradient result 
for the smallest values of x, affords a good check on the 
experimental technique. 

Also shown in Fig. 4.7, by the flagged symbols, are 
data for negative angles of attack (expansion flow) 
reduced in the same manner—i.e., the reference state 
(1) is that behind the leading-edge Prandtl-Meyer ex- 
pansion fan for inviscid flow over the plate. There 
appears to be no reason to expect x, to correlate the 
data in this case, and the degree of correlation obtained 


is somewhat surprising. 


Large Re,, a = 0 

At the other limiting extreme of very large Re, for 
a = 0, the general similitude (Section 2.3) indicates 
that MW*C;,/x should depend only on A, (for given 
values of y, Pr, H,/H..) regardless of whether the shock 
is strong or not. In this case, x serves only as a scaling 
factor. For the strong shock wave and y close to one, 
the zero-order result in this limit gives 1/*Cy propor- 
tional to K,!/*x (Section 2.2). 

Fig. 4.8 shows typical experimental data for 2,000°K 
stagnation temperature plotted in the form J/*Cy 
versus K,'/*x. Results are given for the highest value 
of Re, obtained, namely 15,000, as well as for the value 
3,700 obtained at two different values of ¢ differing by a 
factor of about 4. Also shown are the present zero- 
order result, and a self-similar solution employing the 
plane blast-wave pressure distribution’ with y = 1.4 
and Pr = 1. 

The correlation obtained in Fig. 4.8 is considered 
fairly good. The trend of the results for the higher 
density level, i.e., for the higher value of Re, per inch, 
is to shift to the right with increasing KRe;, as expected. 
At a lower density level but the same value of Re, = 
3,700, however, the value of M/*°Cy is reduced. A con- 


tributing factor here is that the strong shock condition 
is met best in the case of large leading-edge thickness, 
Regarding comparison with the theory, the experi- 
mental values of 1/*Cy, for the highest values of Re, 
are some 20) to 25 percent higher than the zero-order 
result, except in the leading edge vicinity where experi- 
ment and theory cross. Here the flow deflections be- 
come unduly large, and the heat transfer is presumably 
reduced by the associated shock-wave heating effect. 
Away from the leading-edge region the slope of the data 
points is quite close to the theoretical slope. At Ke, 
= 15,000, the shock slope is everywhere larger than the 
arbitrary limit of 0.30 previously discussed in connec- 
tion with the schlieren results of Fig. 4.5. However, 
the resulting shock-wave heating from this effect would 
tend to reduce the heat transfer rather than increase it. 
The strong-shock-wave requirement of the theory is 
met reasonably well for ¢ = 0.203 in. The error in- 
curred in the zero-order theory from requiring y close 
to one is such as to predict a heat-transfer level which 
is too high, as discussed in Section 2.2 and illustrated 
by the similar solution for the plane blast-wave pres- 
sure distribution, y = 1.4, and Pr = 1, shown in Fig. 
4.8. For the values of A, appropriate to the Re; = 
15,000 case, the plane blast-wave pressure distribution 
has been shown to be in very good agreement with that 
calculated exactly for inviscid flow from the method of 
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Fic. 4.9. General correlation of experimental heat transfer for 
arbitrary B (a = 0°, 7) ~ 2,000°K). 
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Fic. 4.10. Correlation of experimental heat transfer for small Re, at angle of attack (7) ~ 2,000°K 


characteristics.°*? The discrepancy between theory 
and the correlated experimental data may be attributed 


in part to the assumption of Pr = 1 in the theories. 


General Correlation for Arbitrary 8 = x-K,.~?*, a = 0 


Provided that the strong shock-wave condition is 
satisfied, the simplified hypersonic similitude shows 


that for a = 0, the quantity 


«(0.664 + 1.73 /He) yo! 2(K,/ x2)? M Cx 


should be correlated in terms of the variable 8 = 
x.A,-*/8 (for given values of y, Pr, H,./H») when both 
bluntness and viscous-interaction effects are impor- 
tant. A general presentation of the experimental heat- 
transfer data for zero angle of attack and 7) ~ 2,000°K 
is given in this form in Fig. 4.9. The zero-order result 
applicable for a strong shock wave and y close to one is 


also shown. 


In general, the overall correlation obtained in Fig. 
4.9 is considered to be reasonably good. The data 
appear to correlate somewhat better at the lower values 
of 8 (large bluntness effect) than at the higher values 
(large viscous interaction effect). As previously dis- 
cussed, at the higher values of 8 the strong shock-wave 
condition is not well satisfied in the experiments. With 
weak shock-wave conditions, an additional dependence 
of the ordinate of Fig. 4.9 on K,x,~? is to be expected 
This additional dependence on A,x,~* is observed in the 
data of Fig. 4.9 and generally accounts for the data 
separation at intermediate and large values of 8. 

As regards agreement with the zero-order theory, the 
experimental heat transfer is generally higher than pre- 
dicted. At intermediate and large values of 8 this dif- 
ference may be accounted for by failure of the experi- 
ments to meet the strong-shock condition, by the fact 
that the theory is subject to an error of order e, and by 
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the fact that Pr ~ 0.71 for the experiments. For the 
data points of any one experiment at intermediate and 
large values of 8, it is seen that the best agreement with 
the theory is obtained closest to the leading edge where 
the strong-shock condition is best satisfied. The agree- 
ment here is typically in the range 10 to 15 percent. 
At low values of 8, the experimental heat transfer is 
typically some 20 percent to 25 percent higher than 
the theory, as discussed in the previous section for large 
Re,. For the largest values of ¢ in this range of 6, the 
experimental heat-transfer points cross the theoretical 
curve in the vicinity of the leading edge. As previously 
discussed, this behavior probably results because of the 
large flow-deflection angles at a fixed station x in the 
leading-edge region when ¢ becomes large. 


Small Re, at Angle of Attack 
For sufficiently small values of Xe, and a strong shock 
wave, the present viscous similitude predicts that the 


quantity 
[0.664 + 1.738(4H,,/H..) |\Cu 
€ | | 
a| IV y 


should be correlated in terms of J/°a*x,~! for small 
angles of attack a. Fig. 4.10 shows a plot of the ex- 
perimental data for small Re; 7) ~ 2,000°K, in this 
form along with the theoretical zero-order result for y 
close to one. The angle-of-attack range includes data 
fora = +5, +10, +15, —5, and —10°. In addition, 
data obtained at nominal zero angle of attack have been 
plotted for an arbitrary, small value of a of 0.01 radians. 

The correlation obtained in Fig. 4.10 is fairly good. 
For a positive a of 5° there is a suggestion of bluntness 
effect at the higher value of Re,. This is understand- 
able in terms of the parameter [ = K,.Ma/2x? ~ 
Rea/ M? (Section 2.2) which controls the viscous simili- 
tude for a blunted plate at incidence. The correlation 
of the negative angle-of-attack data is better than 
anticipated in view of the fact that the strong-shock 
condition is rather poorly met for expansion flow over 
the plate. 

The experimental heat transfer in Fig. 4.10 is gener- 
ally higher than predicted by the zero-order theory. 
The agreement with the theory is generally best for 
positions close to the plate leading edge where the 
strong-shock condition is best satisfied. It will be noted 
that at a = 15° the experimental data agree quite well, 
as expected, with the zero-pressure-gradient result of A 
D. Young calculated for Pr = 0.72. Young’s value is 
about 20 per cent higher than the limiting, zero-order 
result of Section 2.2 for large M?a?x,7"'. 


Large Re, at Angle of Attack 


The data plot analogous to Fig. 4.10 for large Re, 
and dominant bluntness effect at angle of attack is 
shown in Fig. 4.11 along with the zero-order prediction. 
The similitude indicates that the quantity (1/°/x,) X 
V ck iy lal ® C, should be correlated in terms of 
M'a?K,~*/8 in this case. The boundary-layer dis- 


placement parameter x, = MVC V Re, appears asa 
scaling factor in the ordinate. 

A correlation in Fig. 4.11 requires small flow defiec- 
tions as well as the strong-shock condition and large 
Re, The correlation obtained and the agreement with 
the zero-order theory are fairly good, particularly for 


a < (), in view of the experimental limitations. As in 


the case of small Re,, the experimental heat transfer is 
greater than predicted. The marked falloff in heat 
transfer, relative to the theory, close to the plate leading 
edge for negative angle of attack is undoubtedly due to 
the large flow deflections in this region and local failure 
of the small perturbation assumption. This effect is 
evident at a = 0, as illustrated in Fig. 4.9 for the 
smallest values of 8. The falloff is increased for a < 0) 
and decreased for a > 0. This behavior results from 
change in the effective bluntness of the square leading 
edge due to change in the stagnation-region flow pattern 
with angle of attack. The effective bluntness de- 
creases for a > () and increases for a < 0 because of 
shift of the flow stagnation line on the body. 


(5) CONCLUDING REMARKS 


A theoretical and experimental study has been carried 
out on combined leading-edge bluntness and boundary- 
layer displacement effects in hypersonic flow over thin 
bodies. 

A theoretical approach has been developed based on a 
continuum flow model consisting of three adjoining 
regions: an inner laminar boundary layer, an outer 
(detached) shock layer, and between these two, an en- 
tropy layer. Solutions for this flow model have been 
obtained by considering the mutual interaction of 
these three regions under the requirements (@) to (c) 
stated in Section 2. In the boundary layer, a small 
perturbation procedure in e = (y — 1)/(y + 1) is de- 
veloped which yields the local flat-plate similarity as a 
leading approximation. In the entropy layer, the low 
density leads to the simple pressure-area relation Eq. 
(2.11). In the detached shock layer, the results of 
existing shock-layer theory are applicable. The zero- 
order theory in this approach corresponds to the classi- 
cal result of Busemann. 

The zero-order theory corresponding to the leading 
approximation under requirements (4) and (c) has been 
developed in detail; its simplicity permits solutions to 
numerous specific problems of interest. In general, the 
solutions reveal the combined effects of bluntness, vis- 
cous interaction, and angle of attack in hypersonic 
flows over thin bodies. Continuous transitions are 
exhibited between limiting solutions which agree with 
existing blast-wave, strong-viscous-interaction, and 
hypersonic-wedge theories. One interesting feature 
found for inviscid flow over a blunted wedge is 
an oscillatory pressure decay over the afterbody sur- 
face. This result may be interpreted as one limiting 
form of shock-layer reimpingement at glancing inci- 
dence. 

Hypersonic similitude has been examined in Section 
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2.3 for the present problem where both bluntness and 
boundary-layer displacement effects are present. The 
classical similitude of Tsien and Hayes can be extended 
to include these effects for correlation of pressure and 
velocity fields, as well as flow quantities of interest 
within the boundary layer, provided the vicinity of the 
blunt leading edge or nose is excluded. For flows in- 
volving strong shock waves, a simplified similitude is 
obtained. The 
strong shock and also small e, may be readily compared 
with experimental data correlated in this similitude. 
The theoretical results of the present paper pertain 
Certain aspects of 


zero-order results, which require a 


chiefly to the case of plane flows. 
the treatment can readily be applied to flows with rota- 
tional symmetry, for example, flows about a blunted 
cone, and to three-dimensional boundary-layer flows 
over wing systems. Results from these applications of 
the theory are discussed in reference 65. 

On the experimental side, studies have been made of 
hypersonic airflows at Mach numbers around 12 over 
sharp and blunt leading-edge flat plates for conditions 
of large boundary-layer cooling. These studies were 
made in the C.A.L. 11 X 15-in. hypersonic shock tun- 


nel mainly at stagnation temperatures close to 2,000°K, 
giving wall-to-stagnation temperature ratios of the 
order of 0.15. Shock-wave shapes were obtained by 
schlieren photography, and surface heat-transfer dis- 
tributions were measured by thin-film resistance ther 
mometry. At zero angle of attack, the leading-edge 
thickness and stream Reynolds number were varied 
over a sufficient range to encompass the limiting cases 
of dominant bluntness and dominant viscous-interac- 
tion effects. Leading-edge Reynolds numbers varied 
from about 3 to 1.5 X 10‘ for stream Reynolds numbers 
varying from about 1.4 XK 10% to 1.8 X 10° per inch 
Angle-of-attack effects from +15 to — 10° were studied 
Heat- 
transfer measurements were also made on a sharp plate 
at zero angle of attack for air stagnation temperatures 
up to 4,000°K. 

The experimental data have been reduced and plotted 
in forms indicated by the general hypersonic similitude 
of Section 2.3 for viscous flow with leading-edge blunt- 


at low and high leading-edge Reynolds numbers 


ness effects. In addition, comparisons have been made 
within this framework with the predictions of the zero- 
order theory given in Section 2.2. 
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The results of the schlieren studies of shock-wave 
shapes permit an assessment of the degree to which the 
experimental flows met the assumptions of the theory as 
expressed in terms of the local shock-wave angle. In 
particular, the similitude used for arbitrary bluntness 
and viscous interaction, and the zero-order theory, both 
require the condition of a strong shock wave in addition 
to small flow-deflection angles. For given stream 
Mach number, these conditions place upper and lower 
bounds on the local shock angle. 

Generally speaking, when due allowance is made for 
local failure of the strong-shock or small-flow-deflection- 
angle conditions in the experiments for a = 0), the 
schlieren and heat-transfer data for 7) ~ 2,000°K are 
found to be correlated reasonably well by the similitude 
forms of Section 2.3 over a very wide range of varying 
bluntness and viscous-interaction effects. For the 
present experimental conditions, at least, the correla- 
tion obtained suggests that the hypersonic viscous 
similitude is valid even with leading-edge bluntness. 
The general heat-transfer correlation in the bluntness- 
dominated region appears somewhat better than that 
in the displacement-dominated region, probably because 
the former meets the strong-shock condition best. The 
correlation of heat transfer for angle of attack is found 
to be quite good. Good correlation is also obtained for 
heat transfer on a sharp plate at stagnation tempera- 
tures from 2,000 to 4,000°K. 

As regards agreement of the theoretical predictions 
of Section 2.2 with experiment, the zero-order theory 
for a = 0) predicts heat-transfer rates generally lower 
than those measured. However, the discrepancies 
appear to be largely accounted for by the limitations of 
the experiments in meeting the condition of a strong 
shock wave, by the fact that the theory is subject to an 
error of order ¢, and by the fact that Pr ~ 0.71 for the 
experiments. Where the strong-shock condition is 
best satisfied, close to the leading edge, the zero-order 
result for heat transfer is typically some 10 to 15 per- 
cent low. 

For the angle-of-attack results at small and large Re,, 
the theoretical, zero-order heat transfer is typically 
some 15 or 20 percent below that measured. At large 
compressive inclination giving high-density flow over 
the plate, the measured heat transfer checks well, as 
expected, with the classical prediction for zero-pressure 
gradient. The latter value, calculated for y = 1.4, 
Pr = 0.72, and other conditions pertinent to the ex- 
periments, lies some 20 percent above the correspond- 
ing zero-order result. 

Finally, it should be stressed that the present experi- 
mental studies have been made at essentially one Mach 
number only. In view of the parametric importance of 
Mach number, further experiments at substantially 
higher Mach numbers in air are desirable. It is antici- 
pated that further development of the hypersonic 
shock tunnel will provide the high enthalpies necessary 
to carry such studies into the Mach number range of 
15 to 25 in air. In this range it can be expected that 
the bluntness-displacement problem will eventually be 


influenced to a significant degree by complex real-gas 
phenomena. 
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Summary r = m*Dj/L,*De, bending shear stiffness ratio of the 
panel 
The flexural motion of a rectangular sandwich panel with an s = vy2/vr2, Viscous orthotropic index of core 
orthotropic viscoelastic core is investigated. Differential equa- t = time tl 
tions and boundary conditions are derived by taking into account u,oru, = displacement in the x; direction SU 
the rotational inertia of the whole panel, and the shearing de- v, w) q 
formation and viscoelastic damping of the core. The dynamic U, V, W = displacement response to a unit step function ti 
response of a simply supported rectangular panel to a unit step x(or x, = rectangular Cartesian coordinates (see Fig. 1) 
fanction is first obtained, and then the general transient and y, 5) as 
steady-state response of the panel to an arbitrary input function dias = (N,,*m? + N,,*n62){(1 + f)(m? +. 262)? ti 
is expressed in terms of the unit response by means of the con- (1 + K)(1 + »)}~, ratio of the edge compres 
volution integral. It is shown that three families of modal sion to the elastic buckling load of the panel es 
shapes of vibration are possible for each nodal pattern of the corresponding to the wave pattern m and n st 
panel. The frequencies and amplitudes of vibration of these B = L,/L,, length-width ratio of panel : 
three modes are discussed in detail for the case of an isotropic 5 = variational sign 06 
and elastic sandwich panel. The effects of viscoelastic core and ” = r(1 + f)(m? + n%8*), shear deformation param th 
edge compression are also discussed. Two simple physical eter be 
parameters which represent the effects of shearing deformation u = Poisson’s ratio of facing en 
and rotational inertia, respectively, are found very essential in the p’, p", p = mass density of upper facing, lower facing and ra 
analysis. A number of limiting cases are discussed in detail core : 
a . ° . a 
Vr2, Vyz = coefficients of viscosity of the core P 
py = (1 + f)(p’h’ + K2p"h") + (1/3)1 — K 4+ KF? th 
Symbols ph = (1 + K/h)? times the moment of inertia sti 
: : an : per unit area of the sandwich panel with respect ta! 
A, B, C= amplitudes of vibration : 7 : I ‘ 
D EWI +4 2/1 + KV 2) | i tiff to an axis parallel to the xg axis and lying in da 
) = E’h’h* 2/(1 - (1 — p?), bending stiff- . , : ‘ 
setter - as the neutral surface of the panel 
ness per unit width of sandwich panel : : ‘ tet 
fA me ‘ , - ‘ Pr = mass per unit area of sandwich element 
De = hG,A(1 + f)*, shear stiffness per unit width of * aie ; : ene ; is! 
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Fic. 1. Sandwich panel under loading (W/xx, Myy, Wxy and sup- 
ports are all omitted for clarity) 


that this classical theory cannot give satisfactory re- 
sults for the cases of sharp transients or for the fre- 
quencies of vibration of high-order modes. Of par- 
ticular importance is Mindlin’s! flexural theory of plates 
as deduced from the three-dimensional theory of elas- 
ticity. 

A sandwich panel is a composite structure consisting 
essentially of a relatively thick, low-density, and low- 
stiffness core bound between two thin sheets of high- 
stiffness material. Because of the softness of the core, 
the shearing deformation in a sandwich panel under 
bending can be appreciable and because of the heavier 
concentration of mass at the outer facings which are 
rather far apart, the rotational inertia of a sandwich 
panel may well be of considerable importance. It is 
the purpose of this investigation to study the steady- 
state and transient vibrations of a sandwich panel 
taking these factors into account. In addition, the 
damping of a sandwich panel due to the bonding ma- 
terial between the core and the facings is simulated by 
assuming that the sandwich core exhibits viscous re- 
sistance to shearing deformation. 

The problem of the static bending and buckling of a 
sandwich panel has been treated by the authors in a 
previous paper.” The sandwich construction to be 
considered here is similar to that of reference 2, i.e., 
the sandwich panel is rectangular in shape; the two 
facings are of different thicknesses and the core is 
orthotropic. The loadings acting on the panel are 
uniformly distributed time-independent edge forces 
and moments and arbitrarily distributed time-depend- 
ent transverse loads. Fig. | shows a sandwich panel 


under the specified loadings. 


Assumptions 


(a) The forces in the plane of the panel are taken 
entirely by the two thin facings of the sandwich panel. 
Each individual facing acts as a membrane and does 
not take bending. 

(b) The core of the sandwich panel carries the main 
part of transverse shear loading. The distribution of 
this transverse shear is approximately constant through- 
out the thickness of the core. The core behaves like 
an orthotropic viscoelastic continuum, being char- 
acterized by its two shear moduli and two coefficients 


of viscosity. 


w 


(c) The deflections of the panel are small, and are 
constant across the panel thickness. 

(d) The upper facing and the lower facing have the 
same Poisson’s ratio. 

Assumptions (a) and (b) regarding the direct stress 
taken by the facing and the transverse shear taken by 
the core are widely adopted assumptions based on the 
relative stiffness of the core and the facings. The 
assumption that the facings are membranes is equiva 
lent to neglecting quantities of the order of the square 
of the ratio of facing thickness to core thickness in 
comparison with unity. The assumption that the 
transverse shear is constant throughout the thickness 
of the core is a direct result of assumption (a) that the 
core does not take forces parallel to the plane of the 
panel. For homogeneous plates, the transverse shear 
varies considerably across the thickness of the plate. 
Some authors! * introduced a numerical factor in the 
relationship between the average transverse shear 
stress and strain to take into consideration this non- 
uniform distribution of transverse shear. Yu!‘ in- 
vestigated the appropriate value of this numerical 
factor for sandwich plates from a more general theory 
and showed that this factor is very close to unity 
for most common sandwich panels. In assumption (b), 
the sandwich core is assumed to exhibit viscous resist 
ance to transverse shear. This is to simulate the 
damping of the bonding material. The assumption 
of constant deflection across the panel thickness simpli- 
fies our analysis. This assumption, however, excludes 
a possible symmetrical mode of vibration of the panel 
about its neutral surface. For most structural mate- 
rials the Poisson’s ratios differ by a small amount. By 
making the last assumption, that the two facings have 
identical Poisson’s ratio, we can show’ that the lateral 
displacements in the two facings are related by the 
simple relationship 


—Ku', (1) 


u” —(E'h'/E*h")u’. 


Qa 


This simplifies our analysis a great deal. 


Analysis 


The present problem is different from the static prob- 
lem of reference 2 in the following respects: 

(a) The inertia of the panel, both translational and 
rotational, is considered. 

(b) The sandwich core exhibits a viscous shear stress 
which is proportional to the time rate of transverse 
shear strain 

(c) The bending rigidity of each individual facing is 
neglected. 

As in reference 2, the behavior of a sandwich panel 
can be completely characterized by the displacements 
u’, of the middle plane of the lower (or upper) facing of 
the sandwich panel and the deflection u; of the panel 
The deflection “3 is assumed to be constant across the 
thickness of the panel. The displacement components 
U, are assumed to vary with the panel thickness as 
shown in Fig. 2. It should be emphasized that the u, 
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are displacement components associated with the bend- Eh’ iat , pl 
: g , ; ; Quu’,,. — p)( 
ing of the panel. No coordinate is provided for the 2(1 — 2?) (1 + f)2uw'y ya + (1 B)( aap eq 
displacements due to the direct edge compression and 1+K tic 
shear. In contrast to the static problem of reference t's as) | — Ga ( F u', — U -) - fu 
7 : ] 1 
2, the displacement components u; are functions not ex 
only of the space coordinates but also of time. ‘+s, \ = ’ 9 eq 
a4 . ‘ : i ; Va3 ; U at — U3, } = SPrll a,tt (2) 
lhe differential equations of motion of the panel in h / 1+ K th 
terms of these displacements and their space and time . . pl. 
pees : ates ; (No summation on a) 
derivatives can be found by considering the equations pu 
of motion of a differential element of the panel.* A Likewise, the equation of motion in the transverse tré 
differential element is shown in Fig. 2 together with direction is wi 
the forces acting on it. The law of motion applied - i+z ., K’ 
to the differential element for the two rotational degrees q Vas U3,a3— NGas j U aa U3 aa} — eq 
of freedom in the x1-x3; and the x»-xv; planes gives the sa 
following two equations: . l+A , . £. 
hPa U aat— U3, aat " Pi U3 4 (9) T 
h r 
aa P _ ; : th 
Similar to reference 2, the remaining equation of rota- | 
_ ; = tion is satisfied identically by the symmetry of stress 
These equations, together with the boundary conditions, . Va : » ‘ ‘ ! 
were also obtained from the generalized Hamilton’s principle. tensor in the facings and the other two translational 
Because of lack of space the energy integrals are not shown here. equations reduce to Eq. (1) between the lateral dis- 
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Fic. 2. (a) Variation of displacements wa across panel thickness (assumed);  (b) Inertia and internal forces acting on a sandwich ele- + 
ment in flexural motion (forces due to stretching of panel not shown). 
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LOADED SANDWICH 


placement components in the two facings. In these 
equations repeated Greek subscripts indicate summa- 
tion over the range 1 and 2. The notations used are 
fully described in the section on symbols and are not 
explained in the text. In deriving the two rotational 
equations of motion, Eqs. (2), the internal forces in 
the facings are assumed to be attached to their middle 
planes. The rotational inertia of the facings 1s com- 
puted as if the mass of the whole facing were concen- 
trated at its middle plane. This is reflected in our 
writing p, = (1 + f)(p’h’ + K*p"h") + (1/3) —-K+ 
K*)ph. Notice that p,{h/(1 + K)] is approximately 
equal to the moment of inertia per unit area of the 
sandwich panel with respect to an axis parallel to the 
X, axis and lying in the neutral surface of the panel. 
The position of the neutral surface is determined by 
the factor K and is shown in Fig. 2. 
The boundary conditions are 


}[E’h'/2(1 — w2) JAC + Sf) [2un’,,6a3 + 
(1 sn w)(U" a + U's «) | <5 M ap} Ng = 0 (4) 


or 6u’ Ng = 0 


as l+K 
Nas U3 8 hGs (u, * oe h u's) — 
1 


; . > : 
Nes \ U3,at — -  Uat { ta = 0 (9) 
h 


or 6u3; = 0 


These boundary conditions correspond to the three 
boundary conditions necessary for the complete deter- 
mination of the flexural problem of a plate when shear- 
ing deformation is considered.'* As written in the 
subscript notation, these boundary conditions are 
readily applicable to panels of arbitrary shapes. 

The dynamic behavior of the flexural motion of a 
sandwich panel is therefore completely characterized 
by the three displacement components u;(%, X2, ¢) 
which are solutions of the differential equations (2) and 
(3) subjected to the boundary conditions (4) and (5). 


In terms of these new units and dimensionless parameters the three 


r d7Umn™ _ Gta» * 
‘-_ “_ + 7*(] + K) 


ie d,*? d,* 


] 
r1i+f/adt+ w( 2 mn3) San” 


l 
+ (1+ K)} ri + (m: + - 
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Simply Supported Rectangular Panel 


For a simply supported rectangular sandwich panel, 
the displacement components u; can be assumed in the 
following forms: 


, mmx , ny 
» = ru > Umn(t) COS sin 
m=l1 n=1 re in 
mmx nry 
, ° é 
, = p p Umn(t) sin cos 
m=1 n=1 } ae iL. 
_ mmx | ny 
~~ = Zz z Wmn(t) sin - sin = 
m=1 n=1 L. L, 


where m/2L, and n/2L, are the wave numbers defining 
the nodal pattern of the panel and u,,,(¢), Umn»(t), and 
Wmn(t) are the amplitudes of vibration associated with 
the particular nodal pattern. These displacements 
satisfy the alternate set of boundary conditions sug- 
gested by Hoff for simply supported rectangular sand- 
wich panels. 

If we expand the forcing function g in (3) in terms 
of the eigenfunctions of the region as g(x, y, t) = 


i . Manx , 
> > dmn(t) sin I sin 


m=l1n=1 “z 4y 
sion together with the displacements given above into 
the equations of motion (2) and (3), we will obtain three 


nTry , ; 
and substitute this expres- 


coupled second-order ordinary differential equations 
for the amplitudes of motion Upp(t), Umn(t), ANd Wy»(t). 
It is found that these equations can be greatly simpli- 
fied and made dimensionless by the introduction of 
certain dimensionless parameters and a set of new char- 
acteristic units. The new set of basic units is defined 


as follows: 


Time: L, = the time for plane elastic 
VE! p’(1 — pn’) wave to travel the 
length of the lower 
facing 
Length: zh = m times the thickness of 
the core 
Mass: ’h'(xh)? = the mass of a square of 
lower facing with side 
length equal to 7 times 
the thickness of the core 


equations for the amplitudes of motion become 


: Hn's*) +)1 | témn* + 


m\ dWmn* m 
T : - Wnn = 0 (6a)T 
L,*) de (* ) . 


e @0mn™ : | ag a ~ 97° l — M 9 * 
— p,* —— + r*s(1 + K) + 1 + K)Y (Ql +f) 2°8? + —— om?) + 8 | Omn™ + 


1 d,*? d,* 


r(1+f)(1 + K) 


+ uy nB \ dWmn* ng 
5 Omens Umn* — T*S : —g - | Wmn* = 0 (6b) 


L,*j] d,* Ll.” 


+ Starred quantities are dimensionless quantities in terms of the new set of basic units. 
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¢ A *Wnn* — 7*(m* + 50*6*) dting® | k + gn*B° 


ee Pt é,"* L,**(1 + K) d,* . ie l+K 


m \ dumn* m nB \ dimn* np r 
7 * * coms * ion” — r*s a g oo —= _ " ely (6c) 
/ 5 d, ie tae ae : ag T 


Note that in these equations every symbol is dimensionless. 
In order to solve these equations for the response #»*(t), Umn*(t), and Wm,*(t) to an arbitrary gm,*(t) we shall first 


r(N,.*m? + Nyt) | Wmn* — 


find the response of the system to a unit step function 


(0), t< oD 
Hnn*(t) = J * (7 
1, t>0 


by the method of Laplace transformation. If we designate the response to the unit step function by the upper case 
alphabet Unn*(t), Vinn*(t), and Wn,»*(t), replace gmn*(t) by Hmn*(t) in (7), and take the Laplace transformations of 
Eqs. (6) from the physical ¢-space to the transformed p-space, the following set of transformed equations is obtained 
under the condition that the initial displacements and velocities of the panel are zero: 


| —_ 
) - pp? + r*(1 + K)p + (1 + A) ri + f) (wm: + e u's?) + i} + 
“ 


U...* 
mn 2 


Vian” yr + fi + K) ( : . mn — Wan* Fie (r*p + Wt =() (Sa)? 
] y 
Una” yl + f)(l + &) ( . mus) + Vn* } 4 pp? + r*s(1 + K)p + 


Vr 


; vee iu ] Oe ) 
(1+ K){]r(1 + f) {2°82 + : m*}) + g iia Wmn™ 7 = (v"sp + £¢ () (Sb 


we J n * a t oat —_. J np ¥ ¢ g ( W. * J r * 9 r*(m* r snp") 
 ¢ 7.,* (r*p + 1) Ven \L,* LF" sp => ar To wn a Pi*p L,*(1 + K) 
l m” + gn°B? ns Slici 4 
al ae Net vita) | : iz _- oP 


This set of three linear simultaneous algebraic equations in three unknowns can be solved by the well-known 
Cramer’s rule. The solutions for Un,*(p), Vinn*(p), and Wmn*(p) in (8) can be expressed as fractions of polynomials 
in p. The singularities of them are in general simple poles and the inversion of these Laplace transforms to get 
Unmn*(t), Vnn*(t), and W,,,,*(t) can be performed easily when geometrical and elastic details of the sandwich panel 
are given. 

When Unn*(t), Vinn*(t), and W,,,*(t) are known, the response of the panel in a particular nodal pattern to an arbi- 
trary input ¢»,(¢) which starts at zero at ¢ = (), can be expressed by the convolution integral as 


t 1 ag | Py a are 
Sina (t) = a Oma" (t — T)d7 (9a) 
J, dl 
ets). of Pape 
Una (t) = = Vinn*(t — T)d7 (9b) 
1, @i 
bd j */( 7 
ee aes 
Wmn*(t) = } —— Wan*(t — T)d7 (Se) 
Jo dl 


By summing up the contributions associated with each nodal pattern we obtain the total response of a sandwich 
panel to an arbitrary forcing function 
, _ mmx , my 
q*(x, y, ) = . Jmn*(t) sin sin 
1 


m=l1n “7 +y 


ft A bar under a quantity denotes the Laplace transform of that quantity;e.g., Umn = Umn(p) is the Laplace transform of Umn = Umn(t 








6c) 


irst 


ase 


- of 
1ed 


LOADED SANDWICH PANEL RESPONSES 387 





at ~ marx , nry (* ddnn*(1 
as ua" = y ® 5 cos SI ; Uma" (t — T)dI 10a) 
m=1n=1 ai dl 
; — << |. MAX nry (*' dg,,,*(1 ; 
y” > z sin cos - ; Vinn' (t — 1 )d7 Lb) 
ln=1 , — dl 
> * 
A Ty... mx. ney Can (i — 
u Za La SD sin fig ( 1 )d1 lQe 
m=! dl 


e 


Sandwich Panels with Isotropic Core--A Simple Example 

When the core of a sandwich panel is isotropic, the system of Eqs. (S) is simplified, and certain general properties 
of the solution can be obtained and interpreted in terms of the parameters of the panel without having to know the 
explicit numerical values of the parameters. 

Using Cramer’s rule, we obtain the solutions to the system of Eqs. (S) as follows: 


| " 1/3 Ai(p) 
| n'(p) = (lla) 
(pr/pr)(1 + K)? A(p) 
, ) Ao(p) 
| ma"(p) = cans See (11b) 
(p:/pr)(1 + K)? A(p) 
, r ? A:;( >) 
Winn" (Pp) = - f (llc) 


(p:/pr)(1 + K)? A(p) 


where A(p) = (Pp) (2 ) + (?) Twy(1 + o) + (2 ) (: Tits cnn + 
Wy Wy Wy + U] 
(C) mow (0-5 a) 
TWy( ND) | — Am) + 10(1 — am, | (11d) 
wy l+n 
l 
Ai(p) = yior(;" 3) (2 ) (lle) 
L 
+- 
Ax(p) = vio) (™ 5) (2 * ‘) (11f) 
= 40) (0) (ee) LCS) + (C) m+ + ag 
(Dp) (Pp) p 1+K a a TUN) ( n) § 
v(p) Tey oC)ee (er 51) (11h) 


. G.- 
= (1 + K) Vj. 
1p, 


The frequency «) and the nondimensional parameters 7, ¢, and a,,, 


(111) 


have very important physical interpretations 
They are to be discussed in later sections. 

The common factor ¥(p), which could have been cancelled froni the A's, is kept here because it represents a possible 
pure shear mode of vibration to be discussed later. 


The singularities of U’,,,.*, Vm,* and IW’,,,* are in general simple poles at different values of p, 


p () la 
Pp pi Yi T 10 P Pi Yi — 10) 2b 
P pe Y2 T le Pp py Y2 — lw I2¢ 


where p;, /1, P2, and /» are roots of the quartic appearing in A(p). 


The residues at the pole p = () are 


mi." - 
; (l3a 
w’(1 + f)(m? + n’B")?(1 — am 
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“aGl.* 

R: zoe a = ae - (13b) 

ie m?(1 + f)(m? + n?B?)?(1 — amn) 
W..,* i + K)Q + 7) (13c) 
mn: ae a Ge vs 209\9 -_) 

w(1 + f)(m? + n7B?)7(1 — amn) 

; 1 f)(m? n*B?)2) —! 
where mn = (N,-*m? + N,,*n?B?) lial baal l 


l (1 + K)(1 + 7) 


f 








is the ratio of the edge compression to the elastic buckling load of the panel corresponding to the wave pattern m,n. 


The residues at the other poles are 


: r/m m Pa { 
r..*: >a ss x 1} 44 
(pt/pr)(1 + K)? = [rn (2) ad | ( 
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= cnn (2) + (raw) (nd) (; = “=| (?=)t (13d) 


I+y 


‘ r/ x np Pe j ( 
Sue’: err i ~e 0 z ] 7 4 
(p:/p,)(1 + K)? (*) E (2+) r | \ 


) (?) 
+ 3rao(1 + o) ~ 


Wo 


aa Pa , : p _ Amn Pa t 2. 
L+n ome) = + ry (: i + :) (2)} a 
r/ p\? *) 
y = - Sa 3 " a l 
W mn - (pr p)(I rs K)2 |(2) + TH (2 + ( + n) 


) (2) + (zw) (nd) (: — ) (2+) | (13f 
Amn TU) —_ r ot) 
n wo otis l+n wy J J , 


\4 (22) + tli +o (2+) + 


The inverse transform of U;,,*(p), Vinn*(p), and Wmnn*(p) can be expressed in terms of these residues as® 


Umn*(t) = Residue at p = 0 + > (Residue at p = PalePat (14a) 
Vinn*(t) = Residue at p = 0 + bi (Residue at p = p,)ep,t (14b) 
Wrmn*(t) = Residue at p = 0 + > (Residue at PD = PadlePat (14c) 


(Pa = Pr, 1, pe pz) 


The first term on the right-hand side of each equation 
is independent of time and represents the deformation 
due to a statically applied load of unit intensity. The 
other members of the equation represent the transient 
response of the panel owing to the sudden application 
of the load. The real part y represents the damping, 
and the imaginary part w the damped natural fre- 
quency of vibration of the panel. 

It is seen that for each nodal pattern of the panel 
there are two families of modes of vibration represented 
by the frequencies w; and w:. For these modes, the 
amplitude ratio 


Une (#) Vua t = Mm np (a@ = l, 2) ( 15a) 
is the same for both modes, while 


Umn™ (t) 


J Vm n (*) (t ) 


(m/L,*) [two(po/ao) + 1] (15b) 
rere "> roe ae" a meaaneRS rer aear 5b. 

(1 + K)[(po/w0)? + rwo(pe/wo) + (1 + n)] 
depends on the value of ~,. When damping exists 
(r ~ 0), p, is a complex quantity, and the ratio U,,,(t)/ 
Wan(t) as given by (15b) is generally also a complex 
quantity. This indicates that there is generally phase 








difference between the displacement components u 
(or v) and w because of the different amount of damp- 
ing of these displacements. 

There is also another possible mode of vibration which 
is not excited in an isotropic panel by the loading in a 
direction perpendicular to the panel surface. It is 
described by 3, where 


bs he i 
_* 5 | =r + Yeon)? — 4 ] +> e é (15c) 


are the roots of the equation ¥(p) = 0. This mode of 
vibration is characterized by the amplitudes 
Wmn(t) = 0, Umn(t)/Umn(t) = —nB/m_ (15d) 


and represents a pure shear vibration without trans- 
verse deformation. 


The Three Modes of Vibration of an Isotropic 
and Elastic Sandwich Panel Under No Edge 
Forces 


In the preceding section, we have shown that three 
different modes of vibration are possible in a simply 
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supported isotropic and viscoelastic sandwich panel 
under quite general loading conditions. In this sec- 
tion, we shall give a detailed description of these three 
modes of vibration for the case in which the sandwich 
panel is not subjected to any edge forces and the core 
of the sandwich panel is purely elastic, i.e., there is no 


damping. 


(a) Frequencies of the Three Modes 

The frequencies for the three modes of vibration can 
be obtained from Eqs. (11d), (12b), and (12c). For 
the present case these frequencies, designated as «, we, 
and w;, for each particular value of m and n, are given 
in the following: 


baja ~ 1/9 +<4+6 = 
V(1 ++ ¢)? — 4n¢] (16a) 
(w2/wo)? = (1/2)[(1 + 9 +) + 
V(i+n+¢)? — 4n¢] (16b) 
n/(w/wo)? = (1 +n +) - 


(wi/wo)? (16c) 
W3 2 l — : 
(“) = i +4 5 n (16d) 
where w = (1 + K)V G.. ho, 


can be interpreted as the frequency of pure shear vibra- 
tion of an infinite panel. 
The three frequencies satisfy the inequalities 


We F W3> (17) 


It is of interest to see that when shearing deformation 
and rotational inertia are neglected, the only finite 
frequency as given by Egs. (16) is 


(w1),=¢=0 = (4?/L,”)(m? + n28°)V Da, (1 + fp. (18a) 


We shall call this frequency w,. It differs from the 
classical flexural frequency for homogeneous plates 
only by the appearance of the small parameter f. It 
can easily be shown that the following relationship 


exists: 
w,/wW = WV nb (18b) 
W/O = We/ Wo (18c) 


Notice that 7 and ¢ are functions of m and u which 
depend on the nodal pattern. Therefore the fre- 
quencies w;, we, w3, and w, also depend on the nodal 
pattern. The formulas for these frequencies in this 
section are for each particular wave pattern m and n. 
For all the possible patterns there is a total of triply 
infinite frequencies. 


(b) Amplitudes of the Three Modes 

The amplitudes of the motion in a particular nodal 
pattern excited by a unit step function perpendicular 
to the panel surface are obtained from Eqs. (13) and 
(14) as follows: 


Uh mL,/ mh 
rh w?(1 + f)(m? + n°B*)? 
> Amn” COS tig (19a 
Vinn(t) nBL,/ mh 
i y . - ae a 
ah w?(1 + f)(m? + n*6?)? 


bo 


> 2D Pan “OOS Wal (a (19b) 
a 


Wnt) (z,.” wh?) (1+ K)(1 + n) 


— = =< + 


wh w2(1 + f)(m? + n?B?)? 


> Con OS ae (19¢c) 


a 


where 
= 1M pr/ pt 
Ans =-— —— ; . (19d) 
w(1 + K)? (L,/mh) [(w,/an)* — no] 
fal rn p,/ pr 
Rag = me ; (19e) 
w(1 + K)? (L,/mh) [(w./wo)* — no] 
. te rpr/pill + 9 — (we/an)? 
Cn = (19f) 


r(1 + K) [(wa/a)! — no] 


In these equations A pm», By, Cn“? are the ampli- 
tudes of Mode 1, and An,, Brn, Cron are the 
amplitudes of Mode 2. There exist the following 
simple relationships among these amplitudes: 


Aga" Baa" (we/wy)' — no 
_ <a < U0 (20a) 
Ana i e (@1/a@)* — nod 
Ams" Ana ms 
” ~~ U (20b) 
Bun |: ae n3 ' 
Cen Rios ; (wy wo)” — (1 + n) : 
acre : rn >U (20c) 
Cus™ Ams *) (w2/wo)” — (1 + n) 
fon mhm/LAi + K) 
’ _ - > 0 (20d) 
on (1 + ~~ (@y We) ~ 
pe rhm/LA1 + K) 
= <@ (20e) 


Coa” (1 + 9) — (w2/wo)? 


Mode 3, which is not excited by loads perpendicular 
to the panel surface, is characterized by the amplitudes 


Wmn(t) = 0 { 


ae ; > 20 
Umn(t)/Umn(t) = —nB/m < Of (20f) 


With the above information about the amplitudes 
and frequencies, these three modes of vibration can 
be readily described. Fig. 3 illustrates the distinctive 
features characterizing these three modes. Among the 
three modes, Mode 1 has the lowest frequency of vibra- 
tion. The displacement components u, v, and w are 
in phase, as can be seen from Eqs. (20b) and (20d). It 
can also be shown from Eq. (20d) that A mnp/CnrQ? < 
1, indicating that Mode 1 is predominated by flexural 
deformation, even though some extent of shearing de- 
formation exists. It is this Mode 1 that corresponds to 
the classical flexural vibration. 
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Mode 2 has the highest frequency of vibration. Eqs. 
(20b) and (20e) show that the displacement components 
u and v are in phase, while w is out of phase with w and v. 
It can be shown from Eq. (20e) that | Amn / Cn | can 
be quite large for small ¢, which is generally the case. 
In other words, the motion of Mode 2 is predominantly 
shear in the plane of the plate. This type of motion is 
called thickness shear by Mindlin.’?’ Because of this 
predominance of thickness shear the concave face of 
the panel is stretched while the convex face is shrunken. 
This is contrary to Mode 1 and to the classical bending 
of plates. It can be shown that the shearing deforma- 
tion of Mode 1 and Mode 2 are in phase, even though 
the u (or v) displacements of these two modes are out 
of phase. This is because the shearing deformation 
of Mode 1 is mainly due to the deflection w, while the 
shearing deformation of Mode 2 is mainly due to the 
displacements uw and v. 

Mode 3 has the intermediate frequency of vibration. 
The deformation associated with this mode is pure shear 
in the plane of the plate, as indicated by Eq. (20f) that 
Wmn(t) = 0. Eg. (20f) also shows that the two non- 
ranishing displacement components u and v are out of 
phase. This shows the distinction between the shear- 
ing deformation of Mode 2 and Mode 3—that is, if the 
shearing deformations of these two modes are in phase 
in the x-z plane, they are out of phase in the y-z plane, 
and vice versa. In other words, in Mode 2, if one of 
the faces of the sandwich panel is stretched, it is 
stretched in all directions, while in Mode 3, if that face 
is stretched in the x-direction, it must be shortened in 
the y-direction. Fora simply supported isotropic panel 
the third mode is not excited by loads perpendicular to 
the panel surface. For other edge conditions, how- 
ever, all three modes are coupled together and the 
third mode will be excited by loads perpendicular to the 
surface of the panel. 

It should be noted that the above description of the 
three possible modes of vibration holds exactly only 
for isotropic panels. For orthotropic sandwich panels, 
Mode 3 is also associated with some extent of flexural 
deformation and the forces in a direction perpendicular 
to the panel will excite Mode 3, even though its ampli- 
tude is small. These results will be discussed in a 


later investigation. 


The Nondimensional Parameters 7, ¢, and K 


We have shown that the frequencies and ampli- 
tudes of the flexural vibrations of a sandwich panel can 
be expressed in a very simple way in terms of two 
parameters of the panel, 


n = (Do/Dg)(1 + f)(m?/L,")(m? + np?) 


1 ( ) _ (m? + 7B? 
é = 2 = AOE == O-) 
unc p wli+k 


“7 


The factor common to 7 and @¢, (m°/L,")(m* + n?6°), 
is inversely proportional to the square of the wave- 


length; Dy/Dg is the ratio of the bending rigidity to 
the shearing rigidity of the panel; (p,/p.)[i/(1 + A)? 
is approximately the ratio of the mass moment of 
inertia to the mass per unit area of the panel. There- 
fore the parameter y is directly proportional to the 
ratio of the bending rigidity to the shearing rigidity of 
the panel and inversely proportional to the square of 
the wavelength. The parameter @ is directly pro 
portional to the ratio of the mass moment of inertia 
to the mass per unit area of the panel-and inversely 
proportional to the square of the wavelength. The 
parameter 7 measures the relative importance of the 
bending rigidity to the shearing rigidity of the panel 
and is designated as the shearing deformation param- 
eter. The parameter ¢ represents the relative impor- 
tance of the rotational inertia to the translational 
inertia of the panel and is called the rotational inertia 
parameter. For the case of the classical analysis of 
plates, shear deformation and rotational inertia are 
neglected, and therefore n = @ = 0. The greater the 
values of » and ¢, the greater are the effects of shear 
deformation and rotational inertia. Notice that for 
a particular panel the ratio n/¢ is fixed, independent of 
the wavelength. For familiarization with these param- 
eters, numerical values of 7 and @ for a common sand- 


wich panel are given in the following. 


Description of panel 


Size: 1s & 1840. 
Skins: 0.01 clad 2024-T3 aluminum alloy 
Core: 1/S8-in. cell size, 0.007 3003-H19 aluminum 


foil, 0.375 thickness, density 3.1 Ibs./ft. 
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LOADED SANDWICH 


For this panel 


= 1.8 X 10-* (m* + n*) 
= 8.4 X 10—* (m? + n?*) 


2632 


[he corresponding fundamental frequencies are 


w/w, = 1.29 K 10° 
@/0. = 0.98 

_ w/a. = 1.31 & 10° 
w3/w, = 1.29 K 10? 


In general n is much greater than ¢. Notice that the 


expressions for w)/wo (Or @1/w,) and w»2/w) (OF we/w,) 





[Eqs. (16a) and (16b)] are symmetric in 7 and ¢. In 
other words, 7 and ¢ play the same part in affecting the 
frequency of Mode 1 and Mode 2. Since for sandwich 
panels in general » > ¢, it can be concluded that the 
effect of shearing deformation is much more important 
than the effect of rotational inertia; in fact, it is more 
important than in the case of homogeneous plates. 
Since 7 and ¢ are inversely proportional to the square 
of the wavelength, it can be seen that the effects of 
shearing deformation and rotational inertia are more 
important for shorter wavelength pulses. This fact 
is well known. 

The introduction of these parameters is very con- 
venient and of practical value. Suppose for a particu- 
lar panel we want to know whether the bending rigidity 
of the facings, say D,’, is important. We simply 
add Dy’ to Dy, to obtain a new value of 


Dy + Dy’ 
n = +e yD 


rT 24 wg 
(m* nB?) 
Do | PR 


and by substituting it into the frequency equations 
get a very good estimate of the frequencies for which 
the bending rigidity of the facings is included.* The 
same procedure can be applied to estimate the effect of 
the rotational inertia of the facings. 

In addition to y and @¢ there is another important 
parameter, K = E’t'/E"t", which is the ratio of the 
stiffness of the lower facing to the upper facing. A 
determines the position of the neutral surface of the 
panel, and takes care of the effect of the difference of 
thickness and material property between the upper 
facing and the lower facing. For a symmetrical panel, 
i.e., a panel with identical upper and lower facing, A = 
|. For a large class of practically important applica- 
tions, however, A ~ 1; for instance, due to aerodynamic 
heating, one face of the panel may be weakened con- 
siderably (E’ ¥ E"). 


Behavior of the Three Frequencies for Limiting 
Values of the Panel Parameters 7 and ¢ 

(a) n <1, @¢ < 1~—ie., the panel is very rigid in 

transverse shear and has very small rotational inertia. 


* This does not predict the appearance of additional modes of 


vibration. See reference 4. 
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In this case 


w/a, = 1 — (1/2)(n + @) 

we/we = [1 + (1/2)(n + o)]/V no 
ied. / 

s/t ~[1+ , n| V ne 


a = w. = (m?/L,?)(m? + n23?)V Dy p(l +f 


The result that w, and w; go to infinity as 7 and ¢ go to 
zero gives us another way of looking at the classical 
theory. Instead of saying that these two modes of 
vibration do not exist, it is more appropriate to say 
that their frequencies of vibration go to infinity when 
shearing deformation and rotational inertia are neg 
lected. 

(b) » < 1—1.e., the panel is very rigid in transverse 


shear: 

W | ] 2 

} ) | 4 ¢ Q 
l+@ 
We W — 5 (: es ; ) 
np zi+¢ 
l— yp y 

@®3/W- = 1 + ; n V nD 

When 7 = 0, 


GS &, V1 + ce) 


w, differs from the classical flexural frequency only by 
the factor 1/W1 + @. In general, unless for higher 
overtones (m,n very large), ¢ is a very small quantity, 
and w, can be very closely approximated by ow, => 
w.[1 — (@/2)]. It can be further inferred that the 
effect of rotational inertia is important only for higher 
overtones, or what is equivalent, for short-wavelength 
waves. This fact is well known. 

(c) ¢ < 1—1.e., the panel possesses very small rota 


tional inertia: 


- 


W3/We = y( + , 5a) no 


When ¢ = 0), 
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As in case (b), w, differs from the classical flexural fre- 
quency a». only by the factor | VI + ny. Since 7 is 
in general a much larger quantity than ¢, the effect of 
shearing deformation on the flexural frequency is also 
much larger than the effect of rotationalinertia. Again, 
it can be seen that the effect of shearing deformation is 
greater for higher order modes. 

(d) » > 1—.e., the panel is very soft in transverse 


shear: 


IIe 


[ oo Vite 


a - n+ ¢ 
wie =[1+5 15] —. 


W3/ W, =¥ ] + 
26 n(l — p) 
When 7» > &, 
oa, = 0 


w = w,/V¢ 
w3 = avi = 2) m 


In the present case, the flexural frequency w, vanishes 
asn— ©. This is to be expected since in our theory 
we have neglected the bending rigidity of the individual 
facings; therefore when the core no longer transmits 
transverse shear, the bending rigidity of the whole 
panel vanishes. The two predominantly shear modes, 


w2 and w;, do not vanish even though the core no longer 


(2 rol _ h? 24 92g? 2 1 
oe Boat i ™ n?B?) yore _ 


hye 


carries transverse shear. There is no contradiction in 
this case, because essentially in these two modes the 
two facings execute extensional vibrations which result 
in the shear of the panel as a whole. The following 
table summarizes the limiting behavior of these fre- 
quencies. The other limiting case ¢ ~ © is not con- 
sidered, since p,/p; does not become very large. 


TABLE 1 
wo @) we ; Ww; 
nO) 
We x 
o—> 0\ 
n—>OU x a: V1 + ti) x 
¢—~0 00 w V1 +7 ~_ = 
n— « O 0 we/a/¢ we/ (1 — n)/o 
we = (r?/Lz?)(m? + n2B2)~V 'D, pl +f 


Application to Homogeneous Plates 


It can be shown that the frequencies obtained in the 
previous section are applicable not only to sandwich 
plates but to homogeneous plates as well. For homo- 
geneous and isotropic plates, the physical parameters 
have the following values: 


D, = Eh®/12(1 — w2), K=1, f =0 


p,/pr = 1/3 within the assumption of first-order 
theory; 

G,. = (m®/12)G where 7/12 is the shear coefficient; 

n [2h?/L,?(1 — w)|(m? + n?B?) 

and @ = (7h?/12L,")(m? + n*6?) = [w?(1 — w)/24]n 


The three frequencies now become 


h? p = x w\ |? x h? : x 
1+ (m> + n’B") + - 3 D2 (m?> + n*3*) ¢ (21a) 


we\? 1 h? r ae ») gt 
at ae ) 1+ L (m? + n?B") | = T +. 





9 


z 


(2 = 1+ h?/L,? (m? + n°B?) 


Wo 


where wy = 


re 
mv 
he . a 9 r\ |? - — 
1 + — (m?+ ng? —+ + - (m? + n?B*)- = (21b) 
a l—u oe, f 


(21c) 


rVG/hp; (21d) 


With a change of notation, these frequency ratios can be shown to agree exactly with those obtained by Mindlin, 
Schacknow, and Deresiewicz.* It is interesting to note that these frequency ratios depend on only a single param- 
eter, (h?/L,")(m? + n°8?), besides Poisson’s ratio. Physically this is quite clear. Since the portion of the plate en- 
closed by nodal lines behaves exactly like a simply supported plate, the frequency of the overtone corresponding to 
this nodal pattern would be the same as the fundamental frequency of a plate with geometrical dimensions similar 


to the portion of the plate bounded by the nodal lines. 


Application to Sandwich Beams 


The results obtained in the previous sections are for sandwich panels. In the case of sandwich beams the third 
mode does not exist. Except for this difference, the results obtained for the frequencies and amplitudes of Mode |! 
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and Mode 2 for sandwich panels are exactly applicable to sandwich beams provided we interpret properly some of 
the parameters. All we have to do is replace the bending rigidity of the panel by the bending rigidity of the beam [let 


uw = Oin Dy = E’h'h*(1 + f)?/(1 + K)(1 — w?)] and let the wavelength in the y direction go to infinity (y8 > 0). 


Effects of Viscoelastic Core and Edge Compression 


In the previous sections, the application of our general theory to the flexural vibrations of a simply supported rec 
tangular panel was discussed in detail. For simplicity, the effects of viscoelastic core and edge compression were not 
considered. In this section we shall consider these effects in the manner of a first-order correction to the previous 
results. A more exact and more general treatment of these effects will be the subject of a later investigation. 

When the effect of viscoelastit core is considered, Eq. (lle), which gives the frequency and damping coefficient of 
Mode 1 and Mode 2, becomes 

(p/ao)* + (p/wo)*rao(1 + 6) + (p/we)?(1 + 9 + ) + (p/ao) TH0(nd) + nd = OV 
This is a quartic equation which could be solved by the standard method. The result, however, would be cumber- 


some. We shall find an approximate solution by utilizing the fact that tw) is small (small damping) so that we can 


expand (p/wo) as a power series in Twy—.€., 
(p/ao) = (y + 1w)/ao = (p/w) + (p/n) (two) + (p/w) (Two)? +... 


It can be shown that » is purely imaginary and gives us the frequencies obtained previously; p"" is real and nega- 
tive and gives the damping coefficient; p“ is again purely imaginary and gives the correction to frequency due to 
the existence of damping. 

Correct to 0[(7w))?] the damping coefficient is 


y* = (p/an) tay < O 
MUO Che) et j— Mode | 
= on a - 7 
4 V (1 ++)? — 4nd ge Mode 2 


{—Ta(n?¢/2) for. <1 (Mode 1) 
” or 
(—soo(1 + o)/2 @¢<1 (Mode 2) 


™ 


This shows that Mode 2 will be damped out much faster than Mode 1. This is to be expected since Mode 2 is pre- 


dominantly shear. The frequency is changed by the amount 
Pa" Po" Pa‘? /wo 16 (pa /wo)? + 1 + 9+ ¢)| + 3(1 + ) (pa /ws)? + no 


(rw)? = — (Tw)? Ee Ee a I 


w ; Pa® 2 [2 (w/a)? + (1 ++ 9)] 
( _ 1, Mode 1 
a 2, Mode 2 


For 7 <1, ¢ < 1 the change in the frequency of Mode 1 is negligible. 


The nondimensional frequency of Mode 2 is 


reduced by the amount (1/8) (zw). 
As to the amplitudes of vibration, it was already pointed out after Eq. (15b) that because of the different amounts 
of damping, the displacement components u, in the plane of the panel and the deflection u; are not exactly in phase 
or completely out of phase. 
Next we shall consider the effect of edge compression. 


frequency becomes 


{ 3 2 
() + (2) tao(l +o) + (2) (; +r7t+o- . ; cnn + (2) Two( nd) 
ay Wo Wo 


wo 


In this case the equation for the damping coefficient and 


x (1 — ve ) + 16 0 — Qmn) = 0 
n 


where a», is the ratio of the edge load to the corre- = ‘ — 
, : age : : Damping coefficient 
sponding buckling load. This quartic equation can be 
solved approximately by the power-series expansion Mode 1: proportional to (2 
method used above. The result for the case of 7 < 1 Mode 2: almost independent of a,,, 
and ¢ < 1—1i.e., for common sandwich panels in their 

Frequency 


fundamental modes or lower overtones—is summarized roe : 
Mode 1: proportional to Vi - Can 


— Amn) 


as follows: 
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Fic. 4. Normalized frequency of Mode 2, @)/a» (=w-/we) vs 
V¢ at different n/¢ for the cases (a) simply supported rectangular 
sandwich plate and (b) isotropic elastic core. 








0.2 04 


Mode 2: almost independent of a,,, 


Therefore the frequency of Mode | decreases with the 
increase of edge compression, a fact to be expected. 
The result also suggests a method of determining ex- 
perimentally the elastic buckling load of a sandwich 
panel by measuring the change of the frequency of Mode 
| with respect to the change of edge compression. 
Mode 2 is not much changed by the presence of edge 
compression, again because it is predominantly a shear 


deformation. 


Graphical Presentation of Results 


The frequencies for the three modes of the flexural 
vibration of an elastic and isotropic sandwich panel 
are computed, and the results are presented in Figs. 
3,4, and 5. The numerical values of the panel param- 
eter range from @ = 0) and n/¢@ = 0 to @ = 1.96 and 
n/d = 500. In each figure the frequency is normalized 
by dividing it by w = (1 + K)VG,. hp,, the fre- 
quency of pure shear vibration of an infinite panel. 

In Fig. 3 the normalized frequency of Mode 1, w/w», 
is plotted against the rotational inertia parameter Vo 
for different values of n/¢@. For a particular panel 
this ratio assumes a definite value, in particular 7/¢@ = 
24/[r?(1 — uw)] for a homogeneous plate. Therefore 
each curve of this figure represents a particular panel. 
As we move along that curve in the direction of in- 
creasing ¢, we obtain frequencies of higher overtones of 
that panel (m? + n°@? increases, while the panel param- 
eters remain unchanged). This is the reason that the 


ratio n/@ is used as a parameter instead of 7 itself 
When we want to see the effect of increasing (decreasing) 
the shearing deformation on the frequency of vibra- 


tion, we simply pick out a fixed abscissa 1/¢ and move 
from a curve of smaller (greater) 7/@ to a curve of 
greater (smaller) n/@. In order to see the effect of 


rotational inertia, we can follow a similar but slightly 
more tedious procedure. It should be kept in mind, 
however, that as we vary 7 and @¢, the characteristic 
frequency w = (1 + K)VG,. hp, also changes 
Therefore it may seem to be more significant to use 
w. = (2?/L,")(m? + n?2)V Dy (1 + f)p:, the classical 
flexural frequency, instead of wy, the shear frequency, 
as a basis of comparison. This is done in Fig.4. How- 
ever, the plot in. Fig. 3 does show certain advantages. 
It is seen that the curves for different values of the 
ratio n/@ show very small scattering. In fact for n/¢> 
50 the curves approach the straight line /¢ 
(w/w = Vv ¢) very rapidly. For designer’s use, this 
is very convenient. Also favoring the use of wy) as a 
basis of comparison is the fact that wy is fixed for a 
particular panel and independent of the wave pattern 
Fig. 4 is a plot of the normalized frequency of Mode 
2, we/wy against V ¢. The numerical values of the 
parameter »/¢ range from (0) to 500. Since it has been 
pointed out that a»/wy = w,/a, this figure also gives 
the ratio of the classical flexural frequency, w,, to the 
frequency of Mode 1, a. It can readily be seen from 


this figure that 
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Fic. 5. Normalized frequency of Mode 8, w3/wy vs. Wn for the 
cases (a) simply supported rectangular sandwich panel and (b) 
isotropic elastic core. 
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(a) Mode | 
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(c) Mode 3 
Fic. 6. The modal shapes of the three modes of vibration of a 


sandwich panel. 


(a) wy S we. 

(b) w/w, decreases with the increase of 7 and @. In 
other words, the greater the effect of rotational inertia 
and shearing deformation, and the higher the order of 
modes, the smaller is the frequency of Mode | in com- 
parison with the classical flexural frequency. 

(c) For large values of ¢, a /w, = | V nd. 

It can be seen by comparing the scale for w;/wy and 
w/w in Figs. 3 and 4 that w, can be much smaller 
than we. It should also be pointed out that w,/w. = 
w;/wo, or the ratio of the classical flexural frequency to 
the frequency of Mode 2, is also given by the ordinate 
in Fig. 3. Again this is very convenient for the rapid 
estimation of the frequency of Mode 2. 

Fig. 5 shows the normalized frequency of Mode 3, 
w;/w. In this figure w/a) is plotted against V 7 in- 


stead of V 4, because 


m 


| 
| — 
ws/wo Vit 7 


depends on 7 only. For large values of 7 the curve 
approaches the straight line w;/a) = V(l — p)/2 Vn 
asymptotically. 

It can also be seen from Figs. 3, 4, and 5 that in 
contrast with the classical flexural frequency, over- 
tones of a, we, and w; are spaced closer than the order 
of (m? + n°B*). In fact, for very high overtones 
(m,n >> 1) they are approximately proportional to 
V ‘m? + np". 

The modal shapes of these three modes of vibration 
are shown in Fig. 6. The distinction in the deforma- 
tion associated with each mode is evident, as has been 
discussed in a previous section 

The ratios of the amplitudes of vibration of Mode 1 
and Mode 2 are also computed. The ratios of the dis- 
placement components in the y and g directions, 4 ,,,“! 
Amn = Bryr™/Bn,™, are the same and are shown in 


mn 
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Fig. 7. For large values of ¢, this ratio approaches the 
asymptotic value 7/¢. The ratios of the deflection in 
the z direction, C,,,°"/C, , are shown in Fig. 8 
These figures seem to indicate that the amplitude of 
Mode 2 is negligible compared with the amplitude of 


Mode 1. 
to a unit step function, and.it should be understood 


However, these are the amplitude responses 


that Mode 2 cannot be disregarded under general load 


ing conditions. For loading conditions in general, the 
coefficients of the Fourier expansion of the load form 

rapidly decreasing sequence. Therefore, for sustained 
high-frequency excitation the fundamental vibration of 
Mode 2 will be more important than the higher over 


tones of Mode 1. 


Conclusion 


(a) In this paper the flexural motion of a rectangular 
sandwich panel with an orthotropic viscoelastic core 
is treated. It is found that three families of different 
modes of vibration are possible. The first mode cor 
responds to the classical flexural vibration. The 
motion of this mode is dominated by flexural deforma- 
tion coupled with a small amount of shearing deforma- 
tion. The second mode is dominated by thickness 
shear deformation with some extent of flexural deforma- 
tion. The frequency of the second mode is much 
higher than that of the first mode. This second mode 
may become very important under sustained high- 
frequency excitation, inducing large shear stresses in 
the core. The third mode is pure shear. For simply 
supported isotropic panels as treated in the example in 
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the preceding sections, this mode is not excited by 
loads perpendicular to the panel surface. For other 
boundary conditions, and for orthotropic panels, in 
general all three modes are coupled together and hence 
ail are excited by loads perpendicular to the surface 
of the panel. 

(b) Two fundamental-parameters, 7 and ¢, are de- 
fined. They represent the relative importance of rota- 
tional inertia and shear deformation, respectively. 
In terms of and @¢, the three modes are expressed in 
simple forms which are convenient for engineering 
applications. When 7, ¢ — 0 the second and the 
third modes disappear, and the first mode reduces to 
that of the classical theory with shear deformation and 
rotational inertia neglected. The greater the param- 
eters 7 and @¢ the smaller is the frequency of the first 
mode as compared with the classical frequency. In 
conventional sandwich panel, n/¢ > 1. 

(c) The theory is applicable to unsymmetrical sand- 
wich panels. The difference in the thickness and the 
material of the two facings of the panel is very con- 
veniently taken into account by the parameter K, 
which is the stiffness ratio of the two facings. K = ] 
for identical facings. 

(d) The effect of edge compression and viscous shear 
in the core are represented by the parameters a,,, and 
t, respectively. As should be expected, the frequency 
of the first mode decreases with increase in the edge 
compression. The damping coefficient is greater for 
the second mode than for the first because the second 
mode is essentially a shear motion. 


Remark 


After the completion of this paper, some papers by 
Yu*° on sandwich plates were referred to the authors. 
In these papers a one-dimensional theory of an infinite 
symmetrical sandwich plate is developed from the 
theory of elasticity, taking into account very general 
core and facing properties. The theory has also been 
generalized to the two-dimensional case.'! It is be- 
lieved that the work of Yu is very important in estab- 
lishing the limit of validity of the present theory, which 
is more restrictive but more amenable to engineering 
applications. 
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Some Considerations on the Laminar Stability 
of Time-Dependent Basic Flows' 


S. F 


. SHEN* 


U.S. Naval Ordnance Laboratory 


Summary 


As a stability criterion for infinitesimal disturbances in an 
incompressible, parallel but time-dependent basic flow, it is pro- 
posed to introduce the concept of “‘momentary stability,’’ which 
is said to prevail at the instant if the kinetic energy of the dis- 
turbances, as a fraction of the kinetic energy of the basic flow, 
tends to decrease. The significance of such a criterion is briefly 
discussed. For special time-dependent basic flows which are de- 
scribed by similar velocity profiles at all times (except for changes 
in amplitude), im the inviscid limit only a change of the time scale 
is needed to reduce the solution essentially to that for the steady 
The disturbances may be of either the transverse-wave 
The result indicates a very 


case. 
or the longitudinal-vortices type. 
strong destabilizing influence of deceleration, which is likely to 
overshadow that of the velocity profile under normal circum- 
stances. The observations of Fales (on flow over the flat plate) 
and Coles (on flow between rotating cylinders) are believed to be 
largely due to the deceleration. At finite Reynolds numbers, 
the usual procedure of calculating the stability solution on the 
basis of the instantaneous profile is further shown to be valid 
only for extremely slow acceleration or deceleration. Even when 
the solution is acceptable, the condition c; = O for neutral sta- 
bility may not be used without reservation. To calculate 
momentary stability properly, a procedure for a slowly varying 


but more general profile is also described. 


Symbols 

( = wave speed of the disturbance 

ae = real and imaginary parts, respectively, of the wave 
speed 

t = time 

to = time constant of the changing basic flow 

u,v = disturbance velocity components in the x- and y- 
directions, respectively 

x,y = Cartesian coordinates along and normal to the basic 
flow, respectively 

y = thecritical point where U* = « 

E = kinetic energy ratio, Eq. (5) 

Eo, E; = defined in Eq. (42) 

G, = ‘growth rate”’ defined as (1/2E)(dE/dt) 

Re = Reynolds number 

Rec, = minimum critical Reynolds number 

r = time dependence of the special velocity profile, Eq. 
(6) 

T, = time dependence of the velocity profile, Eq. (25) 

U = velocity profile of the basic flow 

i = maximum velocity of the basic flow 

U*(y) space dependence of the special velocity profile, 
Eq. (6) 

U,*’ = velocity gradient of U*(y) at y = y, 

a = wave number of the disturbance 
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B time constant parameter for certain special basic 
flows 

6 = thickness of the basic flow in y-direction 

€ = small parameter for slowly changing profile, Eq 
(25) 

ny = wavelength of the disturbance 

® = disturbance mode defined by Eq. (17) 

®, = the eigenfunction for the disturbance 

Ps = defined in Eq. (31) 

Vv = disturbance stream function 

Ws = defined in Eq. (29) 

¢ = defined by Eq. (19) for Eqs. (20) and (21), or by 
Eq. (36) for Eqs. (37)-(39) 

T = variable of stretched time scale, Eq. (8) 

Superscripts 


(0), (1) refer to terms in the expansion Eq. (26 


Introduction 


| IS NOW GENERALLY ACCEPTED that the theory of 
laminar instability provides the explanation for the 
initiation of large disturbances which may lead to the 
transition to turbulent flow when no other causes are 
Heretofore the basic ‘‘mean’’ flow has been 
Where the basic 


present. 
assumed to be steady and parallel. 
flow is neither of these, calculations have been made 
by applying the theory to the instantaneous local ve- 
locity distributions 
or quasi-parallel treatment. 
to assume that such treatment probably gives the right 


i.e., resorting to a quasi-steady 
It is reasonable, of course, 


answer when the changes are indeed sufficiently slow 
with respect to the variable concerned. But likewise 
the range of validity of this procedure in actual appli- 
cations—i.e., the maximum permissible time- or space 
rate of change Thus, G. I. 
Taylor! voiced an early criticism regarding the correla- 
tion of the laminar stability theory of parallel flows 
with the transition of the boundary layer, for which 
there is always a growth in the downstream direction, 
This objection has been counter- 


is by no means clear. 


however small. 
balanced, at least from the experimental side, by the 
results of Schubauer and Skramstad? which confirm in 
detail the theoretical predictions, based upon the local 
flow, in the case of the boundary layer over a flat 
plate. 

On the question of unsteady ‘‘mean”’ flow, no sys- 
tematic investigation seems to exist in the literature. 
In estimating the transition Reynolds number of the 
flow between rotating cylinders, Schlichting’ applied 
the laminar stability theory for transverse-wave (Toll- 
mien-type) disturbances to the instantaneous profiles 
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taken when the outer cylinder is suddenly set into 
motion. The minimum critical Reynolds number 
based upon the displacement thickness turned out to 
be about one third of the experimental transition Reyn- 
olds number. Since one can usually rely on the sta- 
bility theory only for qualitative trends of transition, 
Schlichting regarded the agreement to be as close as 
might be expected. However, the probable effect of 
unsteadiness of the mean flow was left completely un- 
touched. More recently, two experimental observations 
of the generation of the vortices of the transverse-wave 
type in decelerating flow have come forward. Fales* 
demonstrated this phenomenon by gradually stopping 
a moving flat plate in a water tank. The slowing 
down of the laminar boundary layer caused distinct 
transverse vortices to appear, made visible by coloring 
the fluid with dye. Coles,® in a series of observations 
of the Taylor stability of the flow between rotating 
cylinders, included a start-stop operation of the outer 
cylinder alone. As in the boundary-layer case, after 
the outer cylinder was stopped, vortices with axes par- 
allel to the circular cylinders were shown to form, to 
parade and to grow in intensity before eventually burst- 
ing into turbulence. The only theoretical attempt to 
explain this phenomenon of vortex creation by decelera- 
tion, to the author’s knowledge, was that due to von 
Karman and Lin® on Fales’ experiment. They noted 
that the instantaneous velocity profile of the boundary 
layer after the plate was stopped should resemble that 
of a free-convection boundary layer. The presence of 
a point of inflection of the velocity profile suggests that 
the profile should be quite unstable with respect to 
the transverse-wave disturbances. To support the 
explanation, a calculation of the critical Reynolds 
number of such a velocity profile was carried out by 
use of the standard theory. Under conditions typical 
of Fales’ experiment, the comparison was shown to be 
qualitatively satisfactory. Still, the emphasis was on 
the shape of the instantaneous velocity profile alone, 
so that these authors, like Schlichting, implied the 
validity of the quasi-steady approach. To be stressed 
is the fact that thus far no detailed comparisons have 
been made of the behaviors of the disturbances in an 
unsteady basic flow with those in a steady basic flow, 
such as were provided by Schubauer and Skramstad 
in the case of slightly nonparallel (boundary-layer) 
flow over a flat plate. A more critical examination of 
the disturbances in an unsteady basic flow appears to 
be of fundamental interest, even if from the theoretical 
viewpoint only. 

The present paper represents a step in this direction. 
In the first place, the criterion for laminar stability of 
an unsteady basic flow requires some modification from 
that of the usual steady basic flow. It is proposed to 
adopt the concept of ‘“‘momentary stability” in dealing 
with time-dependent basic flows; such stability is said 
to prevail when the disturbances tend to become less 
disturbing at a giveninstant. The concept of ‘‘momen- 
tary stability,’ of course, is a weakened version of the 
usual concept of stability of a steady mean flow, which 


requires the distu bances to vanish as t > ©. The 
more restrictive definition is necessary since within 
the period of interest the unsteady basic flow never 
attains a condition of f > ©. Furthermore, as the 
basic flow is changing with time, a disturbance tends 
to be unstable only if it grows more rapidly than the 
basic flow. It is proposed to measure the disturbance 
on the basis of the ratio of the kinetic energy of the 
perturbations to that of the mean flow, both averaged 
throughout the boundary-layer thickness and over 4 
wavelength. 

The physical model to be investigated is an incom- 
pressible, parallel but unsteady basic flow, upon which 
are superposed two-dimensional infinitesimal perturba- 
tions of the transverse-wave type.t The general dis- 
cussions are centered around a special ‘‘similar’’ basic 
flow which preserves the velocity profile except for a 
time-dependent magnitude. The advantage of the 
special flow is that it permits a transformation of 
the time scale in the inviscid limit, so that the per- 
turbation equation may be reduced to the same form 
as for a steady basic flow. Deceleration is found to 
have an almost overwhelming destabilizing influence 
regardless of the shape of the velocity profile. Fales’ 
and Coles’ experiments should perhaps be explained 
from this viewpoint rather than on the basis of any 
quasi-steady consideration. 

The quasi-steady method of calculation turns out 
to be acceptable only for extremely small acceleration 
or deceleration—the time constant of the basic flow 
must be much greater than Re’ times the character- 
istic time of the disturbances. Even then the quasi- 
steady ‘‘neutral curve’ cannot be accepted at face 
value with regard to its significance for stability. 

A formulation of the procedure for calculating the 
“momentary stability’? boundary of slowly changing 
but otherwise completely arbitrary basic flows is out- 
lined briefly. Essentially, the eigenvalue problem is 
taken as identical with the quasi-steady approach, but 
a time-dependent part of the ‘‘mode’’ is determined. 
With the time-dependent mode, the momentary sta- 
bility or instability of, say, points corresponding to 
the quasi-steady neutral curve may be determined. 

It should, perhaps, be noted that investigations in 
this paper are directed only toward the types of dis- 
turbances normally considered for steady basic flows. 
Whether other types of disturbances may not also be 
of importance for unsteady basic flows is a question that 


is left completely open. 


The Momentary Stability Concept and the 
Growth Rate Criterion 

Let us first recall the classical stability theory for a 
steady, two-dimensional, parallel basic flow. The 
fundamental linearized differential equation for two- 
dimensional disturbances superposed on the steady 
velocity profile U(y) is 

+ Needless to say, little basic difference should be encountered 
in dealing with disturbances of the longitudinal-vortex (‘‘Taylor- 


Gortler’’) type. 
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Here WV is the disturbance stream function such that 
the disturbance velocity components (“, v) in the x- 
direction (downstream) and the y-direction (normal to 
the flow) are given by 


u = OW/Od7, v= —OwW/Ox 


The boundary conditions on u and v are homogeneous 
and apply at the extreme values of y. Thus, the prob- 
lem is defined independently of the x- and ¢-coordinates, 
and one is led to assume the elementary solutions 


ta(x—ct) 
¢ 


Vv = P(y) (2) 


with the constant a real and the constant c possibly 
complex. The homogeneous boundary conditions now 
are prescribed for ®(y), so that an eigenvalue problem 
results. The equation for ® is the well-known Orr- 


Sommerfeld equation” 


(b’’ — a®’bh)(U — c) - U"@ 


(1 laRe)(e'”’ + a'® | (3) 


2arp’’ 


where Ke is the Reynolds number, the reference ve- 
locity and the reference length being both taken as 
Reynolds number Ke and wave 
The disturbance 


unity. For given 
number a, one gets the eigenvalue c. 
will be either amplified without bound as the time 
{— ©, or eventually damped out, depending on the 
sign of the imaginary part of the eigenvalue c. 

The above résumé of the usual theory for steady basic 
flow has been presented in order to draw contrast to the 
present case of interest. For the unsteady parallel 
basic flow, the fundamental equation, Eq. (1), remains 
unchanged even though U’ = Uy, t). Furthermore, 
let the boundary conditions remain homogeneous and 
prescribed at the extreme points of y. While it is still 
reasonable to assume solutions which are periodic in x, 
the time-dependence becomes a different matter. In 
this case the assumed elementary solution, analogous 


to Eq. (2), is given by: 
V e'™* Wy, t) (4) 


so that uw = e'™ (O%,/Oy) 


wr 
v => —lae VY; 


For given Ke and a, an eigenfunction V(y, ¢, a, Re) is 
defined. It is necessary, however, to determine which 
values of a and Re represent a “‘stable’’ solution. 

One must now reflect on what is meant by “‘stable”’ 
in discussing the possible disturbances represented by a 
given solution W;(y, f, a, Re). It is obvious that the 
asymptotic behavior of the disturbances as tf > © is 
no longer meaningful. In the first place, in considering 
unsteady basic flow it is already implied that behavior 
at finite time is of primary interest. In fact, the un- 
steady basic flow U(y, t) may be described only for 
finite times. Even if the velocity does have a limiting 
profile U(y, ©), the analysis for behavior as t ~ © 
would reduce to an examination of the profile U(y, ©) 
by the usual formulation for the steady basic flows. 


F 
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The result would be significant, however, only if the 
prescribed unsteady basic flow is maintained in spite 
of all the disturbances at all times prior to the estab- 
lishment of the limiting flow L(y, ~). Consequently, 
one must have some idea of how the basic flow is being 
distorted by the given disturbance W,(y, f, a, Re) at 
every instant. As a matter of fact, for unsteady basic 
flows, the tendency toward stability or instability at 
any moment is all that can be stated definitely. 

As the ultimate goal of the laminar stability theory 
lies in its possible correlation with the laminar-turbulent 
transition phenomenon, it is recognized that logically 
the increase or decrease of a disturbance should be 
measured only by a comparison with the basic flow 
itself. If the disturbance is decreasing but the basic 
flow is decreasing at a faster rate, the disturbance rela- 
tive to the basic flow would appear amplified at a later 
instant. 
the basic flow increases even more rapidly, then the 
From such 


Conversely, if the disturbance increases but 


disturbance would appear to be decaying. 
elementary considerations a yardstick for measuring 
the disturbances immediately suggests itself.* First 
form the ratio of the instantaneous kinetic energy of the 
disturbance to that of the basic flow, both averaged 
over one wavelength in the x-direction and throughout 
the thickness in the y-direction. The time 
change of this ratio then indicates whether the relative 


rate of 


kinetic energy of the disturbance is increasing or de- 


creasing at any instant. The energy ratio is given by 


5 fh 75 
[ | (u® + v*)dxdy | Udy (5) 
0/70 Ji 


Then for convenience one may define a ‘‘growth rate”’ 


E(t) = 


« 


G, = (1/2E)(dE/dt) 


“momentary sta- 
lf G, > 0, 
relatively at the 


and use it as the criterion for the 


bility or instability’’ of the basic flow 
the disturbances are 


moment, and therefore the flow may be said to be 


increasing, 


momentarily instable, whereas if G, < 0 the flow may 
be said to be momentarily stable. The term stability 
is thus used in a special sense, reflecting only the tend- 
ency of the basic flow to become more and more, or 
less and less, distorted. For the steady case, how- 
ever, the “growth rate’’ G, reduces to the amplifica- 
tion rate of the disturbances in the usual theory. 

What must be emphasized here is that the momen- 
tary stability defined on the basis of the growth rate 
carries a more ambiguous significance in its probable 
correlation with transition than does stability as defined 
in the usual sense (instability for unbounded disturbance 
as t > o) for the steady basic flow. For the latter, 
although the mechanism of transition is not yet fully 
clear, it is quite evident that conditions will change 
as the disturbance keeps on growing, and that even- 


* The same yardstick has been used previously by Lin® in 
discussing the stability of the decaying turbulence spectrum 
In a private communication to the author, Prof. Lin suggested 
that others might have used it under other circumstances. It 
certainly is a very natural, intuitive choice 
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tually turbulent effects must appear. However, for 
unsteady basic flow, the ‘‘growth rate’’ indicates no 
more than a tendency. To illustrate, suppose that 
the ‘growth rate’ is first positive (implying unstable 
disturbances, according to the definition) and then 
negative. As the disturbance becomes more and 
more pronounced during the time interval of “‘insta- 
bility,”’ does this mean that turbulence will set in dur- 
ing this interval? If turbulence sets in, it is not cer- 
tain that the flow may become laminar again at some 
later time, even if the “‘growth rate’’ calculated on 
the basis of the original laminar basic flow is, say, 
moderately negative. If turbulence does not set in, 
then the flow remains laminar in spite of the unstable 
period, since the disturbance will again be smoothed 
out in the later stable period. Therefore, in such 
cases the significance of the “unstable” period is rather 
uncertain. The question obviously is identical, in 
a way, with the difficulty of predicting the transition 
point for a steady basic flow, and can only be answered 
when the mechanism of transition is more clearly under- 
stood. In fact, a quantitative criterion of transition 
will be needed. When a theory to predict the transi- 
tion point for steady basic flow is established, then 
and only then will it be possible to calculate the transi- 
tion for unsteady basic flow according to the concept 
of ‘‘momentary stability’? as proposed here. Other- 
wise one cannot say as much about the unsteady case 
as about the steady case. Whether transition will or 
will not follow at a later time can be answered only in 
special cases by examining the ‘“‘momentary stability’ 
of the basic flow. However, on the basis of the (mo- 
mentary) stability criterion as discussed above, it will 
at least be possible to interpret the solution of Eq. (4) 
(assuming that a solution can be found) for every given 
a and Re. The procedure for solution is more tedious 
than for the steady case in that one must calculate not 
only the eigenvalue but also the eigenfunction itself, 
in order to evaluate the ‘‘growth rate’ from Eq. (5). 
Then, if the basic flow is ‘“‘momentarily unstable’ for 
given a, Re, and for all ¢ > t, the basic flow will be 
“unstable” for t > t, although transition will have to 
come at an even later time. On the other hand, if for 
Re < R.,, G, is less than zero for all a > O and all ¢, 
then Re,, will be the minimum critical Reynolds number 
for the unsteady U(y, ¢) concerned, in the same sense 
as for a steady flow. These are two examples for which 
at least the sufficient condition for stability can be 
written down by adopting the viewpoint of this paper. 

Even before actually determining the solution, one 
is led to suspect that accelerating flows very probably 
would have higher critical Reynolds numbers than those 
estimated on the basis of the instantaneous velocity 
distribution, and that the reverse would be true for 
decelerating flows. At least for nearly steady flows, 
when the time constant of the basic flow is much 
longer than the period of the wave disturbance, it seems 
permissible to calculate the disturbance solution on the 
assumption that the instantaneous basic flow will per- 
sist forever. Such was, of course, the procedure used 


by Schlichting*® in his calculation of the critical Reyn- 
olds number for instantaneous profiles of the flow be- 
tween rotating cylinders following the sudden start of 
the outer cylinder. In this case, the reference velocity 
(that at the outside cylinder wall) is constant, and the 
disturbance amplitude and hence its kinetic energy 
would temporarily remain constant along the neutral 
curve determined by the quasi-steady procedure, while 
the kinetic energy of the basic flow evidently increases 
with time. According to the viewpoint of the present 
paper, these quasi-steady neutral disturbances must 
actually be regarded as decaying relative to the acceler- 
ating basic flow. In fact the argument is even stronger 
than suggested by the reasoning applied to Schlich- 
ting’s example. To clarify the point, consider the 
special basic flow 

U = T(t/to) U*(y) (6) 


where é) is the time constant. The velocity profiles are 
now identical at different times except for the scale 
factor 7. By the quasi-steady procedure, the instan- 
taneous profile would be 7U*(y) and the neutral dis- 
turbances would have amplitudes proportional to 7. 
It would then appear that the energy ratio E, Eq. (5), 
would be independent of time along the quasi-steady 
neutral curve. The latter would seem thereby to 
satisfy the condition for neutral momentary stability 
regardless of whether 7’ represents an accelerating or 
decelerating flow. The following closer examination 
shows, however, that such a conclusion is erroneous 
the accelerating flow would actually remain stable 
for the quasi-steady neutral disturbances. 


Analysis of a Special Unsteady Flow 

For simplicity, let us first analyze the special un- 
steady basic flow, Eq. (6). Following the classical 
theory, we may first adopt the ‘‘inviscid’’ approximation 
to obtain some information regarding the behavior of 
the disturbances. Such information is valid as Re—> © 
and often can guide the treatment for finite Reynolds 
number. To determine the ‘“‘viscous’’ solution at 
finite Reynolds number, however, is rather difficult 
for the type of profile considered. Nevertheless, the 
limitations of the quasi-steady approach can be checked 
without much complication. 


(1) The Inviscid Approximation 

Consider now the abbreviated equation obtained by 
dropping the viscous terms from Eq. (1): 
0?U ow 


a) 
" Vv — = (7) 


fe) . 
Pr rs dy? Ox 


fe) 
where U = U(y,t). For the special flow 
U = T(t/to)U*(y) 
the mathematical problem thus defined can be reduced 
to the classical inviscid equation by a transformation 
which is essentially a stretching of the time variable. 
That is, if 
r= bf T(t/to) d(t/to) (8) 
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then WV satisfies 


ve + UG) — ve = oP" nd =0 (9) 
2 ) ay a ) = 
T ‘9 Ox ¥) Ox } 


) 
an equation which is identical in appearance to the 
classical steady case, except that 7 replaces ¢. One 
can immediately borrow the classical solution for U*(y) 
and the corresponding boundary conditions, and write 


Vv = e'* *—Cr) ?(y) (10) 


Therefore the eigenvalue c for a given a as Re > ©@ is 
the same as that in the classical case with the distribu- 
tion U*(y).F 

In this special type of the basic flow, the “‘growth 
rate’’ may be defined in a simpler manner, since the 
time dependence appears only as a separate factor, in 
both the functions U and ¥. Thus one may take 
the energy ratio equivalent to Eq. (5) as 


E = const. X (e*’ /T)? (11) 
and the growth rate is given by 
G,; = ac;l aad (7* 7) ty) (12) 


7’ being the derivative of 7 with respect to t/fy. Con- 
sequently, G, < Oif 
ac; < (1"/T*)(1/to) (13) 


It is seen that for truly ‘“‘neutral’’ disturbances, deter- 
mined by G, = 0, ac; will indeed be proportional to 7” 
(which is positive for accelerating flow and negative 
for decelerating flow) and inversely proportional to the 
time constant of the basic flow. (It is trivial to show 
that the same time-scale stretching, Eq. (8), works 
also for the inviscid case in the stability of disturb- 
ances of the Taylor-Gortler type for flows over a curved 
surface.) Since for the basic flow [Eq. (6)] 


(1/U)(OU/ot) = (T°/T)(1/to) (14) 

a more explicit form of Eq. (13) is 
a ¢, < (1/U)(0U/dt) (15) 
where ¢; = c;J. The quantity ac, clearly represents 


the amplification rate as computed by the quasi-steady 
procedure, the factor 7 entering here because of the 
change of the reference velocity with time. In other 
words, for the special profile Eq. (6), one may in the 
inviscid limit calculate the eigenvalue c for a given a 
by the quasi-steady procedure, but the “‘momentary 
stability’ is mot given by the neutral disturbances so 
obtained. The acceleration or deceleration has a very 
profound effect. With accelerating flow, c; can be 
positive and the disturbance will remain ‘“‘momentarily 
stable.’’ With decelerating flow, c; has to be suffi- 
ciently negative to maintain such stability. 

The experience with available computed results for 
steady parallel flows suggests that generally the quan- 
tity ac, in the unstable region of the stability diagram 

+ Note that the quasi-steady procedure is obviously inappli- 
cable in this special case, since, by Eq. (10), the disturbance 


amplitudes are not proportional to T. 
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Fic.1. Typical stability diagram for steady basic flow with pres 
sure gradient. 


is rather small, say, of the order of 10 Thus, as a 
rule of thumb, when the time constant of the changing 
basic flow is much shorter than 10° X (reference 
length/reference velocity of the basic flow), an acceler- 
ating flow will probably be always momentarily stable 
and a decelerating flow (if the Reynolds number is not 
too small) always momentarily unstable—whatever 
the velocity distribution may be. 

To illustrate this point further, a typical stability 
diagram including contours of constant amplification 
rates, as calculated by Pretsch,* is presented as Fig. 1. 
The velocity distribution corresponds to one of the 
Falkner-Skan profiles with adverse pressure gradient. 
Suppose now in an unsteady basic flow one has the 
same spatial velocity distribution but a time dependence 


of 


The condition for neutral momentary stability, in the 
inviscid limit, is then 


: —pt 
ac; = Be '' 


An accelerating flow will have 8 > 0, and the time con- 
stant f will be simply 1/8. An examination of the 
asymptotic behavior (at large Reynolds numbers) in 
Fig. 1 shows that up to the time @/ ~ 0(1) no disturb- 
ance is unstable enough to be disruptive if 8 > 0.013. 
Up to such a time, therefore the basic flow would be 
completely stable.| The case of decelerating flows may 
be discussed in a similar fashion; the conclusion is 
naturally in the opposite direction. Although these 
results are derived strictly only for the special flow of 
Eq. (6), the criterion in the form of Eq. (15) might be 
of qualitative significance even for more general flows. 

A corollary may also be drawn regarding the pro- 


7 Of course, inviscid consideration will not determine whether 
at finite Reynolds numbers there might not be an imbedded in- 
stable region, according to the concept of momentary stability 
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file shape U/*(y): The existence of a point of inflection 
in U’*(y) is much less significant for the stability, at 
infinite Reynolds number, of accelerating or decelerat- 
ing flows. In the accelerating case, the inflection 
point is necessary for instability but not sufficient, 
whereas in the deceleration case it is sufficient for in- 
stability but not necessary. 


(2) A Critical Examination of the Quasi-Steady Approach 


To determine the minimum critical Reynolds num- 
ber according to the concept of ‘‘momentary stability,” 
the full Eq. (1), including the kinematic viscosity 
It seems clear right from 


v, must be considered. 


the beginning that the special case 
U = T(t/te)U*(y) 


which permits a stretching of the time scale to reduce 
the inviscid differential equation to the classical form, 
no longer has the same advantage at finite Reynolds 
numbers. However, it is convenient to use it as an 
example to demonstrate the restrictions of the quasi- 
steady approach. Consider the complete viscous form 
of Eq. (9): 
(0/07r)V°W + U*(y)(0/dOx)V°h — U*''(0W/dx) = 
(v/T)V°V°" (16) 


One sees qualitatively, besides the time-stretching 
(with 7 replacing ¢), that the effective Reynolds number 
with 

Yer, = ¥/T 

Now let us follow the steps of a 
Guided by the inviscid solu- 


will vary with time. 
quasi-steady approach. 
tions for Eq. (9), we may postulate the solution of Eq. 


(16) in the form: 


Y = o%-*) By, 7) (17) 


The equation for ® reads 


(0/d7)V°@ + ia[(U* — Vb — U*"b| = 
(v/T)V°V2@ (18) 


with V? = (07/d0y") — a. Now for the viscous solu- 
tions, we may use the same classical arguments, focus- 
ing attention near the critical layer y = y. where L’* — 
c = 0 (see reference 7), thus Eq. (18) is simplified to 
(0/O07r)V2® + ta(U* — c)V°b & (v/7)V°V72@ 

In the neighborhood of the critical layer, however, 
a transformation may be used to magnify the y-direc- 
tion in the standard manner: 


6 = (vy — Me)(aU.*'/ vert.) (19) 


Hence the differential equation becomes 
ver (aU .*’)~ 10/7) 02h + iV? = 

(07/0¢?)V°@ (20) 
If the term in 0/07 on the left-hand side is omitted, the 


solutions are known to be 


Vb = AEM) © [2/3(ig)? | + O(ven /aUo*’)'* (21) 
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as long as c; is small. From Eq. (21) an estimate of 


(0/07r)V°*@ may be made: 


a Oln ¢ O 
V*) = E # ve | 


Or Or OK 
But 
lial Ll df az 1 
cs de Tis ' at T? 


so that the first term of Eq. (20) is negligible only if 


;7* 7 


me < v’* (al*’) (22) 
31°" t 


Again, since in this special case 
1 ol’ 7 
U’ of T to 


Eq. (22) may be rewritten in dimensional form as 


1 oU ( ou, - 
: < bv ra (25) 
LU’ ot \™ oy 
Now let 
U, = maximum value of the instantaneous velocity 
distribution 
Re = Reynolds number based upon the instantane- 


ous value L’,, and the thickness 6 


Then a more convenient version of Eq. (23) is 


— 1 


(Unto 6) . < (ad) Re (24) 


as long as (OU’./Oy) ~ O(U,,/6). Physically, the left- 
hand side represents the ratio of the characteristic 
time 6/U,, for the disturbances to the characteristic 
time f) of the basic flow. If we calculate the neutral 
curve by the quasi-steady procedure, we find that 
ad ~ O(1) near the minimum critical Reynolds number, 
and that the ratio 6/U’,,t) has to be less than O(/Re™ “°) 
to justify the solution. In addition, it may be recalled 
that ordinarily the amplification rate ac; in the un- 
stable region of the quasi-steady stability diagram is 
very roughly speaking, smaller than ((Re~“*) near 
the minimum critical Reynolds number (e.g., see Fig. 
1). The disturbances calculated by the quasi-steady 
procedure thus would amplify or decay at a rate prob- 
ably comparable to the rate of change of the basic 
flow that satisfies Eq. (23). From this consideration, 
even if Eq. (23) is satisfied one still cannot justify the 
use of the quasi-steady procedure alone to predict the 
minimum critical Reynolds number for momentary 
stability. The remains the criterion 
for such purposes. 
from the quasi-steady procedure seems to be of very 


“srowth rate” 
Even the neutral curve derived 


limited significance in dealing with unsteady basic 


flows. 


Analysis of the Momentary Stability of Slowly 
Varying Profiles 
profile U 


While illuminating, the special 


T(t/to)U*(y) is too restrictive for most actual applica- 
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tions. Furthermore, the stretching of the time scale, 
by Eq. (8), is no longer helpful when finite Reynolds 
numbers must be considered. Therefore, a slowly 
varying general profile will now be examined: 


U(y, 1) = Uy) + eTi(t/t)U™ (y) +... (28) 


¢ being a small constant parameter. The solution WV 
will be developed in ascending powers of e, 


¥ = YOR, 9, 2) + Or, 9,2 +... (26) 


It is then necessary to solve for the perturbation Vv" 
for given a, Re, noting that Y should be the classical 
solution for U'(y). Upon substituting Eq. (25) into 
Eq. (1) and equating terms of order e, there results 


a) v yi I ov"! 
vw » — yyy — Yor - 
of * Ox Ox 
1 7: [ 71 Q [ l 0 | 
VV = — v9 — po" — ge 
‘ad Ox Ox 
(27) 
A particular 7) (t/t) convenient for treatment is 
7) =e (28) 
for which it is possible to take 
VO = MW, (x, y, £) (29) 


The equation for V4") becomes, upon substitution of 
Eqs. (28) and (29) into Eq. (27), 


( >) 1, 4 Yo 2 gay — yor Me 
3+ ap (0 yo — Yor" — 
dad aie ai a ax 
1 U 1 0 af U ] 0 
IT 27 2! —_> — ap) (1)// wt 30 
vv°V-N oo + i. } (30) 
Let YO = ht*—% GO (y) 


Then one may take 


WV, = o* x t) Py l (y) (31) 


so that Eq. (30) reduces to 


B [ @ 
(0) - , —a-|@®@ <a 
[ ( ’ | BS | ™ 
— | Py , -— 


= d* + 47 
— U? ly? —a? lo +4 UO pl (32) 
dy? 


One recognizes the left-hand side of Eq. (32) as identi- 


cal with the equation which determines the classical 
for U(y), except for the term 7(6/a) 
In addition, there are 


UNG, 0 


solution &° 
added to the wave velocity c. 
the inhomogeneous terms on the right-hand side de- 
pending on U and 6. Let the boundary condi- 
tions on &, ') remain homogeneous and independent 
of time ¢,{ Then the solution of 4“ from Eq. (32) 
follows the classical procedure for the solution 6“, the 


+ Thus the boundary itself must also be independent of time 
so that the conclusion following would not be applicable to 


boundary layers with rapidly changing thickness. 
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inhomogeneous terms on the right-hand side introducing 
a particular solution obtainable by quadrature. 

Little need be said about the inviscid approximation 
(obtained by dropping the viscous term). As a par- 
ticular example, suppose UL’) = LL”, ie., 

Uf = {] + ec”) (Vy) 
hence the time-stretching transformation [Eq. (S)| 
may be used as a check. It is easily verified that the 
inviscid solution for Eq. (32), if the boundary conditions 
permit, is simply 


d, (| —ic(a/B)& (33) 


One has, therefore, to 0(e), 


alx t) 
¢ 


[1 + e(tac B)e™ |@e (34) 


WV 
On the other hand, the earlier method gives 
v =e" ™ ot 
with 


r= tf (1 + e&)d(t/t) = t+ 3 
Developing e” ‘“’ for small €, one obtains again Eq. (34). 
For the viscous solution, as long as 8/a@ is small (of 
order y’’), one may focus attention near the critical 
point y = y, in the usual manner. Eq. (32) is thus 


reduced to 


; ; B d*@,y ;, 
uo — oe ~ — 
a dy* 
vp did,“ d*@ 
~~ rd 
1a dy! a, dy? 
where ©) is the viscous solution for the profile U 
by the classical theory. For points near the quasi- 
steady neutral curve pertaining to U’’(y), 
bb” = @%(¢), F£ = (vy — »)al’/r) (36) 
The perturbation ®,“'’ in such cases may be solved from 


, 


(d?/de*) By  — AF — $1) be" S i[(al.'/p) x 


UO /UMY + OO" /U "je" =—(87) 
where o 1(B/v) (al '/v) 
UO’, U, U™’ are values of dU /dy, U™, dU /dy, 
respectively, at y = y, and the primes indicate differen- 
tiation with respect to y. The solution #4"’’, and 


hence #4“, are thus obtainable by quadrature. In 
fact, if one puts 
69" = —1U)(a/B)OO”’ + Vb,” (38) 


then A®,“!)”’ satisfies the following: 


(d? do7) AD, Lyre u(¢ _ (,) AD, <a 
it(U."'/U ©'@©” (39) 
Hence even for U“ = U—.e., for similar profiles 


described by Eq. (6)—there remains the solution Ad," 
which cannot be absorbed by the scale-stretching 


transformation. It may be recalled that the viscous 
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solutions are important for flows with a solid boundary 
which usually lies at a finite value of ¢ from the critical 
point. The influence of Ad,“ under these circum- 
stances is not likely to be trivial. 

The procedure for the solution of the disturbance 
may now be summarized as follows. First, the two 
inviscid and two viscous solutions of &” are the classi- 
cal ones for the prescribed U(y). Given a, Re, and 
the boundary conditions, the eigenvalue c and the 
eigenfunction ®, are then determined. These are to 
be put into Eq. (32) together with the given 8 and LU”, 
and the solution ®,“? is now uniquely determined by 
the boundary conditions. In practice one may again 
prefer to handle the calculation by splitting it into the 
inviscid and viscous solutions as discussed above. At 
any rate, the disturbance 


V = ¢ tals) Tg 0) + ef p,() | 


is thus completely determined except for the ampli- 
tude which is arbitrary in the linear theory. 
To assess the momentary stability at various times, 


the growth rate must be calculated. By definition, 


u = OW/dy = Realfe*"—” (06,/dy)} L 


PS... (40 
oW/Oy = Real} —tae’*” ™ ®,} f ) 


vA 


U 


where & = © + ee”, is the eigenfunction. 
Thus, with subscripts R and / denoting real and 
imaginary parts, respectively, 
9 —iact Be Dae 9 —tact 
u? = (e~'** @,') |p? cos? ax — (e7'™ ,") p 

e ~'*¢@,’),-sin Zax + (e7'* ,’),? sin? ax 

I e/jl 

Hence, since the integration is to be taken over one 


wavelength, 


ff u-dx dy = x ga “f \®,’ |2 dy 


where c; is the imaginary part of c. In a similar way 


one can evaluate 
SS vdx dy 
Finally, 
E = (c?#"/2) f [|@,’ |? + a?|&, |2]dy , 
S [Uo +4. eT, U™ |dy 
= (e7*/2)Eo(1 + 22,71) + O(e*) (41) 
where 
a f l |, wr le + a? lp, 2 |dy J S(U®)‘dy 
(Sf [®,( DY! 4b gh Ob \dy 
Real - le o’l2 4,2] - ; 
S [/@. 2 + a ©, |? |dy 


Sf U©UMdy) 
(U)‘dy J 


Ey, 


(42) 


Differentiating, one gets 
G, = (1/2E)(dE/dt) = ac; + BIE, (43) 
Thus 


AllV 


ie 
G; 5 Q as ac; 


MAY 1961 


For the special case U = U™ and in the inviscid 


limit, Eq. (24) holds. The function /; turns out to be 

FE, = k—] (45) 
with k = It is easy to verify that the cri- 
terion Eq. (44) is identical to the previous one, Eq. 
(13), when both are applicable. 


—ac;/B. 


Concluding Remarks 


(1) In extending the usual stability theory for infini- 
tesimal disturbances to time-dependent basic flows, one 
must modify the criterion of stability. It is proposed 
to introduce a concept of ‘‘momentary stability,’ which 
is said to prevail at the instant if the kinetic energy of 
the disturbances, as a fraction of the kinetic energy of 
the basic flow, tends to decrease. The time rate of 
change of the kinetic energy ratio is essentially what is 
termed the ‘‘growth rate’’ in this paper. The ‘‘growth 
rate’ is so defined that it coincides with the amplifica- 
tion rate in the case of steady flows. 

(2) It is more difficult to correlate the time history of 
the momentary stability of a time-dependent basic flow 
to the laminar-turbulent transition phenomenon. For 
instance, the significance of an interim period of momen- 
tary instability, followed by momentary stability, is 
rather obscure. Only if the flow turns out to be always 
momentarily unstable after a certain instant, may one 
infer that transition to turbulence must occur sooner 
or later after that instant. Likewise, a minimum 
critical Reynolds number can be defined if the basic 
flow is always momentarily stable, for disturbances of 
any wavelength, below the particular Reynolds number. 

(3) The Orr-Sommerfeld equation for disturbances 
in an incompressible, parallel but time-dependent basic 
flow is examined. In the inviscid limit and for basic 
flows of the type 


U(y, t) = T (t/t) U*(y) 


the disturbance equation is reduced to a form identical 
with that for the usual steady flow, after a stretching 
of the time scale. The conclusion is that momentary 
stability will prevail if 


ac; < (1/U)(0U/dt) 


where ac; is the amplification rate calculated for the 
instantaneous profile of the changing basic flow. From 
examining typical stability calculations, one is led to 
predict that when the time constant of the basic flow 
is much shorter than 10° X (the time constant of the 
disturbances of the instantaneous profile), accelerating 
flows are likely to be always momentarily stable and 
decelerating flows always momentarily unstable—for 
the type of disturbance considered. The shape of the 
velocity distribution loses much of its importance. It 
is believed that Fales’ and Coles’ experimental observa- 
tions are largely if not predominantly attributable to 
the effect of deceleration revealed in this paper. 


(Continued on page 417) 
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Influence of Transverse Shear on Nonlinear 
Vibrations of Sandwich Beams With 
Honeycomb Cores 


HU-NAN CHU* 
The Martin Company 


Summary 


In this paper are derived the nonlinear equations of motion for 


a honeycomb-core sandwich beam. The nonlinearity comes 
from the moderately large deflections 


The derivation utilizes the variational 


the linear stress-strain 


law is still assumed. 
principle and takes into account the effect of transverse shear 
deformation. Solution of nonlinear periods is determined for a 
beam with hinged, immovable ends. 

Three numerical examples are worked out for sandwiches with 
aluminum-alloy facings and aluminum-foil honeycomb core. 
The cores have transverse-shear moduli ranging from 22,000 to 
150,000 psi. The thickness-length ratios used are one tenth 
(1/10) and one twentieth (1/20). 

The conclusion is that in nonlinear vibration analyses of 
honeycomb sandwich constructions commonly used in flight 
structures, the influence of transverse shear deformation can, 


asarule, be neglected. 


(I) Introduction 


| Spensaieon VIBRATIONS OF HOMOGENEOUS BEAMS 
have been investigated by A. C. Eringen,' while 
H. Chu and G. Herrmann? have considered nonlinear 
vibrations of homogeneous flat plates, and E. Reissner*® 
has reported on nonlinear vibrations of homogeneous 
cylindrical shells. 

In these few examples from the available literature, 
“The nonlinearity of the problem,’’ to quote V. V. 
Novozhilov,‘ “arises from the geometry.” 

To be more specific, the rotations of the cross sec- 
tions are sufficiently large to be of consequence in such 


strain-displacement relations as 
€r = (Ou/Ox) + (1/2)(0w/dx)? 


where e, is the lengthwise strain; w and w are the dis- 
placements in the directions of length and thickness 
respectively; and x is the lengthwise coordinate. 
Terms such as (1/2)(Ow/Ox)? arise from large rotations 
which are absent in the linear theories. Note that the 
results of references 1, 2, and 3 apply only to wave- 
lengths large compared with the thicknesses—i.e., the 
effect of transverse shear deformation has been neg- 
lected. 


Received March 21, 1960. Revised and received June 20, 
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In sandwich constructions, owing to the weaker core, 
transverse shear deformation affects a much wider 
range of wavelengths® than in the corresponding homo- 
The wavelengths in question usually 
vibration 


geneous cases. 
include those for which finite-amplitude 
analyses are of interest. Therefore, it is necessary to 
determine the possible influence of the transverse shear 
deformation on the period of nonlinear vibrations. 
This is done here by analyzing the finite-amplitude 
vibrations of a honeycomb-sandwich beam, and _ in- 
cluding transverse shear deformation in the strain- 


displacement relations. First a linear distribution of 











‘ } 
: -, 
lan 
i Oo — 
—_ = 
-<~ fF 
Fic 1. Cross section of a sandwich beam 


the displacements across the thickness is assumed, and 
then the equations of motion are derived, with the use 
of the principle of minimum potential energy. A form 
of solution is then assumed and periods are found which 
are then compared with periods found by neglecting 
transverse shear deformation. 


(II) Equations of Motion 


(a) Strains and Displacements 


The three-dimensional strain expressions for moder- 
ately large rotations are given, for example, in Novo- 
zhilov’s book. When applied to the sandwich beams, 
the three-dimensional strain expressions reduce to the 
following: 
= Crzri . (1 2) (€ze1 - 


€rzi 


wy)? | 
f (1) 


where €,;; = OuUz;/OX, 2¢:4 = | 


(Ou,;/O2) + (Ou,;/Ox)> (2) 
2w,, = (Ouz,/02z) — (Ouz, ox) ) 
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The subscript 7 takes on the values 1, 2, and 3 denot- 
The letter 
e refers to the strain expressions in the linear theory, 


ing the various layers, as shown in Fig. 1. 


while w, is the rotation. 

Eqs. (1) are only approximately correct, in that 
€,,; and e€,,; are not considered. Accordingly, this dis- 
cussion must be limited to frequencies that are suffi- 
ciently low to render ¢,,; and ¢,,, unimportant. It 
should also be noted that the strains in Eqs. (1) are 
not necessarily small compared with the rotations. 

The displacements are then assumed to be as follows: 


U(x, 3, t) = u(x, t) + s(x, t) 


U(x, 2,t) = u — h(w’ + p) — sw’ (3) 
Up3(x, 2,1) = u + h(w' + vy) — sw’ . 


Uzi(x, 2, t) = w(x, t) 


where the prime notation (’) indicates differentiation 
with respect to x, the lengthwise coordinate of the beam. 
y may be called the shear angle of the core. Trans- 
verse shear deformation of the faces has been neglected, 
since the faces are assumed to be thin as compared 
with the core. 

The strain-displacement relations are then deter- 
mined by use of Eqs. (1), (2), and (3): 


€rn1 = U’ + (1/2)w’? + sy’ 
ep =u —h(w +y’)+ 
é73 = u' + hw" + yp’) + 


(b) Stress Equations of Motion 


To derive the equations of motion, we shall use what 
is most commonly known as the principle of minimum 
potential energy, in which variation is performed only 
on the displacements. To save space, however, the 
details of the derivation will be omitted. For these 
details, interested readers may consult a research re- 
port of the same title and by the same author as the 
present article, issued by The Martin Company, 
Denver Division, March 1960. The procedure is out- 
lined as follows. 

First, the total strain energy W is calculated with 
the use of Eqs. (4), and the total kinetic energy 7 and 
the total work done II’, by the external forces are found 
with the aid of Eqs. (3). Then the variational prin- 
ciple states that for all arbitrary time intervals (f2 — f)) 
the true state of stress corresponds to that in which 


The result is the following set of stress equations of 


motion: 


CO , 
(Ni2 + Ny) + FF, = 2(hy p; ++ hyp) ut 
Ox 
) ; —? h 
[Mu + n(Nes — Nez)] — On + : m, 
“ 37 , Ow 
2 pity = Iyhope | 2h — hy | | 
(ft) 
Oo .. 
~ (Map + Mas) — hy (Ni3 — Nz)) + 
Ox” 
ra) [ Vo +N ) Ow 10 | a de hy Om, 
ee iV 72 LV 73 Dw rl qd h i = 
(J + hake a ie oY 2 4,08 
- 2p2)W UiN2* p2 = 2h” 
_ ” + Ox? 
where 


—h, h 
Ma = [ Orr. 2 dz, M,;3 = f Orr3 2 dZ 
Jh h 


, - 785 yeh ee 
F, a Srls am mM, = fay h 7; = Sele —h (S) 
u, W, ¥, etc.) indicate differentiation with 
In Eqs. (7), the letter o denotes 


and dots 
respect to time. 
stress, while in Eqs. (8) /, and f, are the components of 
surface traction in the x- and the z-directions respec- 
tively. 

In deriving Eqs. (6), the fact has been utilized that 
resistance of the honeycomb cores perpendicular to the 
honeycomb flutes is small compared to that of the faces. 

The appropriate initial conditions are the specifica- 
tions for initial displacement and initial velocities. The 
appropriate end conditions are the specifications for 
one member of each of the following products: 

(Nie + Nz), [Mir + y( Nis — Ni) |W 
(Ay(Nes — Niaz) — (Maz + M,s) |(Ow/Ox) 


[Qu + Qe + Ors — (he/h)m,\w 


—h h 
where Q,o = [. fas, 03 = [s dz 


However, these conditions do not insure that the solu 
tion will be unique; they merely establish that energy 
will be conserved® within the approximations imposed 
earlier in this paper. 

For vanishing cores, /; > 0. Eqs. (6) reduce to 


(ON,/Ox) + F, = 2hopoti 


(0? /,/Ox?) + (0/0x)(N,Ow/Ox) + g + (Om,/Ox) 
Qheport’ — (2/3) poho*(0*zi/Ox") 


where N,-= Nw + Nz 
and M, = Mn + Ma + Mz 


These equations are the same as those shown in refer- 
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SANDWICH 


ence 6, except that here the notation 2h. replaces h 
and pe 1s p. 


(c) Stress-Displacement Relations 

[here are six three-dimensional stress-strain rela- 
tions for the orthotropic, transversely isotropic’ core. 
é,,, and ¢... can be eliminated algebraically among the 
six, while €,,; and €,. are zero when applied to a beam. 
Thus, this reduces the number from six to two for this 
case, namely: 


Orr = Fieéen = VAC yyI + (ky Ea)vaoin 


Gen = WnGexn (Ya) 
where the subscript z indicates quantities associated 
with the z-axis, the axis of elastic symmetry. Thus 
y, is the isotropic Poisson's ratio in the transverse plane 
(the x, y-plane), while v. is Poisson’s ratio associated 
with orthotropy. / and /£., and G; and Ga, are re- 
lated similarly. 

Likewise, for the two identical isotropic faces, we 


have 
Orr? Feoe + VveGyy2 + V2O222 
7 >) a 
Or = 2Gr€, 2 
Bee ; (9b) 
Orrs = Lv€zr3 + v2Gyy3 + ¥20223 
Or23 = 2Gorez23 





Next we integrate Eqs. (9a) and (9b) according to 
Eqs. (7), employing the strain-displacement relations 
of Eqs. (4). 
lowing are observed: 


In carrying out the integration, the fol- 


Integrals containing ¢,,. and o,;-3; are zero because the 
transverse shear deformation of the faces has been 
neglected. This is justified if the faces are both thin 
and heavy as compared with the core. 

Integrals containing o,,; and o-,; are neglected, which 
will necessarily limit our computations to frequencies 
somewhat below those of the thickness-stretch and 


width-stretch modes. This is justified because we are 


Pheu" + w'w") + F, = 
kiyh,? + 2hyheke) yp" =“ Nh.*k2w’"’ — 2h.Galp mt ww’) 4 
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concerned primarily with the interaction of transverse 
shear and large amplitudes on frequencies of the first 
branch which describes primarily the flexural modes 
The lowest thickness-stretch mode, for example, is 
associated with the extensional mode. Although, in 
nonlinear vibrations, the extensional and flexural mo 
tions are coupled through the nonlinear terms,® both 
the extensional mode and the thickness-stretch mode 
have frequencies much higher than the first-branch 
flexural mode.* Incidentally, this would furnish a 
justification for the neglect of the longitudinal inertia. 
However, this neglect will be justified by reasons that 
appear later in the paper. On the other hand, it is 
easily seen that the width-stretch modes have fre- 
quencies still higher than that of the thickness-stretch 
modes. 

The result is the following set of stress-displacement 
relations: 


Ny = Inky[u’ — hw" + W’) + (1 2)” 4 
(1/2)(h + hy)w" 
N33 = hoks{u’ + hy(w" + y’) + (1 2)w’ = 
(1/2 (h 7 hi ine’ F) 
Man = (2 yA hy yw’, On 2G i(y + w’) { 10) 
M,» —(F2/2)ho(h + hy) {u’ — hy(w" +’) + 
(1/2)w’?] — (12/3)(h? — hy*)w" 
M3 = (E2/2)holh + hy) [u’ + hy(w" + wv’) 4 
(1 2)w’? | — (Ek 3)(h? — hy 2” 


Ordinarily, a shear coefficient « should be associated 


with the expression for Q,. However, since the 
of this 


thin and heavy faces, and since « is nearly unity for this 


scope paper is limited to sandwiches with 


type of sandwich,® Eqs. (10) will remain unmodified. 


(d) Displacement Equations of Motion 


If we substitute Eqs. (10) into Eqs. (6) we obtain 


2( hyp, + hops) tt 


(h/h)m, = 


(7 spyhy® + 2pehy*he)W — hyhe* p.(dw Ox (11) 


, 


|—*- 3(h? = hy*) + hy hz? |E(0tw Ox?*) =_- 2hyhokey'"’ +(O Ox) [2hoEew’ (u’ i oW *) | ee 
2nGa (Wy + w’) +g + (Ae h)(Om, Ox) =2(Mipr + heop2)t + Iyhe*p.(OY/Ox —? 3 pehe*(O7z Ox?) 


If flexure of the core is to be neglected, then the term * ;/,/,*y” in Eqs. (11) must be put to zero. 
Next, if the nonlinear terms are put to zero in Eqs. (11), then the longitudinal motions will become uncoupled from 


the flexural motions, while the last two equations of (11) will describe two branches of flexural motions 


The first 


(lower) branch, which is the one presently of interest, has in-phase relationships between y and —w’ in free vibra- 


tions. 
analogous to a double pendulum. 


The second (higher) branch has out-of-phase relationships between y and —w’ 


The situation is somewhat 


(III) Solution for Free Vibrations 


In free vibrations, F,, m,, and g of Eqs. (11) are zero. Also, for primarily flexural vibrations, the longitudinal 
inertia term is small to a higher order than other terms and can therefore be neglected. 


We may therefore write, to a first approximation, the first equation of (11) es 


This is made clear by per- 


turbation procedures in reference 2. 
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QheEs (u” + w’'w"”) = 0 (12a) 
The remaining two equations of (11) are written 
(4 3fyh,° ot. 2hyhok2) yp" —_ Ihe? Eow'"' —_ 2hyGaly + w’) = (2 3PM," + 2 pohy*ho) a yh? p2( Ow Ox) (12b) 


[—? 3(h® — hy’) + Ihe? |E2(0'w/Ox*) —2hyheEw’’’ + (0/0x) [2heEow’(u’ + '/ow’?)| + 2nGaly + w’) = 
2(hypr + hops) + hyho*p(OW/Ox) — 2/3 p2h23(O%w/Ox?) (12c) 


If we impose the boundary condition that u = 0 at x = 0 and x = /, then Eq. (12a) can be integrated as follows: 


a I a72 9/2 13 
i= 2 ] a W ~ dx =a (15) 


where / is the length of the beam. If Eq. (12b) is differentiated with respect to x and then added to Eq. (12c), the 


following is obtained : 


[—?/3E2(h® — Iy*)](O4w/Ox') + 2/3 hyp’ + (0/Ox) [2heEew’(u’ + '/2w'*)] = 
2(hips + hoo2)w + (7/spihy? + Zpohy*he + hyhe?p2)(Ob /Ox) — (2/spehe* + Iyhe®p2)(0°w/Ox?) (14) 


If y is set equal to —w’ in Eq. (14), then the transverse shear is totally neglected, and we have 


| I 
—*/;([E\hy? + E.(h® — hy*)|(0*w/oxt) + | i [ w’? ax (0°?w/Ox?) = 
/J0 
2(hypy + heop2)t — 7/3[ ply? + po(h® — hy*)|(07%/Ox?) (15) 
which makes use of Eq. (13). If, on the other nand, the transverse shear of the core is considered, Eqs. (14) and 
(12c) become, by virtue of Eq. (13), 


;' hk, 
—?/3F(h® — hy*)(04w/dx4) + 2/3k hy Ww’? + | ] f wax (0°w/Ox?) = 
0 


2(hyp1 of hop2)w +. 3piMy° + 2 pohy*he a hyhz? p») (ov Ox) —(? 3p2he® + Iyhs* p2) (O07 Ox") (l6a ) 


hoF2 
[—?/3(h® — hy*®) + Iyho?| E2(0'w/Ox') — 2hPhokoy’’’ + | : 


2(hip, + hops) + Ihe? p2(OW/Ox) — ?/3poh2*(072/Ox?) (16b) 


{ w? ax (0?w/Ox?) + 2yGu(y + w’) = 
J 0 


The form of the nonlinear term suggests that we should do well to begin with an assumed solution for w. If we 


assume that 


w= AT(t) sin (2x/l) (17) 


InE, AT\? 
= f w’? dx = 2hek» (* ) 
l 0 2/ 
w, then y takes the 


If we also assume a solution for y and require that the time function in y be the same as 7‘(f) in 


The nonlinear term is then 


form 
vy = BT(t) cos (1x//) (18) 


from Eq. (12b). Furthermore, 7(t) would have to be simple harmonic in ¢. To explain this inconsistency, we ob- 
serve that for very small nonlinearities, Eqs. (12b) and (12c) describe two modes, one of which is ¥-motion, primarily, 
and the other primarily w-motion. We see that Eq. (12b) describes primarily the y-motion. Since Eq. (12b) is 
linear, the influence of large amplitudes is therefore essentially on the lower branch, i.e., primarily on the w-motions. 
Therefore, we turn to Eqs. (16a) and (16b) which describe primarily w-motions while taking into account the in- 
fluence of y-motions. 

If we substitute Eqs. (17) and (18) into Eqs. (16a) and (16b), we obtain 


a, are | 3(h2/h) | gr = 0 (19) 
b,? b,? B77* = 
Pi Pi k(E,/Es)(hy/h)? + 1 — (hy/h)? 


rae oe {| = ss — ia/h)r 26" es fer -0 
Pr Pi (2/m)?(ty/h)(1 — k) + re} (1/3)[1 — (hy/h)*] — (hy/h)(he/h) (a /A)R + (hte/ 2h) ]f 2 
(20) 
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where 
- (3* 3)6reCsi [1 wie (hy h)?*| 
>? = : 
Pi (m4? /3)6?} (hy /h)*® + 7,[1 — (hi/h)*] — 30. — R)r,(hi/h)(h2/h) (i /h) + (he 2h)\f + 40h) 1 + rr,) 
SCL — Rk) (hy /h) Csi? + reCsi?} (2/3) (1 — (hi/h)*] — (ai /h)(he/h)(2(y/ hk + (heh) |} 26? 1 
(Shy/h)(1 + rar,) + 1,(he/h)?[(hi/h)k + (2/3) (he/h) |r?é* 
and 
8 = A/Ge. 6 = 2h/l, 7, = he/hh, re = Ex/Gal (on 
Tp = P2/ Pi, Csi? = Gi/ pi, k = —Bl/Ar f ad 


Eq. (21) serves to determine k. We hope that for most practical purposes k as determined by (21) will not differ 
perceptibly from the k found by equating the coefficients of T* from Eqs. (19) and (20). If the difference is appre- 
ciable, then new solutions for w and y must be assumed to eliminate the over-determinacy ink. For example, if w 
and y can be expanded into series of product of functions of ¢ and x, then the new solutions may be had by taking 
more terms in the series, instead of the one-term approximations of Eqs. (17) and (18). The examples in the next 
section, however, will show that Eq. (21) is adequate. 

On the other hand, if transverse shear deformation is completely neglected, we then have, from Eq. (15), 


ae { 3(h2/h) Jer ” a 
70" 0" = = (23 
/ t (E, /Fe)(hy/h)*? + 1 — (hy/h)? 
where 
m (x? 3)6°rpCsi [1 —_ (hy h)*| 
Por - 9 /a\eo 2 ‘ail (24) 
(m?/3)6?) (hy /h)*® + 7,[1 — (hy/h)*}§ + 4(/h)C + rar,) 


For most honeycomb sandwich constructions, /,/ 2 is 
of the order of 10-4 to 10~-%. Eqs. (19) and (23) will 


therefore have the same form: 
(IV) Numerical Examples 


~ 
— 


T+ pT + p?F(g)T* = 0 (25 


The data in Table 1 are averages from references 9 


3(he/h)B? and 10 for aluminum honeycomb cores. 


where F(8) = . (approximately) (26) 
1 — (h/h)? 
TABLE 1 
From this we can immediately write down Comm 8 tevnaike) Core 2 (intermediate) Core 3 (strong 
il . ee . Ga 2.2 X 10‘ psi 6 X 10* psi 15 X& 10% psi 
* IK (hb * IK(b ‘al : ideal z . 
ri oa 2K (k) pr and ly ” 2K(k) po (27) E. 10 X 10* psi 25 X 10% psi 50 X 10% psi 
c ’ ‘Tl ~ . aad ~y 2 be Le Ti : 4 4 
T . 5, * | - bp) * E, 80 psi 200 psi _ 500 psi _ 
: pi . Pi pig 2.7 lb./cu. ft. 6 Ib./cu. ft 9 Ib. /cu. ft. 
where 
haa : : : The faces are made of aluminum alloy with -, = 10? 
T\* is the nonlinear period, transverse shear con- ' el ._ we. ous j 
: si and pg = 172.8 lb./cu. ft. The thickness ratio 
sidered ; oy 
oF , a r, = 2/23 is assumed for all three cases—that is, the 
T) is the linear period, transverse shear con- : i. . ; : 
pwn core accounts for 92 percent of the total thickness of the 
sidered ; sas os tae ; 
iia : : beam. The & values in Table 2 were obtained by use 
T,* is the nonlinear period, transverse shear not f Eq. (21) ; 
; Or Eq. (a1). 
considered : 
T, is the linear period, transverse shear not 
° Ty >) 
considered ; TABLE 2 
K(k) is the complete elliptic integral of the first Core 1 Core 2 Core 3 
kind ; k 
: aye (6 = 1/10) 0.9218(0.9204) 0.9698(0.9693) 0.9877 (0.9875) 
k = V F(@)/2[1 + F(8)] b 
and pi pi* ee Po po* = V1 ae F() (6 = 1/20) 0.9790(0.9789) 0.9922(0.9922) 0.9969(0.9969) 


linear frequency 

"nails frequency The values in parentheses are obtained by equating the 

coefficients of the nonlinear terms of Eqs. (19) and 

In other words, 7\*/7, = 7;*/7T>, and hence the con- (20). This shows that the transverse shear effect 

clusion: The effect of transverse shear deformation does increase with the amplitudes, since the values in 

on the nonlinear vibrations is negligible if it is negligible 
for the corresponding linear vibrations. 


the parentheses are generally smaller. The increases, 
however, are negligible over the ranges usually of 
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interest. Therefore the method of approximation and the manner of determining k are adequate for most practical 
purposes. Thus, the importance of transverse shear deformation is determined entirely by the corresponding 


linear vibrations. 


The influence of transverse shear on the linear periods can be seen from 


_s \: 
T, Pi T° 


The data in Table 3 were calculated from Eq. (28). 
If we take the case of Core 1 where 6 = 1/10 and the 
transverse shear is about 8 percent of the rotation of the 
cross section, we see that the transverse shear influence 
is less than 0.1 percent with regard to the frequency. 
However, it must be realized that metal honeycomb cores 


TABLE 3 


Core 1 Core 2 Core 3 
T/T; 

(3 = 14/10) 0.999286 0.999766 0.999917 
Yr a & 

(6 = 1/20) 0.999952 0.999985 0.999995 


are generally stronger than many other types of cores. 
For example, cores made of cellular cellulose acetate 
are generally less than one tenth as strong as the aver- 
age aluminum honeycomb cores. Therefore, the effect 
of transverse shear in cellular cellulose acetate-core 
sandwiches can become considerably important. 

In designing honeycomb sandwich constructions for 
flight structures, a rule-of-thumb assumption would be 
a core about 10) times as thick as the total faces and 
accounting for only one third of the total weight, ap- 
proximately. Of our examples, therefore, Core 3 
would seem most appropriate, with its density of 9 
Ib./ft.? Therefore, in nonlinear vibration analyses of 
honeycomb sandwich constructions used in flight struc- 


+ + 


~ 


hy hy (7 hy ) 
k)r, + 
“ih & ME 2h 


(2S) 


(") 4 (: "] 24 4 (‘) (1+ ' 
r = 6° in) 
3 L\h , he ii * 


tures, the influence of transverse shear deformation 
can, asarule, be neglected. 
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A Study of Cotangential, Elliptical Transfer 
Orbits in Space Flight’ 


WILLIAM LI-SHU WEN* 


Ryan Aeronautical Company 


Summary 


The purpose of this paper is to determine suitable cotangential, 
elliptical transfer orbits between a given pair of coplanar ter 
minal orbits in space flight. The condition of tangency of two 
noncoaxial, confocal conics is established first. This condition 
is then used to determine the possible cotangential transfer ellip- 
ses between two coplanar, noncoaxial, confocal conic orbits of 
practical importance. These terminal orbits may both be 
ellipses or one may be an ellipse and the other a hyperbola 
They may intersect with each other or make no contact. What- 
ever the case may be, it is found that there are many such transfer 
ellipses possible between a given pair of terminal orbits. The lo- 
cus of the centers of these transfer orbits can be described by an 
auxiliary conic completely defined by the initial conditions of the 
problem. By use of the coordinates of this auxiliary conic, 7 
and m, aS parameters, the variations of the transfer-orbit ele- 
ments (semimajor axis @ and linear eccentricity c) with its 
orientation angle / can be determined by either analytical or 
geometrical methods. Once the functions a(/) and c(/) are deter- 
mined it is easy to calculate, with reference to orbit orientation 
1, such functions of a and c as the variations of points of departure 
and arrival, as well as transfer energy and transfer time. With 
this information available, the transfer-orbit optimization (or the 
selection of a suitable orbit to meet a specific requirement) can be 


easily made. 


Symbols 
A = apoapsis 
a@ = semimajor axis 
b = semiminor axis 
c = linear eccentricity 
E = eccentric anomaly 
é€ = numerical eccentricity 
F = common focus of conics 
G = universal constant of gravitation 
/ = orientation angle (Fig. 3) 
M = central of mass 
M = mean anomaly 
m = mass of the moving body 
O = center of conic 
P = periapsis 
pb = parameter of conic orbit 
ry = radial distance from center of attraction 
r4 = apodistance 
rp = peridistance 
S = separation of centers of two Confocal conics 
7 = period of orbit 
t = time 
V = velocity in flight path 
y = angle between major axes (Fig. 1) 
A = discriminant 
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n = true anomaly 
= angle between major axis and line connecting orbit 


centers (Fig. 3) 


Subscripts 
1 pertaining to or associated with terminal orbit No. 1 
2 pertaining to or associated with terminal orbit No. 2 
0 pertaining to auxiliary conic 


(I) Introduction 


The problem of cotangential elliptical transfer be- 
tween two coplanar, noncoaxial, terminal ellipses was 
first investigated by D. F. Lawden.' From purely 
geometrical argument, Lawden established the relation- 
ship between the elements of the transfer orbit and 
those of the terminal orbits, and found that between a 
given pair of terminal ellipses there are an infinite 
number of possible cotangential transfer ellipses. To 
find the optimum transfer orbit (with respect to energy 
or time), he proposed a geometrical-graphical method be- 
cause, in general, closed-form solution is not feasible. In 
treating the problem of interception via tangential 
ellipses for confocal elliptic orbits, Gedeon® has derived 
some analytical expressions which can be used for 
digital computer computation of the polar coordinates 
of the intercept point and the characteristics of the 
transfer orbit when the coordinates of the release point 
are given. 

In this article, we treat the whole problem by an 
alternative approach. At the same time, the scope 
of discussion is expanded to cover more general cases, 
i.e., ellipse-ellipse-ellipse transfer and hypertola 
ellipse—ellipse transfer. 

By analytical reasoning, we first establish the condi- 
tions of tangency, intersection, and noncontact between 
two coplanar confocal conic orbits in terms of the orbit 
elements a (semimajor axis), c (linear eccentricity), 
and / (orientation angle). The conditions so formu- 
lated express more clearly the physical situations and 
the geometrical relationships of the problem. These 
conditions are then used to determine the elements (a, c, 
/) of the cotangential transfer ellipses between two termi- 
nal orbits. Both analytical and geometrical methods 
of solution are formulated. Other important items 
during transfer maneuver (such as points of transfer, 
transfer energy requirement, and transfer time) are 
then determined from a and ¢ as functions of transfer 
orbit orientation /. Thus, all the information needed 
for transfer orbit selection is available. 
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(II) The Geometrical Relation Between Two 
Coplanar, Confocal, Noncoaxial Orbits 


Consider two coplanar, confocal, noncoaxial orbits 
| and 2, orientated as shown in Fig. 1 where F is their 
common focus, 0; and O, are their centers, and y is the 
angle between the major axes, always measured counter- 
clockwise at periapsis end. If these two orbits have a 
point y in common, then it can be expressed by the 


polar equation of either orbit as 


Pi pr 
1 + ¢; cos m 1 + é2 COS ne 


pe 
(II-1) 
1 + @ cos (m — y) 


Eqs. (11-1) yield a quadratic equation in cos m 
(B? + C*) cos? » + 2AB cos n, + (A? — C*) = 0 
and its solution is 


—~AB = CW CR + C) — A’ 


— (B? + C2) 
where 
(II-2) 
I I e| €2 
A= —-— B= - cos 
Pi pz Pi 2 


C = (e?/po) sin ¥ 


The geometrical relation between two confocal orbits 
is then expressed by the following conditions: 


A = (B? + C’) — A? 


pi Pipe pe pi pe 
2 1 9 


rip a EO a ee 
py pe" pipe 


— ¢)€2 COS y) 


> (0 twoconices intersecting at two and only 
two points i 
“ ' .} (if) 
=() two conics tangential at only one point 
<0 two noncontacting conics 


From the relations 


9 


b? = a*(1 — e?) = a? — c? 
and p = a(l — e*) = B/a 
we have 


I 


= [c:> — 2cyc2 cos y + C2? — (a1 — az)? | 
b,*b.? 
Since b? is positive for ellipses and negative for hyper- 
bolas, the conditions (II-3) can be reduced to the 
following two cases (see Fig. 2): 
(i) If the two conics are of same kind, we have 


S? = ¢" + 2c\c2 cos y + C2? 


= a" = 201C2 cos + ce” = (dy — Ge)" 
< 


intersecting 
for tangential 
(II-4a) 


noncontacting 
(11) If the two conics are of different kinds, we have 


(a, -+ ay)" 


VIIA 


intersecting 
for < tangential (II-4b) 


noncontacting 


where a and ¢ are positive for ellipses and negative for 
hyperbolas, while S, being the separation of orbit cen- 
ters, is always positive. The conditions for the special 
cases of one or both orbits being circular (c = 0) or 


9 


1) can easily be derived from Eqs. (II-3) 


parabolic (e = 
and (II-4). 
From Eq. (II-3) the condition of tangency is B? + 
C? = A*. Then from Eq. (II-4) the true anomalies of 
the point of tangency are 
B (@2/p2) cos y — (e1/p1) 


cos 7 = — — 


A (1/pi) — (1/po) 
or cos nN. = COs (m1 —_ Y) (II-5) 
(€2/pe) — (e:/pi) cos 


(1 Pi) —_ (1 pz) 


The radial distance 7 to the point of tangency may be 
obtained by substituting the first equation in (II-5) into 
Eq. (II-1); by virtue of Eq. (II-3), we have 


9 (1/pi) — (1/p2) U aT-¢ 
r=2 -O) 
[(1 — e1?)/pi?] — [1 — e2)/po? 






A. ELLIPSE AND ELLIPSE 


B. HYPERBOLA AND ELLIPSE 


Fic. 1. Orientation of two coplanar, confocal, noncoaxial conics 
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C. HYPERBOLA - ELLIPSE 


Fic. 2. Geometrical representation of two confocal conics 
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P 
P. F = COMMON FOCUS 
A = APO-APSIS 
P = PERI-APSIS 


0 = CENTER POINT 
¢ = ORBIT ORIENTATION ANGLE 


Fic. 3. Geometrical representation of a cotangential transfer 
ellipse between two terminal ellipses. 
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Fic. 4 


Fic. 5. 





Linear eccentricity and orientation angle of transfer 
ellipse between two intersecting ¢ llipse Ss 





Characteristic conics of two intersecting terminal ellipses 
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— ——— AUXILIARY ELLIPSE 


Fic. 6. Auxiliary conics of noncontacting terminal ellipses. 


It is to be noted that the second expression in Eqs. 
(II-3) is the same as that derived by Lawden! through 
geometrical argument, whereas Gedeon,* using a dif- 
ferent analytical approach, has derived the intersecting 
condition which is equivalent to the third expression 
of Eqs. (II-3). 


(III) Cotangential Transfer Between Two 
Coplanar Orbits 


For a given pair of coplanar, confocal terminal orbits 
characterized by the parameters (a, ¢, /;) and (do, C2, /2) 
(see Fig. 3) that may both be elliptical or hyperbolic, 
or one elliptical and the other hyperbolic; let (a, c, /) 
be a transfer orbit sharing a focus with the terminal 
orbits and tangential to both. Applying condition 
(II-4), we have the cotangential conditions: 


pO +Vo? — 2qccos (1 —h) +e = +(Q4 —- a)| 


I 
I 


(ade — a)| 
(III-1) 


So $V 6.2 — 2o0c cos (lz — 1) +c 


We also have 


S= +Vo. — 2¢;c2 cos (le — hh) + c? = + (a, — Qe) 


(III-2) 


depending upon the nature of and the relation between 
two terminal orbits. Bearing in mind the trigonomet- 
ric relations in Fig. 3, we obtain from the cotangential 
conditions (III-1) the expressions 


St = So = + (a — a2) od Ss 


or Si + Ss = +(a1 — a)>S (III-3) 


It is noted that the second part of Eqs. (III-3) speci- 
fies the relation between terminal orbits according to 
(III-2), while the first part of Eqs. (III-3) says that 
the locus of the centers (0) of all cotangential transfer 
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orbits is a conic with terminal orbit centers (O; and O,) 
as loci. For briefness, we shall refer this locus as the 
auxiliary conic (ao, Co, lo) with 


2\ao| = |a, — asl, 2\co| = 2|ao|eo = S| 
(III-4a) 
h=1 + o, or Lit ey ck € (4 = 1, 2)4 
sin ¢ SiN ge sin |/, — |; 
C2 Cj S 
: ay(l1 — ey") | 
S == (III-4b) 


1 + &9 COS no 


ro = radial distance of O from O;,, 


ny is measured from /) 
The cotangential transfer orbit (a, c, /) is obtained as 
@=4,; 2 3, | 
c? = ¢c;?7 + 2c,S; cos (mo + gi) + S? (III-5) 
ccos (l; — 1) = c; + S; cos (mo + ¢i) 
The type of the transfer orbit is determined by the 
following criteria: 
a —c=a;? — cP? + 2a,S;{1 + e; cos (yo + ¢:)] 
> 0 elliptic 
= (0 parabola 


< 0 hyperbola (IIT-6) 


It is easy to apply the general criterion (III-1) to 
(III-6) to the specific cases as follows. 


(1) Ellipse-Ellipse-Ellipse Cotangential Transfer 


In this case let us assume a > c > 0, @ > C2 > 0, 
a, > a. Then from Eggs. (III-1), (III-3), and (II-4a), 
we see that: 

(A) When terminal orbits are intersecting, the auxil- 
iary conic is a hyperbola with two branches (see Fig. 4). 
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AUXILIARY ELLIPSE 
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» § 
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ri=~— XS 
1+e,cos(, + 
\ 
\ 
P. \ REFERENCE LINE 
S=+ VC 2+2IC,1C, cos (¢,—¢ )+02 
la,l+a 7 _S 
a= a: Co=aceo= 5 


Fic. 7. Characteristic conic of intersecting terminal hyperbola 
and ellipse. 
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For the branch focused at O;, the associated transfer 


orbits are specified by 


a=a+Si>a 


Cc" = C1" _ 2019) Cos (no = ¢1) + 5? 
c cos (, — 1) = a — S, cos (m — gif (I-92) 
a*e— c= ay," - he + 

2a,S) [1 + €; COS (No — ¢1) | 


Since a? — c? > O for all values of m, this whole 

branch represents the locus of the centers of a group 

of exterior cotangential transfer ellipses with a > a). 
For the branch focused at Oo, the associated transfer 


orbits are specified by 


a= a2 — So< ad 

c? = C2” — 22S cos (no + ge) + Se? 

ccos (ly — 1) = c — S. cos (mo + ¢e) (III-8) 
a? — c? = a? — &* — 


2a2S2 [1 — €2 COS (no ss ge) |, 
The type of the transfer orbit is determined by the 
following conditions: 
(a2o/2)(1 — 07) 
1 + e: cos (r + mo + ¢) 


cotangential ellipse 


VIIA 


for (cotangential parabola —_ (III-9) 


cotangential hyperbola 


The right-hand side of Eq. (III-9) represents an 


ellipse (y’, 7’) focused at O, with 


l=h+am, 7 =r+mt+ ee 


We shall call this ellipse the characteristic ellipse (or 
conic). By applying condition (II-b) we see that in 
the present case the auxiliary hyperbola and the char- 
acteristic ellipse are intersecting. Therefore only that 
part of the auxiliary hyperbola which is enclosed by the 
characteristic ellipse Eq. (III-9) is the locus of the center 
of a group of interior cotangential transfer ellipses 
with a < dy» (see Fig. 5). 

(B) When terminal orbits are noncontacting the 
auxiliary conic is an ellipse with O2 as focus (see Fig. 
6). The transfer orbits are given by 
@=aoe+S=—-a- 

Cc? = Co? + 22S cos (yo, + p») + S,? 
c cos (ly — 1) = co + S2 cos (no + ¢2) 


=a? — * + 
2a2S>[1 =— €- COS (no ¢») | 


(III-10) 


Since a? — c? > 0 for all values of 7, all the possible 
cotangential transfer orbits are elliptical. 


(2) Hyperbola-Ellipse-Hyperbola Cotangential Transfer 


In this case let a; < 4 < 0,@2> > 0. From Egs. 
(III-1), (III-3), and (II-4b) we see that: 


(A) When the terminal orbits intersect, the auxil- 


ELLIPTICAL 


TRANSFER ORBITS 415 





AUXILIARY HYPERBOLA 
§,—S,=la,l+a,<S 


S, =la,l+a,S,=a—a P 
6, =f +? 
Fic. 8. Auxiliary conic of noncontacting terminal hyperbola and 
ellipse 
TABLE 1 


Geometrical Determination of Cotangential, Elliptical Transfer 


Orbit 
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Points of orbital transfer 


Fic. 9. 
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iary ellipse and the characteristic ellipse are as shown 
in Fig. 7. The transfer ellipse elements are a = az — 
So; c? and / are given as before by the expressions in Eq. 
(III-10). 

(B) When the terminal orbits are noncontacting, 
the auxiliary hyperbola is as shown in Fig. 8. The 
transfer orbit is specified by a = a, + S»; ¢ and / are 
given as before by the expressions in Eq. (III-8). 


a? — c? = a? — Cc” + 


2a2S2 [1 + e2 cos (no + g2)| > O 


Since 


for all values of 7, all the possible contangential trans- 
fer orbits are elliptical. 


(3) Circular Cotangential Transfer Orbit 


It is easy to show that the necessary and sufficient 
condition for the existence of a cotangent transfer 
circle (c = (0) is as follows: 


(A) Ellipse—Circle—Ellipse Transfer: a, > a2 > 0 


Intersecting terminal orbits: 
@=% = fa OF @ = Np = fop 
Noncontacting terminal orbits: 
a= "p = Pea 


(B) Hyperbola—Circle—Ellipse Transfer: a, < 0, 
de - 0 
Intersecting terminal orbits: 


a = hip, > Tor 


Noncontacting terminal orbits: 


c= Tip = Tea 


In all cases the center of the transfer circle is at the 


common focus F. 


(IV) Computation of Transfer Ellipse Elements 


From the above analysis the transfer ellipse elements 
a, c, and / are expressible as known functions of the true 
anomaly (no) of the auxiliary conic. Then the vari- 
ations of a and ¢ with / can be calculated by means of 
digital computer according to the formulas derived in 
the previous two sections. In this way, the calculation 
procedure becomes purely analytical. 

If a high degree of accuracy is not required, as in a 
preliminary study, the geometrical method proposed 
in Table 1 may be used conveniently to obtain quick 
results. 

In Table 1 we specify a given conic by a set of four 
parameters (a, c, /, F), where the parameters a, c, / 
are as previously defined. F signifies the location of 
the focus. Referring to Fig. 3, we are given, for the 
two terminal orbits, (a1, 4, 4, F), center point O,, and 
(dz, C2, le, F), center point O2; (a, c, /, F) and center 
point O for the transfer orbit are to be determined. The 
procedure for geometrical solution is as follows: 
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(1) Draw the major axes of the terminal orbits 
through the common focus F according to their orienta- 
tions, and measure S, ¢;, and ¢». 

(2) Construct the auxiliary conic (70, 79) and the char- 
acteristic conic (7’, »’) according to Table 1. 

(3) Measure c, /, ro for different applicable values of 
No. 

(4) Calculate a according to the formulas in the table. 


(V) Points of Transfer, Transfer Energy, and 
Transfer Time 


The points of transfer (points of cotangency), the 
energy required for orbit change at transfer points, 
and the transfer time can all be expressed in terms of 
the transfer orbit elements a and c. Therefore, they 
are also functions of transfer orbit orientation /, and 
can be determined easily once a(/) and c(/) are obtained. 
The desired expressions for these important items will 
be derived in this section. 


(1) Points of Transfer 


Referring to Fig. 9, let 


(r1,.m) = point of tangency between (a), , /;) and (a, 
c, /), where n; is referred to terminal orbit 
(a1, C1, Li); 

(2,42) = point of tangency between (de, c, /2) and (a, 


c, /), where yp is referred to terminal orbit 
(Qe, Co, I). 


Applying Eqs. (II-5) and (II-6), we have, letting 
i= 1,2, 


ab,;? — a,b* c;b? — cb; cos (1; — 1) 


r, = 2 : 


6° — b 


cos — = . re 
ab;? — a,b? 


(V-1) 
Thus when a, c, and hence } are known, the coordi- 
nates of points of transfer can be calculated by above 
equations. Conversely, if (71,m) and (72,n2) are given, 
the required transfer orbit can be determined from the 


following expressions: 


roky — rks ' " Fe “= Fy 
a= = gud c* = a? — 
2(ky —_ ko) ky — ky 
nN. cos |} — 1 cos I» 
tan/ = — - - 
nm, sin 1; — n, sin /» (V-2) 
where k; = (2a; — 7;)/b? 


a; a Ci 
and n, = b2 Be cos ni + b2 
b b? » ;? 


(2) Transfer Energy Requirement 


From the standard expression of orbital velocity, 
the orbital velocity changes (energy changes) at r; are 


2 1\1r 2 1\I12 
VGM ( - ) - ( = ) (V-3) 
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The total energy required to complete the transfer is 
AV = AV, + ADV, (V-4) 
(3) Transfer Time 


The flight time required in going from 7 to 72 along 
the transfer ellipse is obtained by applying Kepler’s 
equation and Kepler's third law, which yield 
Mt = Va®/GM [E, — E; — e(sin E; — sin E,)| (V-5) 


with cos £; = (a — 7;)/c 


La minar Stability (Continued from page 404) 
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(4) At finite Reynolds numbers, the quasi-steady 
solution, based upon the instantaneous profile, may be 
justified only if the time constant of the basic flow is 
much greater than Re’ X (the time constant of the 
disturbances of the instantaneous profile), Re being 
the instantaneous Reynolds number. Furthermore, 
it is still necessary to compute the “‘growth rate’ to 
arrive at a statement about the stability of the basic 
flows, even if the quasi-steady solution itself is accept- 
able. The quasi-steady ‘“‘neutral curve’ provides little 
information regarding the stability of time-dependent 
basic flows. 

(5) For slowly varying but otherwise arbitrary ve- 
locity profiles, the procedure of analysis (with certain 
restrictions) is outlined briefly. For points along, say, 
the quasi-steady neutral curve, both the inviscid and 
viscous solutions for the time-dependent basic flow may 
be thus constructed. On the basis of these solutions it 
will be possible to establish contours of constant growth 
rate similar to the contours for constant amplification 


or decay in the steady case. 
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of the fuel flow and the wave is stable at a much lower tempera- 
ture than that required to ignite it, was duplicated. Under these 
conditions, a flame near the tip of the hydrogen injector and a 
The 
hydrogen injector was modified to inject fuel at the throat of the 
tunnel. Under these conditions, no flame at the injector and no 
evidence of reaction upstream of the shock were seen. Also, 
although there was a temperature rise in the test section upon 
the addition of hydrogen from the injector, there was no move- 
ment of shock. No evidence of reaction was seen in the test 
section when the inlet temperature was reduced below the ignition 


sodium emission upstream of the normal shock were seen. 


point.”’ 
It seems clear that (1) the varying shock position and hystere- 
sis phenomena do not exist with constant upstream conditions, 


and thus (2) Kovitz’ analysis is not relevant. 
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Proceedings of Fourth Con- 


Transit Angle for Low-Thrust Spacecraft 


Frank M., Perkins 


Staff Systems Analyst, Engineering, Convair, 
A Division of General Dynamics Corporation, San Diego, Calif. 


September 30, 1960 


eben COMPLETELY GENERALIZED SOLUTION for trajectory pre- 
diction of low-thrust spacecraft employing a tangentially 
directed, constant-thrust acceleration and operating in a single 
gravity field is presented in reference 1. The relations between 
altitude, velocity, flight-path inclination, and time are presented 
in parametric form as single-valued solitary curves applicable 
for flight all the way from an initial circular orbit on out to 
infinity. These curves apply to all tangential, constant-accelera- 
tion, low-thrust flight paths starting in circular orbit at any 
altitude about any planet. 

Reference 1 omitted the solution for the angle subtended 
at the center of the gravity source by the travel of the spacecraft 
(range or transit angle). This solution is presented in Fig. 1 
as a number of revolutions, starting at any point on the trajectory 
and ending as the spacecraft approaches infinity, versus the alti- 
tude parameter XY. To the revolutions traversed 
between any two points on the trajectory, it is merely necessary 


determine 


to take the difference between the two points on Fig. 1. The 
altitude parameter X is defined in reference 1 as follows: 
X = a!*%(r/ro) = (f/p)?r (1) 

where r+ distance from center of planet 

f = thrust acceleration = force/mass 

a = gravity constant (distance*/time?) = rg 

g = gravity acceleration (at any altitude) 

a = dimensionless thrust acceleration = f/(u//07) 


and where subscript 0 denotes the initial circular orbit. 

The equation for the circular asymptote of Fig. 1 is obtained 
by putting Eq. (14) of reference 1 into local differential form, 
substituting the circular velocity relationship from Eq. (20) into 
this, and integrating to infinity. Of course, circular flight does 
not apply all the way to infinity, but for the purpose of obtaining 
the circular asymptote, it may be assumed to apply in the 


mathematical sense. 
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The slope of: the infinity asymptote of Fig. 1 may be derived 
as follows. By definition, the rate of change of the transit or 
range angle ¢ and radius distance r with time are, 

¢=( V cos @)/r 
* = Vsin#@ 
where @ is the flight-path angle measured above local horizontal. 
(4) 


Ata sin 6 — 1.0 


Dividing Eq. (2) by Eq. (3) and bearing in mind that Eq. (4) yields 


dg/dr = cos 6/r (5) 
Multiplying both sides by 7r/(¢a — ¢) and then putting 7 in 
terms of the altitude parameter X results in the following: 
din (¢a — ¢) —cos 6 —> 0 ; 
s = = = (6) 
din X Ya — ¢ — (0) 


Since both the numerator and denominator of the right-hand 
side of Eq. (6) approach Zero, the value of this fraction is equal 
to the ratio of the time differentials of the numerator and de- 


nominator at infinity. Performing the differentiation and 
substituting Eq. (4) into the numerator yields 
din (¢a — ¢) 6 a 
3 = (7) 
din X —¢ 


As the vehicle approaches infinity, the gravity which had 
caused the flight path to turn becomes negligible compared to 
the thrust acceleration, and the flight path approaches a straight 


line. Under these conditions, the rates of change of both 
flight-path angle @ and transit angle ¢ approach V cos @/r. 
Therefore, 
6=¢ (8) 
and the logarithmic slope of Eq. (7) becomes minus one: 
din(¢a — ¢) 
Ya : ¢$ a —e (9) 
din X 


It is interesting to note that although the flight-path angle 
changes from 39.2° at the parabolic escape point to approach 
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90° at infinity (reference 1), the vehicle completes only about 
0.17 revolutions during this phase of the trajectory. 

The same type of parametric generalization used for tangential 
thrust in reference 1 is applicable to the case of horizontal thrust 
as well. The parametric curves are omitted here to conserve 
space. There is very little difference in the trajectories up to 
the parabolic escape point, although the difference then becomes 
marked because of the high flight-path angle. There may be 
some navigational simplicity in maintaining horizontal thrust 
when an horizon scanner is used to determine local vertical. 
When the flight-path angle becomes appreciable, however, the 
energy gained by the vehicle per unit of propellant expended is 
much higher by directing the thrust along the local velocity 


vector. 
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Generalized Spectral Representation in 
Aeroelasticity: Part II 


Arthur K. Cross 
Norair Division, Northrop Corporation, Los Angeles, Calif. 
September 28, 1960 


IT A PREVIOUS NOTE,! a formalism was developed by which 
measured modal information can be converted into stiffness 
information. These dynamically measured stiffness data show 
in detail any defects present in theoretical stiffness data or in 
stiffness data determined by static testing procedures. 
for highly idealized wings with elastic axes were then outlined. 
It is proposed in the present note to apply the method to a general 
unrestrained, statically indeterminate structure with the aid 
of a greatly simplified formulation of Eqs. (8) to (12) of reference 
ie 

In terms of the notation of reference 1, the equations of motion 


Examples 


for unrestrained vibration are 
H =H — Long = w*CMH (1) 
—Lo'MH = Lig = —L2L;H (2) 
Operating on Eq. (1) with —Zo'M yields 
—Lo'MH + LotMLoqgo = —w*Lot MCMH (3) 
After elimination of —Zo'MH between Eqs. (2) and (3), there 
results 
—w*Lot MCMH = (Li + Lot MLo)qo (4) 
Consequently qo can ke eliminated between Eq. (1) and Eq. (4) 
so that 
AH = w{I — Lo(Li + LotM Lo) Lot M| CMH = w*KCMH (5) 
Alternatively, K may be calculated by eliminating go between 
Eqs. (1) and (2) with the result that 


K = (I + Loli Lot M)=! (6) 
Then the spectral representation for C is 
3m 
C= K7 >) o:°H; Hit/ HitMH; (7) 
i=1 


The above formulation is somewhat awkward, since it requires 
separate handling of the mode shapes at the reference station 
(the go matrix) and the mode shapes at the other stations (the 
H matrix). The following reformulation of Eqs. (1) to (6) 
is based on a straightforward coordinate transformation. 

Consider a general unrestrained structure having 6m elastic 
degrees of freedom and 6 rigid degrees for which a [6(m + 1) X 1] 
mode-shape matrix is denoted by Y where 


Yt = [H*, qo! (8) 
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Then H and YF are related as follows: 


H = UY (J, —Lo] 9 


where U 
[where J denotes a (6m X 6m) unit matrix and Lo is a (6m X 6 
generalization of the rigid-body mode-shape matrix Lo defined 
in reference 1]. 

The representations of the stiffness and mass matrices relative 
to the new coordinates Y are denoted by the symbols x and x, 
respectively, and are determined by the condition that the 
strain energy V and the kinetic energy 7 be invariant: 


V = (1/2)H'CH = (1/2)V*(U'C1U] Y (1/2) Y¥*xY 10 


T = (1 2)‘ (1/2) ¥uY 11) 


Thus the [6(m + 1) X 6(m + 1)] matrices x and yu are defined 


as follows: 
UCU - ° 12 
x , a O L, - 


The equations of motion in the new coordinates are found by 
varying the energy E, which is expressed in the new coordinates 
as follows: 


E = 7+ V = (1/2)[Y*xY — w? YY] (13) 
Then 
6E = O = 6Y"xY — wy V) (14) 


Since the variations 6Y are independent, 
xY¥ = wuY (15) 


For a system acted on by motion dependent or other forces 
described by the [6(m + 1) X 1] matrix F(t), the variational 
condition (14) is replaced by a condition involving the generalized 
work W = Y'F(t): 


le Pte = 
sf [L + W]dt = 0 = -f, 6V*(uY + xY — F(t)|dt (16) 
1 1 


The independence of the variations 6Y then implies that 
uVY+xY = R(t) (17) 


The above representation is particularly advantageous when 
it is desired to reduce the size of the matrices involved by con- 
verting to a modal description. Then, if the coordinate trans- 


formation 
Y(t) = SC(2) (18) 


is carried out, where S represents a [6(m + 1) X /p] matrix of 
expansion functions and C(t) is a (p X 1) matrix of generalized 
modal coordinates, the following variational condition results 
from the substitution of Eq. (18) into Eq. (16): 


te le se 
5 f [L + Wldt = — f 6CXt)[M C(t) + K,C(t) — StF(t)]dt 


tai 
(19) 
Hence the equations of motion for the forced system have the 
following form in the new coordinates: 


M,C(t) + K,C(t) = F,(t) (20) 
where 


M, S'S, K, S'xS, and F,(t) S*tF(t) (21) 


are the generalized mass, stiffness, and force matrices, respec- 
tively. 

The application of the coordinate transformation (9) to the 
construction of the spectral representation of the flexibility 
matrix of an unrestrained redundant structure is outlined below. 
It is assumed that m external loads described by an (n X 1) 
matrix E are applied to the structure under consideration. 
These loads produce m statically determinate loads and , re- 
dundant loads which are described, respectively, by an (m X 1) 
matrix D and an (ry X 1) matrix R. In practice, the equality of 
the number of external loads and statically determinate loads 
must be achieved by a judicious choice of lumped-mass locations 
so that using the same number of lumped masses (i.e., the 
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number of ‘‘external”’ or inertia loads) and statically determinate 
loads does not result in singular matrices T and V where the latter 
are matrices of the coefficients of the relations between the 


ipplied, statically determinate, and redundant loads. That is, 


D=TR+ VE (22 
The restriction does not apply to T for cases where ry # n, for 
In terms of the local flexibility-property 
determinate loads matrix D as 


then T is rectangular. ) 
matrix B and the statically 
given by Eq. (22), the strain energy may be expressed as 
V = (1/2)D‘BD (23) 
The variational principle (16) implies in the present case that 
6(V — W) = 6[(1/2)D'BD — EtY] = 0 (24) 
Substitution of Eq. (22) into Eq. (24) and explicit evaluation of 
the variation yields 
5R'[(T*BV)E + (T‘BT)R] + 
6E'|(V‘BT)R + (V‘BV)E — Y] = 0 (25) 
Since the variations 6R and 6E are independent, R may be 
eliminated between the two resulting equations with the well 





known result that 
Y = [V‘BV — (V‘BT)(T‘BT)-(T‘BV)|E = CE (26) 


According to Eq. (15) above and Eq. (7b) of reference 1, x can 
be determined from experimental modes Y; by the following 
expansion: 

6m 

x = >, wf@pY ¥en/¥in¥e (6m = n) (27) 

i=1 
The basic flexibility-data matrix B is determined from the 
empirically determined x~! matrix by noting that Eqs. (9) and 
(12) imply that the upper left (6m X 6m) submatrix of x is 
C-!,. Inspection of Eq. (26) shows that B can be determined 
by iteration—.e., if By is taken to be the theoretically determined 
basic flexibility-data matrix, 


Bata = (V*)O[C + (V‘B,T)(T'B,T)-(TB,V)|V- = (28) 


Examination of the elements of the empirical B matrix thus 
determined reveals any defects present in either the idealization 
of the actual structure or the static measurement of the geo- 
metric data or elastic constants. If a flutter model is con- 
structed with internal structure corresponding to that of the 
actual airplane, comparison of the dynamically measured B 
matrices of each can suggest structural changes in the model 
to improve the simulation of the airplane-stiffness properties. 


REFERENCES 
Generalized Spectral Representation in Aeroelasticity, 
26, No. 11, 


1Cross, A. K., 
Journal of the Aero/Space Sciences, Readers’ Forum, Vol 


November 1959 


pp. 766-767, 


A Drag Hypothesis for Jet-Flapped Wings 
G. K. Korbacher 


Assistant Professor of Aeronautical Engineering, 
Institute of Aerophysics, University of Toronto, 
Toronto, Ontario, Canada 

June 28, 1960 


Dopp THE HEADING ‘“‘Jet-Flap Thrust Recovery,”’ Foley and 
Reid! recently announced the “significant finding that sub- 
stantially complete thrust recovery has been obtained in two- 
dimensional-flow tests on a jet-flap wing.’’ This finding is based 
on experimental results (Figs. 3 and 4 of reference 1). The con- 
stant slope of the lines in Fig. 3 of reference 1, which is the same 
as the slope 0.94 in Fig. 4 of reference 1, represents a thrust- 
recovery factor of 0.94. This finding is substantiated by test 
results of NGTE* * and ONERA‘ as shown in reference 5. 

A value of 0.0243 for the slope dcp/dcz? is used in reference 1 in 
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the transformation of the experimental values of Fig. 3 to those of 
Fig. 4. The choice of this value does not affect the slope 0.94 
of the “experiment” line in Fig. 4 (reference 1), but it changes 
its intercept with the vertical axis. As this point of intercept 
can be obtained from the test results of Fig. 2 (reference 1) at 
Cu = 0- the value by which c,? 
has to be multiplied is indirectly fixed 
reference 1, the 


at least as an average value- 


According to value 0.0243 was estab. 
lished by measuring wing forces at 
with the slot sealed (c, = 0) and flap undefiected 
prefer, however, the use of the dcpr/dc 7? of the jet-flapped 
wing (cy, + 0), where the subscript T designates total values 


The use of (dcp/dcxz? ey =0 instead implies the ‘“‘drag hypothesis” 


various angles of attack 
One would 


(formerly suggested in reference 4) as 


(dcp dcx*)c,=0 = de pr dc.r? = ().0243 (1) 


No doubt, this relationship would be of great significance if it 
could be proved to be true. It seems, however, that the test 
results of reference 1 leave some doubt about the validity of this 
hypothesis. Additional information about (dcp/dcz? ey =0 is re- 
quired from researchers who work with two-dimensional jet-flap 
wings, in order to prove or disprove Eq. (1) for jet-flap wings in 
general. 

If one plots the crucial graph of the total drag vs. the total lift- 
squared from the test results of reference 1, Fig. 1 is obtained 
It shows that the drag hypothesis can be considered only as an 
approximation. 

Next, if all two-dimensional jet-flap test results?~* are con- 
sulted in order to establish whether dcp7/dcyr? is a constant or 
may be approximated by a constant over a limited c,-range 
(Cu S 1), one still does not obtain a unique answer.® 

Finally, if established empirical relationships and theoretical 
predictions are used, the following relationship is derived® 
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dcpr/dcyir? = 0.63/K?2 (2) 


Here K? is a specific constant for each particular jet-flap wing 
and can be obtained from* ® 


cir = Kc ‘7? sin @ (3) 
where @ is the jet deflection angle. The constant 0.63 follows 
from? 

Cpr = (co), =0 + a(O) cy (4) 


which after differentiation with respect to c, produces a(@). Now 
if from all available two-dimensional jet-flap tests results a(@) is 
plotted vs. sin? 6, Fig. 2 is obtained® and the result is 


d a(@)/d sin?@ = const. = 0.63 (5) 


Eq. (2) suggests that dcep7/dcz7? is a constant as long as K isa 
constant or can be approximated by a constant over at least a 
limited range of c,-values. It was found® that for c, < 1, K 
can be considered a constant, but that for c, > 3, this no longer 
seems to be the case. 

Spence’ has derived an expression from which it follows that 
K is a constant only for small values of c,; for higher c, = values, 
K becomes a function of c, and increases with Cy. Fortunately, 
the average c,-value for jet-engine-powered aircraft at the moment 
of take-off is as low as 0.5. A preliminary study’ of how their 
take-off characteristics can be improved with a jet-flapped wing 
indicated jet coefficients at take-off still below 2.5. This sug- 
gests that the practical range of jet coefficients coincides well 
with that for which the above relationships are applicable. 

Finally, the predictions of Eq. (2) may be compared with 
those of Maskell’s and Spence’s theory® for three-dimensional 
jet-flapped wings. Since Eq. (2) was derived from quasi-two- 
dimensional test results, this comparison can only be qualitative. 
From theory? follows 


dcpr/dcir? = 1/(rAR + 2Cyu) (6) 
This equation, for rAR > 2c, reduces to 
dcpr/dcyr? ~ 1/7AR = (dep ‘dex*)c, =0 (7) 


Further, 
Thus, 


This brings out the limitation of the drag hypothesis. 
Eq. (6) predicts a decrease in slope with increasing cy. 
the predictions of Eqs. (2) and (6) agree in principle. 


AEROSPACE 


SCIENCES—MAY 1961 


REFERENCES 

1 Foley, W. M., and Reid, E. G., Jet Flap Thrust Recovery, Journal of the 
Aero/Space Sciences, Vol. 26, No. 6, pp. 385-387, June, 1959 

* Dimmock, N. A., An Experimental Introduction to the Jet Flap, NGTE, 
Rept. No. R175, July, 1955. 

3 Dimmock, N. A., Some Further Jet Flap Experiments, NGTE, Memo 
No. M255, 1956 

4 Malavard, L 
and Experimental Investigations of Circulation Control, ONERA TN No. 37 


, Poisson-Quinton, Ph., and Jousserandot, P., Theoretical 


June, 1956 

5 Korbacher, G. K., and Sridhar, K., A Note on the Total Drag of Jet 
Flapped Wings, UTIA Rept. No. 64, May, 1960 

® Dimmock, N. A., Some Early Jet Flap Experiments, Aero. Quart., Vol 
8, Part 4, November, 1957. 

7 Spence, D. A., A Treatment of Jet Flap by Thin Aerofoil Theory, RAE 
Rept. Aero—2568, November, 1955; also, The Lift Coe ficient of a Thin Jet 
Flapped Wing, Royal Soc. Proc. Ser. (A) 238, 46-68, 1956. 

8 Korbacher, G. K., The Jet Flap and STOL, Decennial Symposium at 
UTIA, October, 1959. 

9 Maskell, E. C., and Spence, D. A., A Theory of the Jet Flap in Three 
Dimensions, RAE Rept. Aero—2612, September, 1958; also Royal Soc. Proc 
Ser. (A) 251, 407-425, 1959 


Author’s Reply 


Elliott G. Reid 

Stanford University, Department of Aeronautical Engineering, 
Stanford, Calif. 

June 28, 1960 


f bw SUBJECT OF THIS COMMENT is the final sentence of Pro- 
fessor Korbacher’s first paragraph; the statement therein 
that Foley’s finding ‘‘is substantiated’’ by earlier British and 
French work is incorrect. 

The fact is that the results of the cited experiments are so 
conspicuously inconclusive in their failure to demonstrate any- 
thing approaching complete thrust recovery that they actually 
provided the principal impetus for the Standford investigation 
which was directed specifically toward ‘‘verification or disproval 
of the jet-flap thrust recovery theorem.” 

Foley’s work resulted in the demonstration of thrust recovery 
in excess of 90 percent at jet deflection angles as large as 59.1 deg. 
This is hardly ‘‘substantiated”’ by the recovery of only 64 per- 
cent of the jet reaction at a slightly smaller deflection (58.1 deg.) 
in the NGTE work—the results of which are summarized in 
Fig. 4 of Korbacher’s reference 6. 
point (p. 335, reference 6) is significant ; he says, ‘‘The existence of 


Dimmock’s comment on this 


a pressure thrust is undeniable, but its magnitude relative to 
the full jet reaction was initially rather disappointing.” 

To the best of the writer’s knowledge, Foley’s finding was 
unique when published and has been substantiated only by 
the results of subsequent Stanford work—of which no account 


has as yet been published. 


On a Generalized Solution to a Torsion 
Problem in Linear Elasticity 


Howard E. Brandt 

Physics Student, Massachusetts Institute of Technology, 

Cambridge, Mass. 

August 10, 1960 

eee StaniSi¢é!~* with his students at Purdue University 
published a number of excellent papers concerning the 

torsion of prismatical bars with various cross sections. In 

each of these papers the application of the Galerkin method 

has been dexterously illustrated. The advantages of this 

method have been shown by Krylov and Kantorovich,‘ Stanisié,5 

Sokolnikoff,* and innumerable other scientists in the field of 

StaniSi¢é has developed the problem 


mathematical physics. 
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for the torsion of a prismatical bar with hexagonal cross section. 
However, this problem can be‘generalized, and the stress and 
displacement fields can be obtained for any N-sided polygonal 
cross section, using a transform technique illustrated in this 
paper. This generalization is a result of the topological charac 
ter of the stress function in the Prandtl sense and is associated 
with the properties of the Laplacian operator. Moreover, this 
problem has wide application in the design of transport jet 
engines. Hence, in that which follows, a technique will be 
developed in order to find the stress, displacement, and torsional 
rigidity for any prismatical bar of N-sided regular polygonal 
cross section StaniSié’s solution then becomes a special case 


of this study 


MATHEMATICAL FORMULATION OF THE PROBLEM AND ITS 
SOLUTION 


The solution to the problem will be found by use of the 
Prandtl stress function, applying the Galerkin method. This 
method for determining stress and displacement fields has been 
uutlined by Stanisi¢ and Shirely.2. Therefore, only the mathe- 
matical analysis peculiar to the immediate problem will be 
presented in this paper 

The stress function is assumed in the Prandtl sense, 


> ews ¥), (m 1, 2,.3, a (1) 


W(x, y) 


and the Galerkin form of the stated problem becomes 


ee . 
ff V2[W(x, vy) + 2] n(x, v)dS 0, (7 = 1,32. 3, ~s) 
sS “ . 


where 4,(x, y) are coordinate functions which must be found, 
such that the boundary condition y(T) 0 is satisfied on the 


boundary IT of the N-sided polygon; A,, are constants satisfying 
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polygon into N-congruent isosceles triangular surfaces, S 
So, , Sv-r Sy. The base vertices of S,,, m i ; 
N — 1, N, are the points, Pa-i(X%m-1, Yn-1) and Pn(Xm, Yu) 
A function yx, y) defined on the polygonal surface must be 
found such that the boundary condition is satisfied along the 


lines, L», passing through the points, P,, and P It is 


readily shown that 


we x cos (r/N)(2m — 1) + y sin (r/N)(2m l R 0 


Hence ¥(x, y) will be of the following form 

\ - - 
y(x, v) A, JJ | x cos — (2m — 1) + ysin 2n — | R 
sik oat N \ 


where A; is a constant of first approximation to be determined 
later and II denotes the product of expressions with indicated 


index. Hence Eq. (2) becomes 


Jf fi, iene, » 


2]m(x, vdS ) (9) 


where 
\ . ” 
m(x, V) x cos — (2n 1) + ysin 2n l R 
j ale A A 
(6) 
Solving for A, yields 
\ » 2 \ 
 ® J f. II ¢,(x, y)dS 
m=1 "n=l = 
(7) 


m= 


where the integrals are evaluated over the respective triangles, 


























Eq. (2); and S is the domain of integration. es 

We first find a function, ¥(x, y), which satisfies the boundary : : ; 
condition YT) = 0. For this purpose refer to Fig. 1. The Pnl(X, ¥) x cos (x/N)(2n — 1) + ysin(x/N)(2Qn —1)—R 
regular polygonal cross section has N sides which divide the (8) 

Y 
Lm 
a 6n/N 
2(m-1) T1/N tp, 41T/N 
P, (X3,Y5) 
- oy PR (X,, Yo) 
— Ue i Pp 1X y ~~ \ \ ; 
m-i\Am-i.'m- 7 2n/Nn 
5» \i oR 1S Ri&.y \ \ \\ 

a 2m7/N V 1 \ ‘ 

>> { Pip (Xn. Sm Q > be (Xo Yo) Xx 

ta Zz \ 4 
% * 1 Vv. 
, 4 \ M ent 
+ ---------- 4 ea SSeS ‘\ S S ‘4 ae | ae ee N 
f 2 Sy PNB PPN-2 / 
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r 4 
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Fig. 1 









Geometry of the section. 








424 JOURNAL OF THE AEROSPACE SCIENCES—MAY 1961 
But it is readily shown that ei 
’ with 
\ N j N > & . > 3 2) I 
VF I Oi(X, y) = Il o,( x, v) v2 log, gn(X, ¥) + - > log, or(X, ¥) + > log, on(X, ¥) | \ li 
n= n=l | n=1 OX n=1 OY n=1 J f ' “ 
. s = {ving 
N N N 9 aD N a o 
cos (2n — 1)x/N |? sin (2n — 1)x/N 71 - 4] 
» Been t= Damns+| EAE +] Sine Nl int 
" { n=1 n=1 PrlX, Y) n=1 Pn. X, V) J f is vil 
. : io | 
\ N N N the | 
2 tos An — K N 
or V2? 1 ¢.(x, vy) = —II ¢,(x, y) 2. i Xx, y) — - cos 2(n K)xr/N l (9 can 
n=1 n=1 fas % n=1 k(X, V)bn(x, ¥) f tern 
The limits of integration must now be determined for S,,, — } Sati 
m = 1, 2, , N — 1, N. For such purpose we most con- metric identity. Hence rm, = R sec [@ — (2m — 1)x/N 
veniently change to polar coordinates. The lower limit with The lower and upper limits of integration with respect to 9 are | 
respect to 7 for each S,, is clearly rm_= 0. The upper limit with seen from Fig. 1 to be respectively @mn, = (2m — 1)x/N and | - 
respect to y for S,, is readily obtained from Eq. (3), letting On, = 2xm/N. Eg. (8) in polar coordinates becomes as 
x = Reos @and y = ¢ sin 6, and employing a common trigono- on(0, r) = ros [9 — (2n — 1)x/N] —R 10 | Gnit 
- — — — —- — — | . 
seric 
. bout 
Hence Eqs. (7) to (11) imply Ty 
n 2mr/N R sec [(0@—(2m—1)r/N] n trial 
9 . ‘ » led 
:F { r II ¢n(0, r)drd9 | ns 
m=1 2(m—1)x/N 0 n=1 } 


— (11) at 8 


Ai = ae | | 
] e2 mr / | ? sec [0—(2m—1)x/N] N 2 N N ‘n “ 
y | { r ; IT (6, nt pa \e “a | ee t drdé the 
1 0 n=1 = d 


2(m—1)4#N k n=1 Hi( 7, O)b,(7, 8) 


Therefore, 





This 
N N 2mnr/N R sec [@—(2m—1)x/N] N | 
2 Il ¢,(x,¥) >> f r Il ¢,(0, r)drd9 — 
. n=1 m=1 (2m—1)xr/N JO coal 9 ' 
a ] 2mx/! ? sec [0— (2m—1)zx/! N ; 2 N N <) torsi 
Y(x, y) F y ; v] j r , : : - (1 
An — R)a/ l 
7 [ . II ¢,(0,r) + :® < oe (7,0) — 2 = > drdd accu 
m=1 J (2m—1)r/N JO ; ws , ee ' n=1 O41, O)dn(r, 0) J of tl 
The torsional rigidity is given by : : ind 
} J ; * Kantorovich, L. V., and Krylov, V. I., Approximate Methods of Highe \t t 
N 2in/N R sec [0—(27—1)r/N] Analysis, Interscience Publishers, Inc., New York, 1958. 
D = 2G z. ry’, 6)drdé, > Stanisi¢, M. M., Free Vibration of a Rectangular Plaie With Damping 0 
J=1 «4 2(93—1)x/N 0 Considered, Quart. Appl. Math., Vol. 12, No. 4, pp. 361-367, January 1955 


(13) 6 Sokolnikoff, I. S., Mathematicel Theory of Elasticity, McGraw-Hill whil 


Book Co., Inc., New York, 1956. 
where (7, 0) is obtained from Eq. (12) by substituting « = 


r cos @and y = ¢ sin 9, and G is the shear modulus. + ms Solv 
The stress tensor is given by 
0 0 Gay,(x, y) The Torsion of Prismatic Bars of Regular 
[rij] = 0 0 —Gay(x,y) | (14) Polygonal Cross Section 
Gay,(x,v) —Gay.(x, y) 0 


where a is the twist per unit length and the subscripts x, y denote F. Ww. Niedenfuhr and A. W. Leissa 
The Associate Professor and Assistant Professor, Respectively, T 


the corresponding derivatives with respect to x and y. A : 
' . ' : Department of Engineering Mechanics, The Ohio State University 


displacement field is given by also, Consultants, North American Aviation, Columbus, Ohio side 
=e August 22, 1960 
[u;] = axs (15) 
aF(x, y) ; i : 
HE DIFFERENTIAL EQUATION! to be solved is Ace 
) ; , 
P(x,y) faf(x, y) of(x, y) ( = ee bar 
where F(x, y) = : = dx — dy¢ (16) V*y = —2G6 (1) 
: Po(20, Yo) \ oy ox f ‘ ‘7 
subject to the boundary condition ies 
and f(x, vy) = W(x, vy) + (1/2)(x? + vy?) (17) 
- * ; y = 0, on B (2 

Using the generalized method here presented, it can be easily . : : e 
' s si cenenen : ara : 3 where y is the torsional stress function, G is the shear modulus, : 
shown that StaniSié’s solution previously developed in the mer ager ite tl atone’ Bis the! ; Som 
» . ; ° ° 6 15 ¥ ; y1s ~¥ 1 : ) > De is the boundary 
literature for particular N appears as a special case of this theory. B is the angle o ne aalew ; ane ee ath diieniecaptnataen : T 

of the cross section, and V? is the Laplacian operator in two , 
° . ° e » ° as t 
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shown in the figure. The centroid of the polygon will coincide 
with the coordinate origin 

The .V-fold symmetry of Eq. (3) makes it possible to satisfy 
the boundary condition at all points of B by analytically satis- 
fying Eq. (2) only along that part of the boundary which lies 
in the range 0 < @ S x/N. This means, for instance, if a circle 


is viewed as a polygon with NV = o, it is only necessary to satisfy 


the boundary condition at the single point r = a,@ = 0. This 
can be done by terminating the series of Eq. (3) after the first 
term. One then obtains 

y = —GBr?/2 +A 








Satisfying the boundary condition at r=a then yields 
A = GBa?* 2 


The resulting stress function for the circular boundary is the 
exact solution to the stated boundary-value problem. 

\n approximation? to the solution for regular polygons with a 
finite number of sides can be obtained as follows: Using the 
series of Eq. (3) truncated to contain k, say, terms, impress the 
boundary condition (2) at k points of the boundary. 
order case to consider is that 


The lowest of the equilateral 


triangle, N = 3. In this case an exact solution can be obtained 
by matching the boundary condition at two points only-——viz., 
at 6 = 0, r = a cos (7/3), and 6 = x/3, r = a, for in this case 
the true torsion function is expressible in the form 

y = — GBr? r 4 — Ac me A r® cos 30 
This problem accordingly holds no further interest since its 


solution in closed form is well known.! 


Higher polygons require a larger number of terms in the 


torsion function to approach the exact solution, but surprisingly 
accurate results can be obtained by again using only two terms 


of the series in Eq. (3). In the general case, therefore, write 


2 


—GBr?/2 + Ao 4 


S o- 


y= Ayr’ cos N@ (4) 


ind set Y = Oat the center of a side and at a vertex of the polygon. 
At the center of a side we obtain the condition 


0 = —(GB/2)\a cos (r/N)]2 + Ao + Ayla cos (r/N)]* (5) 
while at the vertex we obtain 
Q = —Gfa?/2 + A, - Awa” (6) 
Solving Eqs. (5) and (6) for Ao and Ay we find 
_) =e sin? (7/N) 
“= QaN-2 | + cos* (r/N) 
GBa? sin? (2/ NV) 
Ao = - - (7) 


2 1 + cos* (x/N) 

The maximum shearing stress in the bar is at the center of a 
side of the polygon, and is obtained from the formula 

oy 
= (8) 


r=acos r/N 
or @=0 


Timaik 


According to Eqs. (4) and (7) then, the maximum stress in the 
bar is given by 


: , N sin? (2r/N) cos ¥~2(2/N) 
Tmar = —GBa cos(2r/N)] 14 7 : 
2 1 + cos” (2/N) 


(9) 

Some results of calculation with Eq. (9) are shown in Table 1. 
The case N = 6, the regular hexagon, is of some interest here 
as the stress values in the Table show that the recently published 
results obtained by StanisSi¢é, Hauck, and Mathias® are in con- 
The 
those 


siderable error, being of the order of 30 percent too low. 
particular formulation of Galerkin’s method 
authors appears to be excessively tedious as well as inaccurate 
when compared with the present analysis. 


used by 


The reference values of (tmagz/GB8a) given in Table 1 were 
computed by keeping five terms in the series of Eq. (3) and 
satisfying the boundary condition at five points of the boundary 
segment 0 S @ S w/N. 


Convergence studies of this method 


READERS’ 


FORUM 425 
TABLE | 

\ ( tmazr/GBa) gq (tmez/GBG )reterencs © Error 
3 0.7500 0.7500 (exact 0 
} 0.9898 0.9549 (exact } 
5 1.0838 1.0309 5 
6 1. 1230 1.0623 6 
7 1.1388 1.0753 6 
ta) 1.1437 1.0801 6 
9 1.1433 1.0808 6 
10 1.1404 1.0796 5 
exact) 0 


] l 


have shown that these reference values are accurate to within 
less than one tenth of one percent 
The torque carried by the bar can always be obtained from 


the well-known formula 


T=2Sf yA (10) 


where the integration is carried out over the region of the cross 


section, but it is simpler in the case of a regular polygon to us¢e 


an approximation. It is easy to show‘ in this case that if we 
write 
T = KGBa (11) 
where a is the radius of the circumscribing circle, that the func 
tional K is bounded above and below by 
NV 2r 2r N? 2r 
sin 1 + cos Ss ks sin? ( ") (12) 
8 N \ Si \ 


Table 2 shows values of these bounds and accurate reference 
values of K for the first few polygons. From this table we see 
that there is little point in making further calculations for the 
stiffness of the higher polygons, since the difference in the upper 
N27 
The reference value of K for N = 5 was obtained by truncating 


and lower bounds on K is less than seven percent for 


“IN 








= 





Fic. 1. Typical polygonal sector. 
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TABLE 2 





N Kiower Kreference Kupper 

bound bound 
3 0.162 0.1949 (exact) 0.268 
4 0.500 0.5624 (exact) 0.636 
5 0.777 0.8461 0.896 
6 0.976 1.0359 1.076 
7 1.110 1.189 
8 1.205 1.271 
9 1.275 1.327 
10 1.330 1.380 
a w/2 w/2 (exact) 4/2 


the series of Eq. (3) at five terms, matching the boundary condi- 
tion at five points, and then carrying out the integration indicated 
in Eq. (10). The reference value for N = 6 was taken from 
reference 4. Note that the value of K for a hexagon given in 


reference 3 is incorrect. 
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A Note on the Motion of a Spinning Rocket 
With Eccentric Thrust 


A. G. Bennett* 

Instructor, Division of Engineering Sciences, Purdue University, 
Lafayette, Ind. 

September 1, 1960 


= is an approximate analysis of the deviation of the axis of 
a spinning rocket due to a transverse moment arising from 
an eccentric thrust that rotates with the rocket. It is assumed 
that there is symmetry about the spin axis, no aerodynamic force, 
constant mass, and that the deviation of the spin axis is small 
enough to allow a certain linearization in the equations of 
motion. The apex (a point on the spin-axis unit distance from 
the mass center) traces an overall circular path with superim- 
posed lobes (nutations), the details of which depend upon the 
ratio of spin-axis moment of inertia to transverse moment of 
inertia. When this inertia ratio is unity, the gyroscopic terms in 
the equations of motion cancel and the apex moves in an un- 
bounded spiral path. If the transverse moment is suddenly re- 
moved, the apex moves in a circular path tangent to the direction 
of motion at the instant the moment was removed. The diameter 
of this circular path can be as much as twice the average diameter 
of the path with the moment acting. A most interesting result 
is that there is an optimum ratio of inertias which minimizes the 
maximum deviation. 

The equations of motion are set in body-fixed axes x, y, z with 
origin at the mass center and with x along the spin axis of the 
rocket. The angular position of x, y, 2 with respect to nonro- 
tating axes Xo, Yo, 20 is specified by the three Euler angles 3 the de- 
viation of the spin axis, y the procession of the line of nodes, and 
g the angle between the line of nodes and the y axis. Quantities 
A, B, C are the mass moments of inertia about x, y, z, respect- 
ively. The moment M due to eccentric thrust is applied about the 
body-fixed y axis. With this and the assumption of symmetry so 
that B = C, the spin velocity about the x axis is a constant, 0. 
Defining the inertia ratioa = A/B and with w, and w, represent- 


. Presently Summer Faculty Associate, Boeing Airplane Company. 


ing the rotation about the y and z axes, the equations of motion 
are 
wy — (1 — a)Qw, = M/B, @ + (1 — a)Qw, = 0 
With y = z = Oat¢ = Oand assuming J is constant, the solution 
of these equations is 
w, = [M/(1 — a)QB] sin (1 — a)at 
w, = [M/(1 — a)QB) [cos (1 — a)at — 1] 


From the coordinate geometry it is known that 


Q = y cos v+o¢ 
w, = y sind’sing + J cos ¢ 
®: = y sind cos ¢g — Jsin ¢ 


Together with the solution for the w’s above, these equations 
phrase the problem for the determination of the orientation of 
the rocket. It is here, as in other fixed-point motion problems, 
that difficulties appear. The complicated and unmanageable 
nature of the equations reflects the basic fact that finite angles do 
not add as vectors. Since infinitesimal angles do, we are led to 
look for some sort of a linearization. One may not presume y or g 
to be small. They must be allowed to increase indefinitely. We 
first eliminate y and find 
w, = Qtandsin ¢g + [3d cos gy — ¢tan B sin ¢] 


w, = Qtandcos ¢g — [8 sin g + ¢ tan’ cos ¢| 


With the assumption that J is small so that tan 0) = J, the terms 
in the square brackets become perfect differentials. Defining 


new variables a = J sin gand y = J cos ¢, we have 
w = Qa + 7, o, =Qy—-—a 


Taking 3 = ¢ = Oat? = 0, the initial conditions for these equa- 


tions area = y = Oatt = 0. With m = M/AQ? the solution 
may be written 

a/m = {sin(l — a)Qt]/(1 — a) — sin Qt 

y/m = [eos (1 — a)Qt]/(1 — a) — cos Qt — [a/(1 — a)] 


To investigate the solution, it is convenient to transfer back to 
the nonrotating coordinates x», yo, Zo. Omitting the geometric 
details 

a sin Qt 


cos aat a cos Qt sin aQt 


/m = - — 1, v.,,/m = 


vd a 
—a l-—a l-—a l-—a 


vo 
The angles J,, and #,, are the components of 3 about the x) and 
yo axes. For small 3, the angle 3 is approximately equal to the 
linear displacement of the apex. That is, 3,, corresponds to y, of 
the apex and #,, to negative z) of the apex. With this interpreta- 
tion, we see that the path of the apex is an epicycloid with center 
at z)/m = 1. This path may be generated by rolling on a cylinder 
of unit radius, a cylinder of radius 


b =a/(1 — a)for0 <a<1 or 
b= 1/(a — 1) forl<a<2 


Other details of the motion may be obtained from the familiar 
geometry of an epicycloid. Note that the number of lobes (nuta- 
tions) in the overall rotation about the point zo/m = 1 is just the 
reciprocal of b. If we denote the distance from this same point 
to the apex by r/m, we have 

1+a| 


Ymin/m = 1, Tmar/M = 
l-—a 


Note that 7maz/m is unbounded as the inertia ratio approaches 
unity. This is the situation in the nonlinearized solution as well. 
When A = B = C the gyroscopic terms cancel and the motion 
is a simple acceleration under the action of the applied moment. 
The solution for this case shows the path of the apex is an un- 
bounded spiral. This fact leads one to believe that there may 
be an optimum value of inertia ratio for minimum deviation. 
Such is the case if we set the question appropriately. We wiil de- 
termine the minimum of the maximum deviation as a function 
of inertia ratio for a given mass of rocket, given rotational speed, 
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and given transverse moment. The maximum deviation deter- 


mined previously may be written in the form 


27 1 l+a 


m Ail-a 


For a general body of revolution A may be written as a function 
of the inertia ratio alone. We find that the minimum of the 
maximum deviation occurs for an inertia ratio 
a=(3 — Y7)/2 [0.177 
This value of inertia ratio corresponds to //d = 2.78 for a right 
circular cylinder. This optimum of inertia ratio may have useful 
practical application. 
Now consider the rocket with no external moment but with an 
This corresponds to the 
The solution 


initial transverse rotational velocity. 
situation when the moment is suddenly removed 
can be carried out as previously. It is found that the apex traces 
a circular path tangent to the direction of motion when the 


moment as removed. The radius of this motion is 


9 
r/m = sin 1/1 — a)Qt* 
l—a 
where ¢* is the time at which the moment is removed. If the 


moment is removed precisely at the cusp of the epicycloid, 7/m is 
zero and motion of the apex stops. If the moment is removed at 
the center of a lobe of the epicycloid, the radius of the remaining 
motion is 2/(1 — a), twice the average deviation under the in- 
fluence of the constant moment. 

The results of the analysis are self-explanatory and little com- 
ment is necessary. The main point of the analysis is the method 
of linearization which allows the simple solution. The lineariza- 
tion is not objectionable for some problems. In the case of the 
spinning rocket, large deviations eliminate the problem—per- 
haps the rocket. 

In reference 1, an analysis similar to that presented here is 
made. By utilizing automatic computing equipment it was 
possible to obtain particular solutions including the nonlinearities 
of the basic equations such as varying moments of inertia. How- 
ever, the numerical solution does not lend itself to an easy in- 
terpretation of the general nature of the motion. It is quite 
correctly pointed out in this reference that in a rocket the mass 
may vary rapidly and this fact must be taken into account. 
Variable mass causes an obvious difficulty with the idea of opti- 
mum inertia ratio. However, we know that the laws of motion 
apply instantaneously to a system of varying mass (if the forces 
due to expulsion of mass are included). At each instant we have 
just a new initial value problem with the same basic equations 
The nature of the motion is not grossly different. We expect that 
maintaining the inertia ratio in the vicinity of the optimum for 
fixed mass would insure minimum deviation. 


REFERENCE 
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On the Rolling Motion of Low-Aspect-Ratio 
Delta Wings 


Peter J. Mantle* 
Research Engineer, Canadian Armament Research and 
Development Establishment, Valcartier, P.Q., Canada 


September 6, 1960 


OME PRELIMINARY RESULTS obtained from free-flight tests 


on low-aspect-ratio delta wings are presented. These 


results are experimental evidence of the damping in roll (C;,) 


* Now at the Guggenheim Aeronautical Laboratory, California Institute 
of Technology, Pasadena, Calif 


FORUM 2 





























Comparison of theory and test. 


Fic. 1. 


and nonlinear sideslip' (Cigg) obtained in an aeroballistics 


range at supersonic speeds (1.20 < AJ < 2.50) and with aspect 
ratios 0.50, 1.0, and 2.0. These Mach numbers and aspect 
ratios permitted the flight conditions to cover the region bounded 
by near-normal shock and nearly supersonic leading edges 

The application of the aeroballistics range (as a complementary 
tool to the wind tunnel) in dynamic stability analysis of free- 
flight test data has already been discussed in reference 2. 

From some analog simulations and theoretical studies* based 
on slender-body-theory predictions, it was found that the equa- 
tion of rolling motion which best fitted the experimental tra- 


jectories could be written 


~ p t : . p U*Sl a ; 
pina’ ‘ 7 Cip hee ‘ Cla aB = VU (1) 
is J se 8 
where 
> = angle of bank (rad.) 
a = pitch incidence (rad.) 
8 = sideslip angle (rad.) 
p = air density (slugs/cu.ft.) 
U = wing forward velocity (ft./sec.) 
S = wing-planform area (sq.ft.) 
l = reference length (G.M.C.) (it.) 
Izx = moment of inertia in roll (slugs.ft*) 
Ci = rolling moment/!/2pU?Sl 
and 
C oc; ™ o7C; 
i = 70 la, ' 
; ol pl U) od 0ad08 


To indicate the validity of Eq. (1) to represent the rolling 
motion of delta wings, comparison is made between test and 
analog-simulated trajectories for a typical case in Fig. 1. 

It is to be noted that this particular wing experienced roll 
rates equal to the natural frequencies in pitch and yaw,‘ but 
that any expected resonance was not experienced due to the 
high rates of acceleration in roll through these critical values. 
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Fic. 2. Damping in roll of low-aspect-ratio delta wings, 


Replacing Eq. (1) by its equivalent finite-difference form 
and applying the method of least squares to minimize experi- 
mental errors solutions were found both by hand and on an 
IBM 650 digital computer. Five-point differentiation on 
parabolically smoothed data was used to achieve the necessary 
accuracy. In a certain number of cases, however, certain 
mathematical difficulties were incurred by a weighting in the 
matrices of some elements. This weighting depended upon the 
relative magnitudes of the amplitudes in pitch, roll, and yaw. 
These difficulties were discussed in reference 1 and reference 5 
but only the finalized results are presented here. 

The values of damping in roll and the nonlinear sideslip have 
been normalized for generality in Figs. 2 and 3. These curves 
cover the range for a given aspect ratio from normal shock to 
supersonic leading edges. It is seen that only for very low- 
aspect ratios (slender wings) do the theoretical estimates approach 
the present test values 

Ribner and Malvestuto® modified the slender-body-theory 
results to account for the effect of the proximity of the Mach 
cone to the leading edges of the wing; but this effect has been 
seriously underestimated. 

A probable cause of this marked disagreement is the neglect 
in the theory of the effect of the spanwise variation in Mach 
number due to yawing and also, of more importance, the as- 
sumption of the wrong edge condition by ignoring the leading- 
edge vortex flow. 

As the Mach cone approaches the wing leading edge 
(AR VM?—-1— 4.0), both the damping in roll and nonlinear 
sideslip reduce to very small values. This may be explained 
by consideration of the case of just-supersonic leading edges 
(AR V 2 — 1 = 4.0). For the wing in this condition, if it 
sideslips to port, the effect of sweepback will be to induce a nega- 
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Fic. 3. Nonlinear sideslip of low-aspect-ratio delta wings. 


tive roll rate while the now-supersonic starboard leading edge 
with its lower pressure will induce a positive roll rate. These 
two counteracting effects will thus yield either a zero or small 
residual value of C,_,. 
a3 
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Higher Order Theory of Curved Shock* 


Reuben R. Chow and Lu Ting 

Research Assistant and Research Associated Professor, respectively, 
Aerodynamics Laboratory, Polytechnic Institute of Brooklyn, 
Freeport, N.Y. 

September 1, 1960 


me ON THE ASSUMPTION that the Navier-Stokes equations 
hold inside the shock layer, the stretching or boundary- 
layer technique proposed by R. von Mises! and G. Ludford? 
serves for a systematic approach in the analysis of the shock 
structure. The usual Rankine-Hugoniot shock conditions asso- 
ciated with zero shock thickness appears consistently in the 
analysis as the Reynolds number approaches infinity. 

For moderately high Reynolds number, L. I. Sedov, M. P. 
Michailova, and G. G. Chernyi*® proposed a theory for the 
corrections of Rankine-Hugoniot shock conditions taking into 
account the effect of viscosity, heat conduction, and the shock 
curvature behind a curved strong shock. The effect of shock 
structure was neglected in the theory presented. Probstein‘ 
has attempted to improve Sedov’s solution by taking into account 
the effect of shock structure. In his analysis it is necessary to 
introduce an arbitrary shock thickness and approximate the 
velocity profile by a linear distribution. In the present note, it 
will be shown that this arbitrariness and approximation can be 
avoided by the systematic expansion scheme of references 1 
and 2 which has also been developed for the boundary-layer 
theory® and the mixing problem.® 

It is well known that the zero-order shock-structure analysis 
involves only the local slope of the shock surface, and the shock 
conditions remain the same if the shock surface is replaced 
locally by an oblique shock with the same inclination, because 
the shock conditions involve only the local quantities before and 
behind the shock but not their derivatives. In the next-order 
shock condition the first derivatives of the local quantities will 
appear while the normal derivatives can always be related to 
the derivative along the shock. To preserve the same local 
first derivative along the shock, the two-dimensional shock 
surface can be replaced by an equivalent circular shock with 
the same local curvature and slope. It is, therefore, necessary 
to develop the next-order shock relationships for a circular 
shock only. 

Similarly, the next-order shock conditions across an axially 
symmetric shock are locally equivalent to those across a spherical 
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shock. The higher order shock relationships for a spherical ji T 

shock will be reported later. } + ont : 
In order to demonstrate the technique of stretching of co- ba \ \ 

ordinates, the matching conditions, and how to obtain the higher 





order shock relationships, we discuss in detail the circumferential \\'% 
component v. With the radius of the circular shock equal to \ '3s 
unity, and subscripts 1 and 2 designating the quantities before \ 
and behind the shock, the velocity component v can be expressed % 
as 
vy, = 1% = U sind (1) 
vor, 3, €) = ve (7, BD) + evo (r, I) + Ole) (2) 
where ¢ is the reciprocal of the Reynolds number, defined as: 
€ u/p.U. The symbols without subscripts 1 or 2 indicate 
the quantities inside the shock. The coordinate ry normal to caoet 
the shock will be stretched to s by the factor 1/e, i.e., 7 — 1 = es. 
The velocity component inside the shock can be expressed as: 
v(s, J, €-) = vO(s, 3) + ev(s, 3) + Oe) (3) 
The matching conditions on the boundary of the shock layer 
for the zero-order and the first-order solutions are, respectively / 
Hf 
lim v(s, 3) = 7,(1+*, 89) = Usind; lim v™ = 2,!%(1-, 9) / 
oo a Poy ae a 20 40 60 60 @-deg 
(4) Fic. 1. The first-order tangential component of velocity 
lim v(s5, 9) = 0 lim v(s, 39) = 
=> a ° ” . «4s . 
v(1-, 9) + (0/dr)v2(1-, d)s_ (5) are introduced. The matching conditions are equivalent to 
: ’ . é Eqs. (4) and (5). The normal coordinate inside the shock is 
If dv. /dr does not vanish behind the shock, v™(s, 3) cannot stretched so that the derivatives of a quantity with respect to s or 
his Silos. as s > —© while v2 at the shock is given by the # will not change its order of magnitude. After the substitution of 
equation Eq. (3) and its equivalent into the Navier-Stokes equations, 
vO(1-, 9) = lim [v(s, 8) — 5(2/ar)r(1-, 8) terms of like orders of € are collected and equated to zero. From 
' , ea 5, , a the tangential component of the equation of motion, the zero- 
order differential equation is Ov /ds = 0; hencev = U sind 
In Sichel’s analysis? for a weak normal shock v)(s, #) is finite v(1-, 8) from Eq. (4). The zero-order shock condition can 
as s — — © because behind a normal shock v2 vanishes. be obtained directly by integrating the differential equation from 
For other quantities before, behind, and inside the shock s = —» to +. Similarly, the next-order shock condition 
layer, expansion schemes corresponding to Eqs. (1), (2), and (3) can be obtained by integrating the next-order differentia! equa- 
pe Saree = asl teas i tion and applying Eq. (5) without solving for v™. The result is 
7 dv." ai Hare) o +e PS 
Rp ur = R, 4 ( | (pp — px) [ds + —(p — p.) lds 4 
€ oO” J. Lee 0 ov 
nol £? +e ase 
Yo > | [ (p — p2)Os + er = a, t= [= 280) EZ (m4 — uw) (6) 
i Ss 0 g € e _| ov 


The integral terms on the right-hand side of Eq. (6) are associated with the zero-order shock structure. They can be evaluated either 
4 where we can 


analytically or numerically with the help of the zero-order structure solution. A compact form is obtained for P; = 3 
perform the transformation from the physical coordinate, s, to the velocity coordinate u when evaluating the integrals: 


vol) l j y¥+3 = 1 (7 + 1)M 2 cos? J 
= y()* = tan d - - - a 
U 1+, ‘Wy+1 ++1!1 M 2 cos? 3 (y — 1)M.,2 cos? 8 4 2{ 


= aa ) = ( )i | 1 ( : ‘) tn | 
— log, 2)7 sin J + < og, 2 t - = > sin J 
310 +7)? ( W3s\i+7/-" y— 1) M,*%costd + 2]f 


— 
w 


B = (3/4)(ur/p) 


The result is shown in Fig. 1. Following the same procedure, 
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Note on Three-Dimensional Free-Mixing{ 


Martin H. Bloom 


Professor of Aeronautical Engineering, 
Polytechnic Institute of Brooklyn, Freeport, N.Y. 


September 13, 1960 


_— NOTE describes some quasi-two-dimensional properties 
of three-dimensional viscous free-mixing more fully discussed 
in reference 1. 

The following equations of boundary-layer type can be used 
to describe three-dimensional free-mixing flows (the notation 
of reference 2 is used): 


continuity: 


OR 1 re) 
= J [exesRVi) + 
X1 


re) 
1e2R Ve 
ot + €1€2€3 la: oe [eres 2} + 


y) — 
> [eye RV3 7 pos) ¢ =( (1) 
X3 


x; momentum (neglecting body forces): 


OV; a OV; * OV; a pus Ov; 
R + V; + V2 + {Vs + + 
ot €;0X C20X2 R €30X3 


Vi V2 Oe; V2? Oe2 a oP 
€\€2 OX €)C2 OX, a 0X1 


1 (a) e ( e;0( Vi/e1) R 4 
€17€2€3 Ox: _—— ” €30X3 e} ae 


l re) m ei( Vi /e ) : 
e1°@3 (2) 
€17@2€3 OX2 st €20X2 





xX». momentum: 


OV. _ Ove _ ove - pv; \ OV2 
R +Uu—+hr— +(H+ + 
ot €,0X, C20X2 R €30X3 
V2 V3 Oe2 Vi V2 Oe2 Vi? Oe; = oP 4 
€2€3 ONX3 €1€2 Ox, €1€2 OX» a C20X2 


, u ) ( €20( V2, ‘e2) e e€3 O ( V3 R + 
€102763 Ox e1€2 M e Dx: M i ox» m C2U3 T 
1 oO = é& O ( Vi ] ' 
2 ~ | C2°@3ue (3) 
€1€2°€3 OX} €2 Oxe \ e 


X3 momentum: 


P = P, + 0(6) (4) 
energy equation: 
R oH +y, oH Vy oH (m4 pv3;\ OH = 
ot €;0X} €20X2 R ] €30x3 
1 oes 0H  1—Pr dh RuH\ |. 
€1€2€3 OX3 | aa E €30X3 Pr * | ~ Rost f oe 
and in the inviscid flow: 
Vil OVie/e10x1) = —(1/R.)(OP./e:0%x;) (6a) 
( Vie?/e1€2)(0e:/Ox2) = (1/R.)(OP./e20x2) (6b) 


An orthogonal (within the boundary-layer approximations) 
curvilinear coordinate system for this problem can be set up 

t+ This work was done in connection with AFOSR, Contract AF 49(638)- 
217, AFOSR TN 60-986. 
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as follows: A reference streamline S within the viscous core, 
‘center’ is selected. Coordinate surfaces 


preferably near the 
x2 = constant are defined to contain various inviscid streamlines 
at the outer edge of the viscous layer and S as well. Quasi- 
cylindrical surfaces x; = constant enclosing S form the second 
coordinate family, while surfaces x, = constant approximately 
normal to S complete the system. For the special case in which 
the outer streamlines are straight and parallel, the familiar 
cylindrical coordinates can be used (x; = x, x = VU, Xs = 
é: = 1l,es =4r,e3 = 1) 
Consid 


erations identical to those presented in reference 2* 
show that if « 


of outer streamline curvature in the plane tangent to the oute 


(1/e1e2)( Oe; /Ox2) — O, where x is the component 


most surface x; = constant, the momentum equations exhibit 
homogeneity and admit the solution V2 = yor; = 0. This is 
also a consequence of the homogeneity of the boundary conditions 
which in turbulent flow is based on the assumption »7,; = 0) 
in the neighborhood of S, where V2 = V3; = 0. That is, at 
x3 = Oand x3— &, V2 = wy; = 0. Of course, the initial condi- 
tions (at x; = 0) must be consistent with these conditions 


Therefore, initially swirling flows for which V2 # O (for instance), 
are excluded; for e:ample, see reference 5. 

The two-dimensionality will hold within varying degrees of 
approximation when the term [(Vi/Vi.)2 — (R,/R)|«, which 
appears in Eq. (3) is negligible. This term need not be small 
if « * O and the inviscid Mach numbers are high. 

When cross-flow components are not zero but are small com- 
pared to Vj, the approach used by Carrier’ in dealing with flow 
along a corner can be attempted. This involves solving the 
x; momentum equation together with the continuity equation 
which is arbitrarily divided into two equations. This provides 
three equations for the velocity components, the other two 
momentum equations being neglected in a first approximation 

The determination of the absolute orientation of the three- 
dimensional free-mixing region depends upon the establishment 
of a third boundary condition from which the additive functions 
of the cross-flow components V2 and V; can be determined 
This well-known situation has been discussed most recently by 
Ting‘ for the two-dimensional mixing of two semi-infinite uni- 
form streams. Except in special cases the additional boundary 
condition stems from the requirement that the pressure across 
the mixing region is constant to orders higher than zero. Hence, 
the first-order self-induced pressures due to displacement in 
each of the two streams must be the same. This requirement 
provides the condition by which the first-order cross-flow is 
determined absolutely. 

An analogous but more complicated condition holds for three- 
dimensional mixing. In this case, it develops that the two cross 
flow components at the reference axis must be adjusted so as to 
minimize the mean deviation from the mean of the self-induced 
pressure increments at a given streamwise station. When the 
flow is symmetric about one plane, only one velocity component 


can be adjusted to minimize the pressure deviations. For flows 


symmetric about two planes, the determining condition (V2 
V; = O at x; = 0) is evident and higher order solutions need 
not be considered for this purpose. 

It is noted that the coordinate system used to derive the cross 
flow conditions is not quite the same as that previously described, 
which postulates that x; = 0, coincides with a streamline S 
The reference line S, in the latter instance must be an absolute 
one fixed to the inviscid flow. Since Veand V3 are not necessarily 
zero on Sg, the position of S is finally measured with respect 
60: Se, 
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* But with the further restriction that R and yz are constant and eie2es i 


independent of x2 
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Losses in Subsonic and Transonic Axial-Flow 
Compressors at Optimum Conditions 


Gaetano E. Provenzale 
Head, Compressor and Turbine Section, Curtiss-Wright Corp.., 
Wright Aeronautical Div., Wood-Ridge, N.J 


September 23, 1960 


N*s*) has designed and tested a considerable number of 

axial-flow compressor cascades, as well as subsonic and 
transonic single stages. While the losses of the rotor mean 
and hub radial sections and of the whole span of the stator 
blades have been correlated reasonably well by the two-dimen- 
sional cascade loss parameter @ cos a2/2¢ and a loading factor 
called ‘‘diffusion factor,’’ the rotor-tip radial section correlation 
is obviously unsatisfactory. 

Presented and discussed are the results of an analytical 
investigation based on test data at minimum-loss conditions of 
numerous NASA?~!4 and Curtiss-Wright single stages. The 
results, shown in Figs. 1 and 2, indicate an average improvement 
in the correlation of the rotor-tip test data of about 60 percent 
over the NASA correlation.! Figs. 1 and 2 permit a better 
prediction of the relative total-pressure-loss coefficient at mini- 
mum-loss angle of incidence for flow through subsonic and 
supersonic axial-compressor rotor tips having 65-(A10) or double 
circular-are airfoils. The main conclusion of this paper is 
that higher values of solidity than those employed in most of 
the existing compressor rotors are necessary to obtain higher 
rotor-tip efficiency, especially at supersonic Mach numbers. 

The loss parameter @ cos a2/2¢0 was developed theoretically by 
Lieblein and Roudebush and was shown to correlate well the 
two-dimensional cascade test data as a function of diffusion 
factor.!5 The test data for the rotor mean and hub sections, 
and the stator tip, mean and hub sections were correlated 
reasonably well by the same loss parameter and diffusion factor,! 
probably because the flow at these blade sections is closer to a 
two-dimensional type of flow, while the rotor-tip section ex 
periences various flow phenomena of a predominantly three- 
dimensional type. Test results of some NASA*® and WAD 
transonic rotors, and later analyses'*~'8 had indicated that the 
scatter of the tip loss parameter could be related to local shock 
losses arising from increased blade-surface Mach numbers 

The writer initiated an investigation in order to rationalize 
the tip flow phenomena and improve the rotor-blade design and 
loss prediction in the tip region. The test data of subsonic 
flow were separated from the test data of supersonic flow. The 
subsonic Mach number range available (0.3 to 0.97) covers 
typical values of hub and mean section Mach numbers. This 
separation of the test data should have reduced the scatter of 
the tip loss parameter, but it did not. Fig. 3 shows the still 
All this indicated that the shock 
The idea of a 


very high degree of scatter 
losses were not the main reason for the scatter. 
continuous tip-loss function of Mach number was developed, 
and the loss coefficient became a power function of solidity in 
terms of an exponent » function of Mach number. Values of 
n = 1.0 for M = 1.0 and n = 2.0 for M = 1.1 correlate the 
test data satisfactorily as shown in Figs. 1 and 2. The nature 
of the m function should be determined, as more test data become 
available at Mach numbers greater than 1.1. 

The various flow phenomena present at the tip, such as 
thickening of boundary layer due to radial secondary flows, 
radial clearance, etc., may explain why the tip behaves differently 


FORUM 431 









































10 — 
BL ADE My TIP 7 
REFERENCE SYMBOL SHAPE RANGE SOLIDIT v 
2 c oca 13 | 
3 Oca 64-97 625 
a 4 OcA 65-68 63-6 
os 5 < OCA 65- 65 ; 
6 Vv oca 68- 97 66- 68 a 
7 5 oca 64- 96 e) a 
8 e OCA 4 32 4 
9 r Oca 36- 70 906 | 
« 0 + 65S 39- 72 608 
S o6 65S 50-9) 2 | 
= 12 ° 65S 30- 76 59- 6 
= 13 a 65S 49-75 69-94 
a ia + 655 ?? 28. . 
§ C-wCA6 ‘ Oca 3 984 
a c cas , oca =] 492 
2 04 aces I ° 
° ROTOR TIP (1O% of BLADE HEIGHT 
ad FROM COMPRESSOR TIP a 
My, €10, 9210 
: - o% 
4 
+ 
o2 ; 
& 
- a 4 
Oo sf | 
o i a J 
fe) 2 3 4 6 ? 
DIFFUSION FACTOR, D 
Fic. 1. Total-pressure-loss parameter (@o” cos a2)/2 vs. diffusion 
factor D. 
BLADE 
REFERENCE SYMBOL SHAPE 
08 2 c D.C.A - 
3 DCA 
3 4 OCA 
5 ° OCA 
6 v DCA 
7 rs DCA 
06 }— —___—_____— 
a 
rr 
3 
we - — + —_ + 
< ROTOR TIP (lO% of BLADE HEIGHT 
« FROM COMPRESSOR TIP) 
z My, 71.1, 9*2.0 
04 —! 
” 
” 
°o 
02 + , 
rl 
() 
() 1 2 3 4 5 5 
DIFFUSION FACTOR, D 
> ») hy P > ‘ = n ‘ 4 
Fic.2. Total-pressure-loss parameter (Ge" cos a2)/2 vs. diffusion 
factor D 
BLADE Mr, Te 
REFERENCE SYMBOL SHAPE RANGE SOLIDITY 
: 2 fe OCA 1.32 
3 OCA 84-97 825 
~ OCA 65-88 63-60 
5 OCA 65- 65 ie} 
6 Vv ocaA 68- 97 66- 66 
7 oca 64- 96 1o¢ 
12 8 ® OCA Be 32 
9 Oca 36-70 906 ¢ 
° e 65S 39-7 608 
65S 50- 92 2 
12 * 65S 30-76 59- 60 
13 a 65s 49-75 69- 94 
re) 14 65s 77 126. 
C-wC a6 c oca +] 964 
C-wcas > oca 9 493 
+ a 
a - a 
“~ o8}—— + "] t | 
HA ROTOR TIP (0% of BLADE HEIGHT 
= FROM COMPRESSOR T cee a 
a “,=10 
a *~ 
S$ 06}+— 4- { 
w) e 4 
& 
04 
e 
02 e 
o ‘ | 
8, & F | 
as | 
° = j 
oO 2 3 4 5 _ 4 
DIFFUSION FACTOR, 0 
Fic. 3. Total-pressure-loss-parameter (@ cos a.)/2¢ vs. diffusion 


factor D 


than the mean and hub sections at subsonic flow. High values 
of solidity apparently reduce the negative effects of these flow 
phenomena mostly by providing a better flow guidance. Limited 
available test data indicate that the same effect of Mach number 
on loss coefficient exists over the whole blade at supersonic 


speeds. Wright Aeronautical is engaged in a program of high- 
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work, high-solidity single-stage tests to determine the upper 
limit of rotor-tip solidity above which the frictional losses become 
predominant. 

Concluding, the writer hopes to provide an improved method 
for predicting rotor-tip losses and designing better rotors, and 
also to provide material that would stimulate further the interest 
in the development of a general theory of compressor aero- 


dynamics. 
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— suBJECT of matric structural analysis has been treated 
in two recently published papers in the JOURNAL.'’? The 


MAY 1961 


authors of these papers have made a number of statements about 
the inversion of certain large matrices. It is the purpose of 
this note to bring to the attention of the reader certain facts 
that shed new light on this important problem. It is shown 
here that the situation is not as hopeless as the above-mentioned 
authors intimate. 

In general, the time necessary to invert a large, sparse matrix 
of the type encountered in the method of matric structural 
analysis advocated by the writer* 7 is a function of several things , 
as follows in order of importance: 

(1) The machine, method, and program used. 

(2) The number of redundants in the system, say R. 

(3) The number of nonzero elements present in the matrix, 
say N. 

(4) The order of the matrix, say 1. 

Notice first that the size of the matrix itself is the least im- 
portant item. This premise may be surprising. The simple 
example which follows may illustrate this point. Suppose we 
have two systems of equations in which R and WN are the same 
but ./ is different. Then, actually the inversion time can be 
longer for the system with the smaller 1/, contrary to customary 
thinking. It appears from such reasoning that N/M may be 
a parameter of significance. 

Notice next that the type of machine available is the prime 
factor. Naturally, no problem would exist if we had at our dis- 
posal a very fast machine with very large, high-speed, storage 
space. Present-day machines are limited in high-speed storage 
space. Therefore, one has to resort to auxiliary lower speed 
storage that affects inversion time. The method and program 
we use is important. It is believed that the method presented 
in references 8 and 9 can be made highly efficient because it 
isolates very nicely the redundant part of the solution. The 
size of this subsystem R is an important factor in inversion 
time. Really, the degree of difficulty of an inversion is closely 
related to the size of the redundant matrix. This matrix size 
is the number of equations that cannot be put into a pretri- 
angularized form. It represents the number of redundants in 
the system in the structural sense. Notice that this number 
is determined automatically in the program.*? 

Naturally the smaller the number of nonzero elements in the 
matrix, the less shuffling that has to be done. 

In conclusion, we say that the time for inversion may be a 


complicated function, 
= 7(R, N, M) 


and one cannot draw any conclusions on the basis of JJ alone. 
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